[论文解读] Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

Jacques Masson, G. Chabrier+3 more|Figshare|Sep 18, 2015

Astrophysics and Star Formation Studies参考文献 51被引用 79

一句话总结

本文通过非理想磁流体动力学（MHD）模拟研究了低质量恒星形成过程中的非理想磁扩散效应，结果表明，该效应通过形成磁扩散屏障抑制了磁场通量堆积，从而实现了开普勒盘的形成——而理想MHD则因过度制动角动量而阻止了盘的形成。引入非理想磁扩散后，可形成物理上合理的盘结构，其质量为0.02–0.2 M⊙，半径为20–80 au，具体取决于初始磁化程度。

ABSTRACT

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

研究动机与目标

• 解决理想MHD模拟原恒星核坍缩过程中出现的非物理解磁通量堆积与过度磁制动问题。
• 研究非理想磁扩散在第一 Larson核形成过程中对角动量输运与磁通量重分布的调控作用。
• 确定非理想MHD是否能产生与Class-0和Class-I原恒星观测结果一致的旋转支持盘。
• 评估初始磁场方向与磁化程度对盘形成与稳定性的影响。
• 评估数值方法在模拟磁扩散过程中的作用，避免数值扩散带来的非物理伪影。

提出的方法

• 采用RAMSES自适应网格加密代码与非理想MHD求解器，模拟1 M⊙预恒星核的坍缩过程。
• 通过简化电荷、中性与尘埃粒子的化学网络计算电阻率，引入非理想磁扩散效应。
• 在核附近设置磁扩散屏障（解耦阶段），防止磁场放大超过约0.1 G。
• 采用多种黎曼求解器（HLL、HLLD、Lax-Friedrich）测试数值扩散效应对盘形成与稳定性的影响。
• 通过径向速度剖面定义盘的性质，并与开普勒预期进行比较，同时考虑非开普勒偏差。
• 在不同初始磁化度（μ = 2 至 μ = 5）和磁场方向与旋转轴的夹角下开展模拟。

实验结果

研究问题

• RQ1 在坍缩的预恒星核中，非理想磁扩散与理想MHD相比，如何改变磁通量重分布？
• RQ2 非理想MHD模拟能否产生与Class-0和Class-I原恒星观测结果一致的旋转支持盘？
• RQ3 磁扩散屏障在调控磁场放大与防止过度制动中起什么作用？
• RQ4 数值求解器选择（如HLL与HLLD）如何影响第一Larson核的演化与盘的形成？
• RQ5 初始磁场方向与磁化程度在非理想MHD模拟中对盘的质量与尺寸影响有多大？

主要发现

• 非理想磁扩散在第一Larson核附近形成磁扩散屏障，将磁场放大限制在约0.1 G以内，有效防止了磁通量堆积。
• 盘的质量范围为~0.02 M⊙（μ = 2）至~0.2 M⊙（μ = 5），半径范围为20–30 au至40–80 au，具体取决于初始磁化程度。
• 所有非理想MHD模拟均恢复了开普勒速度剖面，而理想MHD则抑制了盘的形成。
• 理想MHD中的磁化喷流在非理想MHD中因磁制动减弱而变弱或消失。
• 在非理想MHD中，磁场方向与旋转轴的夹角对盘形成影响极小，与理想MHD中显著不同，证实了物理机制的主导作用而非数值效应。
• 理想MHD模拟表现出非物理行为，如喷流中的反向旋转、磁互换不稳定性以及数值磁通量重分布，而这些在非理想MHD中均未出现。