[論文レビュー] Turbulence Can Persist in the Inner Regions of Weakly-Ionized Planet Forming Disks

David G. Rea, Jacob B. Simon|arXiv (Cornell University)|Jan 19, 2026

Astrophysics and Star Formation Studies被引用数 0

ひとこと要約

論文は3D非理想MHDシミュレーションを用いて、傍電解拡散がAmbipolar Diffusion Shear Instability (ADSI)を維持し得ること、さらにHall効果と組み合わせて磁気応力を高め、弱結離散子星形成系の薄いディスクの中平面乱流を生じさせることを示す。約5 AU付近の弱電離プロトラ planetary ディスクにおける乱流を示す。

ABSTRACT

Identifying the mechanisms responsible for angular momentum transport in protoplanetary disks, and the extent to which those mechanisms produce turbulence, is a crucial problem in understanding planet formation. The bulk of the gas in protoplanetary disks is weakly ionized, which leads to the emergence of three non-ideal effects, Ohmic diffusion, ambipolar diffusion, and the Hall effect. These low-ionization processes can in some cases suppress turbulence driven by the magnetorotational instability (MRI). However, it has recently been shown that these non-ideal terms can also affect the dynamics of the gas in fundamentally different ways than simple diffusion. In order to further study the role of low-ionization on disk gas dynamics, we carry out a 3D local shearing box simulation with both Ohmic diffusion and ambipolar diffusion and an additional simulation with the Hall effect included. The strength of each non-ideal term, when present, is representative of gas at a radius of 5 AU in a realistic protoplanetary disk. We find the Hall effect increases the saturation strength of the magnetic field, but does not necessarily drive turbulence, consistent with previous work. However, interactions between ambipolar diffusion and the Keplerian shear lead to the ambipolar diffusion shear instability (ADSI), which can drive the initial growth, not damping, of magnetic perturbations. To our knowledge, this is the first work that explicitly demonstrates the viability of the ADSI in the non-linear regime within protoplanetary disks. At later times in the disk, the MRI (reduced in strength by ambipolar-diffusion), may also be present in regions of weak magnetic field between strong concentrations of vertical magnetic flux and sustain turbulence locally in protoplanetary disks.

研究の動機と目的

• 弱電離プロトプラネタリディスクにおける角運動量輸送と乱流の理解を動機づける。
• ディスク乱流に対する非理想MHD効果（Ohm拡散、傍電拡散、Hall効果）の役割を調査する。
• ADSIおよび関連不安定性が内側のディスクで乱流を駆動または維持できるかを評価する。
• 局所のShearing-box枠組みにおける磁場幾何が応力と乱流に与える影響を特徴づける。

提案手法

• 5 AUのプロトプラネタリディスクのパッチをAthenaで3D局所シェアリング-boxシミュレーションを実施する。
• 垂直層化を伴う等温ガスと垂直磁場を用い、2つのケースを探索する：Hall効果あり（中平面で OA H a = 0.054）、Hall効果なし（OA）。
• Ohm拡散η_O、傍電拡散η_A、Hall効果η_HをElsasser数 Λ、Ha、Amで特徴づける。
• 連続方程式・運動方程式・誘導方程式を非理想項η_O J、η_H J×B/B、η_A (J×B)×B / B^2を含めて解く。
• 1000 Ω^{-1}に渡ってMaxwell応力・Reynolds応力・乱流速度・磁場幾何を解析する。
• 拡散を含む分散関係（Desch 2004）とMRI/ADSIの寄与を評価して線形成長率と比較する。

Figure 1: Maxwell (top row) and Reynolds (bottom row) stresses for OA and OAH. The left column shows the box-averaged values as a function of time, and the right column shows the time- and $xy$ -averaged values as a function of height. The Reynolds stress is highly variable; for easier viewing, a so

実験結果

リサーチクエスチョン

• RQ1 内側ディスクでAmbipolar Diffusionだけ、またはOhmic拡散と組み合わせた場合にADSIを駆動できるか？
• RQ2 Hall効果は磁場強度や応力を高めるか、内側ディスクで乱流を必ず駆動するか？
• RQ3 シミュレーションに現れる不安定性は何で、それらはKeplerianせん断とどのように相互作用して乱流を sustainするか？
• RQ4 大規模な磁場幾何（even-z vs odd-z）は乱流と角運動量輸送にどう影響するか？
• RQ5 弱いイオン化領域の内側ディスクで観測される乱流を長期的に維持できる程度はどの程度か？

主な発見

Name	ζ_CR,0/s^-1	β0	Ha 0	⟨M_xy⟩mid	⟨M'_{xy}⟩mid	⟨R_xy⟩mid	⟨δv/c_s⟩mid	⟨δv/c_s⟩mid*	σ/Ω
OAH	-17	4	5.4e-2	1.11e-3	3.88e-4	3.95e-4	7.77e-2	5.15e-2	0.66
OA	-17	4e4	Infinity	3.47e-4	1.84e-4	1.13e-4	4.94e-2	4.19e-2	0.66

• Hall効果はディスク中平面の総Maxwell応力を約3倍に増加させるが、主に径方向磁場が強化されているためである。
• Hall効果は磁場強度を高めるが、必ずしも持続的乱流を駆動するとは限らず、中平面の乱流はHall有無のいずれでも存続する。
• 傍電拡散はKeplerianせん断と相互作用してADSIを生み出し、MRIが減衰しても磁気摂動の初期成長を駆動できる。
• ADSIは原始的な非線形領域で作用し、特定の領域で持続的な磁気摂動と乱流をもたらす。
• 内平面の乱れ速度は ~5×10^-2 c_s に達し、|z| ≳ 5 H で c_s に向かって上昇し、内側ディスク領域での乱流が持続することを示す。
• 初期磁場増幅の成長率σ ≈ 0.66 Ω で、選択された拡散制約下でMRIとADSIを組み合わせた線形解析と一致する。
• 大規模磁場幾何（even-z vs odd-z）は乱流振幅と電流シートダイナミクスと相関し、長期的な応力変動に影響を与える。

Figure 2: Spacetime diagrams of the azimuthal field $\langle B_{y}\rangle_{xy}$ for each simulation. Both simulations exhibit multiple transient current sheets up to $200~\Omega^{-1}$ ( $\sim~30$ orbits). Thus, we only consider field geometry for $t>200~\Omega^{-1}$ . The field geometry of OAH consi