Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Generalized Nonlinear Ohms Law: How a Strong Electric Field Influences Non-ideal MHD Effects in Dusty Protoplanetary Disks
56
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The magnetohydrodynamics (MHD) of protoplanetary disks are strongly subject to the non-ideal MHD effects arising from the low ionization fraction of the disk gas. A strong electric field induced by gas motions can heat ionized gas particles and can thereby affect the ionization balance in the disks. Our previous studies revealed that in dusty protoplanetary disks, the Ohmic conductivity decreases with increasing electric field strength until the electrical breakdown of the disk gas occurs. In this study, we extend our previous work to more general cases where both electric and magnetic fields affect the motion of plasma particles, allowing us to study the impacts of plasma heating on all non-ideal MHD effects: Ohmic, Hall, and ambipolar diffusion. We find that the upper limit on the electric current we previously derived applies even in the presence of magnetic fields. Although the Hall and ambipolar resistivities can either increase or decrease with electric field strength depending on the abundance of charged dust grains, the Ohmic resistivity always increases with electric field strength. An order-of-magnitude estimate suggests that a large-scale electric current generated by gas motions in the inner part of protoplanetary disks could exceed the upper limit. This implies that MHD motions of the inner disk, such as the motion driven by the Hall-shear instability, could either get suppressed or trigger electrical breakdown (lightning discharge). This may have important implications for gas accretion and chondrule formation in the inner part of protoplanetary disks.
rate research
Read More
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.
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.
Small solids embedded in gaseous protoplanetary disks are subject to strong dust-gas friction. Consequently, tightly-coupled dust particles almost follow the gas flow. This near conservation of dust-to-gas ratio along streamlines is analogous to the near conservation of entropy along flows of (dust-free) gas with weak heating and cooling. We develop this thermodynamic analogy into a framework to study dusty gas dynamics in protoplanetary disks. We show that an isothermal dusty gas behaves like an adiabatic pure gas; and that finite dust-gas coupling may be regarded as an effective heating/cooling. We exploit this correspondence to deduce that 1) perfectly coupled, thin dust layers cannot cause axisymmetric instabilities; 2) radial dust edges are unstable if the dust is vertically well-mixed; 3) the streaming instability necessarily involves a gas pressure response that lags behind dust density; 4) dust-loading introduces buoyancy forces that generally stabilizes the vertical shear instability associated with global radial temperature gradients. We also discuss dusty analogs of other hydrodynamic processes (e.g. Rossby wave instability, convective overstability, and zombie vortices), and how to simulate dusty protoplanetary disks with minor tweaks to existing codes for pure gas dynamics.
Non-thermal acceleration of particles in magnetohydrodynamic (MHD) turbulence plays a central role in a wide variety of astrophysical sites. This physics is addressed here in the context of a strong turbulence, composed of coherent structures rather than waves, beyond the realm of quasilinear theory. The present description tracks the momentum of the particle through a sequence of frames in which the electric field vanishes, in the spirit of the original Fermi scenario. It connects the sources of energy gain (or loss) to the gradients of the velocity of the magnetic field lines, in particular the acceleration and the shear of their velocity flow projected along the field line direction, as well as their compression in the transverse plane. Those velocity gradients are subject to strong intermittency: they are spatially localized and their strengths obey powerlaw distributions, as demonstrated through direct measurements in the incompressible MHD simulation of the Johns Hopkins University database. This intermittency impacts the acceleration process in a significant way, which opens up prospects for a rich phenomenology. In particular, the momentum distribution, which is here captured through an analytical random walk model, displays extended powerlaw tails with soft-to-hard evolution in time, in general agreement with recent kinetic numerical simulations. Extensions to this description and possible avenues of exploration are discussed.
Secular gravitational instability (GI) is one of the promising mechanisms for creating annular substructures and planetesimals in protoplanetary disks. We perform numerical simulations of the secular GI in a radially extended disk with inward drifting dust grains. The results show that, even in the presence of the dust diffusion, the dust rings form via the secular GI while the dust grains are moving inward, and the dust surface density increases by a factor of ten. Once the secular GI develops into a nonlinear regime, the total mass of the resultant rings can be a significant fraction of the dust disk mass. In this way, a large amount of drifting dust grains can be collected in the dusty rings and stored for planetesimal formation. In contrast to the emergence of remarkable dust substructures, the secular GI does not create significant gas substructures. This result indicates that observations of a gas density profile near the disk midplane enable us to distinguish the mechanisms for creating the annular substructures in the observed disks. The resultant rings start decaying once they enter the inner region stable to the secular GI. Since the ring-gap contrast smoothly decreases, it seems possible that the rings are observed even in the stable region. We also discuss the likely outcome of the non-linear growth and indicate the possibility that a significantly developed region of the secular GI may appear as a gap-like substructure in dust continuum emission since dust growth into larger solid bodies and planetesimal formation reduce the total emissivity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
