Ambipolar diffusion in low-mass star formation

Abstract

Angular momentum transport and the formation of rotationally supported structures are major issues in our understanding of protostellar core formation. Whereas purely hydrodynamical simulations lead to large Keplerian disks, ideal magnetohydrodynamics (MHD) models yield the opposite result, with essentially no disk formation. This stems from the flux-freezing condition in ideal MHD, which leads to strong magnetic braking. In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution by including resistive terms in the MHD equations. We focus more particularly on the effect of ambipolar diffusion on the properties of the first Larson core and its surrounding structure, exploring various initial magnetisations and magnetic field versus rotation axis orientations of a 1 M collapsing prestellar dense core. We used the non-ideal magnetohydrodynamics version of the adaptive mesh refinement code RAMSES to carry out these calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral, and grain species. Including ambipolar diffusion leads to the formation of a magnetic diffusion barrier (also known as the decoupling stage) in the vicinity of the core, which prevents accumulation of magnetic flux in and around the core and amplification of the field above 0.1 G. The mass and radius of the first Larson core, however, remain similar between ideal and non-ideal MHD models. This diffusion plateau, preventing further amplification of the field and reorganising the field topology, has crucial consequences for magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened or even disappear when using non-ideal MHD. In contrast to ideal MHD calculations, misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given magnetisation, showing that the physical dissipation processes truly dominate numerical diffusion. We demonstrate severe limits of the ideal MHD formalism; it yields unphysical behaviours in the long-term evolution of the system. This includes counter-rotation inside the outflow or magnetic tower, interchange instabilities, and flux redistribution triggered by numerical diffusion. These effects are not observed in non-ideal MHD. Disks with Keplerian velocity profiles are found to form around the protostar in all our non-ideal MHD simulations, with a final mass and size that strongly depend on the initial magnetisation. This ranges from a few 10 −2 M and ∼20−30 au for the most magnetised case (µ = 2) to ∼2 × 10 −1 M and ∼40−80 au for a lower magnetisation (µ = 5). In all cases, these disks remain significantly smaller than disks found in pure hydrodynamical simulations. Ambipolar diffusion thus bears a crucial impact on the regulation of magnetic flux and angular momentum transport during the collapse of a prestellar core and the...

Magnetically Self-Regulated Formation of Early Protoplanetary Disks

“…The first type of simulations are identical and/or similar to the ones performed in Masson et al (2016). They have an initial core mass of 1 M , a uniform density profile and a uniform magnetic field with a mass-to-flux over critical mass-to-flux ratio of 2 or 5.…”

Apparent Disk-mass Reduction and Planetisimal Formation in GravitationallyUnstable Disks in Class 0/I Young Stellar Objects

“…Even when the magnetic field in the cloud core is relatively strong (e.g., when its mass-to-flux ratio normalized by the critical mass-to-flux ratio is on order of unity), a gravitationally unstable disk can form during the Class 0/I YSO phase. Indeed, recent three-dimensional non-ideal magnetohydrodynamics (MHD) simulations suggest that even with a relatively strong magnetic field, a circumstellar disk forms immediately after protostar formation Tomida et al 2015;Tsukamoto et al 2015aTsukamoto et al , 2015bMasson et al 2016;Wurster et al 2016) and becomes gravitationally unstable Tsukamoto et al 2015aTsukamoto et al , 2015b; for a review of disk formation in magnetized cloud cores, see Tsukamoto 2016). In particular, Machida et al (2011) investigated the long-term evolution of circumstellar disks (until 10 5 years after protostar formation) and showed that gravitationally unstable disks can form even in strongly magnetized cloud cores.…”

Kinematics of a Young Low-mass Star-forming Core: Understanding the Evolutionary State of the First-core Candidate L1451-mm

“…and resolution. The development of a disk might start right after protostar formation (Inutsuka et al 2010;Machida & Matsumoto 2011;Tsukamoto et al 2015), or a rotationally supported structure corresponding to the first core, or surrounding it, might exist before protostar formation (Inutsuka et al 2010;Machida & Matsumoto 2011;Tsukamoto et al 2015;Joos et al 2012;Tsukamoto & Machida 2013;Masson et al 2016). Observations that can test these predictions are needed in order to better understand the formation of protostars and their evolution.…”