Efficient Numerical Treatment of Ambipolar and Hall Drift as Hyperbolic System

Partially ionized plasmas, such as the solar chromosphere, require a generalized Ohm’s law including the effects of ambipolar and Hall drift. While both describe transport processes that arise from the multifluid equations and are therefore of hyperbolic nature, they are often incorporated in models as a diffusive, i.e., parabolic process. While the formulation as such is easy to include in standard MHD models, the resulting diffusive time-step constraints do require often a computationally more expensive implicit treatment or super-time-stepping approaches. In this paper we discuss an implementation that retains the hyperbolic nature and allows for an explicit integration with small computational overhead. In the case of ambipolar drift, this formulation arises naturally by simply retaining a time derivative of the drift velocity that is typically omitted. This alone leads to time-step constraints that are comparable to the native MHD time-step constraint for a solar setup including the region from photosphere to lower solar corona. We discuss an accelerated treatment that can further reduce time-step constraints if necessary. In the case of Hall drift we propose a hyperbolic formulation that is numerically similar to that for the ambipolar drift and we show that the combination of both can be applied to simulations of the solar chromosphere at minimal computational expense.

Constraint on ion–neutral drift velocity in the Class 0 protostar B335 from ALMA observations

Aims. Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal magnetohydrodynamics (MHD) effects proposed to enable the formation of large-scale Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare here gas kinematics traced by ionized and neutral molecular lines and discuss the implication on ambipolar diffusion. Methods. We analyzed the data of the H13CO+ (3–2) and C18O (2–1) emission in the Class 0 protostar B335 obtained with our ALMA observations. We constructed kinematical models to fit the velocity structures observed in the H13CO+ and C18O emission and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. Results. A central compact (~1′′–2′′) component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. With our kinematical models, the infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85 ± 0.2 km s−1 at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km s−1 at a radius of 100 au in B335. Conclusions. The Hall parameter for H13CO+ is estimated to be ≫1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between the ionized and neutral gas could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. This result is consistent with the expectation from the MHD simulations with a typical ambipolar diffusivity and those without ambipolar diffusion. On the other hand, the high ambipolar drift velocity of 0.5–1.0 km s−1 on a 100 au scale predicted in the MHD simulations with an enhanced ambipolar diffusivity by removing small dust grains, where the minimum grain size is 0.1 μm, is not detected in our observations. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0. ′′1) are needed to examine this possibility and ambipolar diffusion on a smaller scale.

A theorem of existence for the equations of magnetohydrodynamics of partially ionized plasmas

The equations of magnetohydrodynamics of partially ionized plasmas have been known for a long time, but rarely studied. Instead, several simplifications have been applied to different physical models, ranging from magnetic reconnection to ambipolar drift. The original system relates the electric field to kinetic magnitudes by means of a non-local law, so that the equations describing it involve partial differentials as well as functional operators. We prove a theorem of existence and uniqueness of solutions for a finite time by means of a fixed point argument in an appropriate functional setting.