[論文レビュー] Gravitoturbulent dynamo in global simulations of gaseous disks

William Béthune, Henrik N. Latter|arXiv (Cornell University)|Jun 8, 2022

Astrophysics and Star Formation Studies参考文献 82被引用数 10

ひとこと要約

本研究では、重力乱流駆動型乱流が弱い磁場を大規模なダイナモ機構を介して増幅できる可能性を示す、全地球的3次元磁気流体力学（MHD）シミュレーションを提示する。局所的シミュレーションとは異なり、このダイナモは全地球的かつ固定された増幅度で作用し、プラズマベータβ ≈1で飽和に達し、磁場レイノルズ数Rm ≲10の低値で最適な増幅を示す。これは、若い原始惑星系や活動銀河核のような質量の大きなディスクで、強力に磁化された安定した状態が形成されることを示唆している。

ABSTRACT

The turbulence driven by gravitational instabilities (GIs) can amplify magnetic fields in massive gaseous disks. This GI dynamo may appear in young circumstellar disks, whose weak ionization challenges other amplification routes, as well as in active galactic nuclei. Although regarded as a large-scale dynamo, only local simulations have so far described its kinematic regime. We study the GI dynamo in global magnetohydrodynamic (MHD) models of accretion disks, focusing on its kinematic phase. We perform resistive MHD simulations with the Pluto code for different radiative cooling times and electrical resistivities. A weak magnetic field seeds the dynamo, and we adopt mean-field and heuristic models to capture its essence. We recover the same induction process leading to magnetic field amplification as previously identified in local simulations. The dynamo is, however, global in nature, connecting distant annuli of the disk via a large-scale dynamo mode of a fixed growth rate. This large-scale amplification can be described by a mean-field model that does not rely on conventional alpha-Omega effects. When varying the disk parameters we find an optimal resistivity that facilitates magnetic amplification, whose magnetic Reynolds number, Rm < 10, is substantially smaller than in local simulations. Unlike local simulations, we find an optimal cooling rate and the existence of global oscillating dynamo modes. The nonlinear saturation of the dynamo puts the disk in a strongly magnetized turbulent state on the margins of the effective range of GI. In our simulations, the accretion power eventually exceeds the threshold required by local thermal balance against cooling, leaving the long-term nonlinear outcome of the GI dynamo uncertain.

研究の動機と目的

• ガスディスクの全地球的3次元MHDシミュレーションにおける重力乱流ダイナモの運動学的段階を調査すること。
• 局所的（せん断ボックス）シミュレーションで同定されたダイナモ機構が、現実的なディスク幾何学において全地球的に作用するかを特定すること。
• 抵抗率、冷却時間、磁場レイノルズ数が磁場増幅の促進および制御に果たす役割を評価すること。
• ダイナモの非線形的飽和状態と、ディスク安定性および降着に与える影響を検討すること。

提案手法

• 有限な半径領域を持つ全地球的球座標および円柱座標で、Plutoコードを用いて抵抗性MHDシミュレーションを実行する。
• 初期に弱いトロイダル磁場を導入し、オーム抵抗を導入することでMRIを抑制し、GIダイナモ領域を明確に分離する。
• 重み付き方位平均を適用して、ダイナモを局所的および全地球的成分に分解する。
• ダイナモの本質的物理を捉え、全地球的挙動を検証するために、平均場およびヒューリスティックモデルを用いる。
• 放射冷却時間と電気的抵抗率を変化させ、パラメータ依存性を調査する。
• 磁場エネルギーの拡散、流れ構造、乱流統計を分析し、ダイナモモードおよび飽和状態を同定する。

実験結果

リサーチクエスチョン

• RQ1 局所的シミュレーションで同定された重力乱流ダイナモが、3次元ディスク幾何学において全地球的プロセスとして作用するか？
• RQ2 抵抗率と冷却速度が、全地球的シミュレーションにおける大規模磁場増幅を可能にする役割は何か？
• RQ3 全地球的ダイナモは、従来のα-Ω機構に依存しない平均場モデルで記述可能か？
• RQ4 GIダイナモの非線形的飽和状態の性質は何か？ また、ディスク安定性および降着に与える影響は？
• RQ5 このダイナモは、磁場エネルギーが支配的となる強磁化状態（β ≈1）をもたらすか？ その場合、ディスクのエネルギー収支や熱的バランスに影響を与えるか？

主な発見

• GIダイナモは、固定された増幅度を示す全地球的プロセスとして作用し、大規模な磁場エネルギーの拡散によって離れた環状領域を接続する。
• 磁場増幅は、磁場レイノルズ数Rm ≲10の低値で最適であり、局所的シミュレーションと比較して顕著に低い値である。
• ダイナモ機構は、局所的シミュレーションで観測された誘導過程と整合的であり、垂直運動とせん断に依存するが、非α-Ω平均場機構を介して作用する。
• 飽和は磁場エネルギーの等エネルギー状態に近づき、プラズマベータβ ≈1の強磁化された乱流状態に達する。
• 飽和状態は、局所的熱的平衡の閾値を超えており、長期的な不安定性またはサイクル的挙動を示唆する。
• 遅い段階で、全地球的振動ダイナモモードおよび密度変調（バンド）が出現し、非線形フィードバックまたは数値的アーチファクトと関連している可能性がある。