Idealized Global Models of Accretion Disks with Strong Toroidal Magnetic Fields

技術要約：強固なトロイダル磁場を持つ理想化された全球的降着円盤モデル

問題提起
角運動量輸送の性質、特に正味の垂直磁束を持たない強固なトロイダル磁場に支配された降着円盤については、依然として議論の対象となっている。ゼロ正味垂直磁束（ZNVF）を用いた局所的なシヤリングボックス・シミュレーションでは、局所ダイナモによって維持される「低 $\beta$ 状態」（プラズマ $\beta \sim 1$）が特定されているが、従来の全球的シミュレーションでは、極めて磁化された状態（$\beta \ll 1$）になるか、あるいは浮力による損失（パーカー不安定性）のために強固なトロイダル磁場を維持できない状態になることが多かった。先行する全球的研究（例：Hopkins et al. 2024a,b; Guo et al. 2024）では、しばしば $\beta \ll 1$ かつ $\alpha \gg 1$ という結果が得られていたが、矮新星やX線トランジェントからの観測的推論は、中程度の磁化を持つ円盤に対応する $\alpha \sim 0.1 - 0.3$ を示唆している。局所的な結果と全球的な結果の間の定量的な不一致の起源、および $\beta \sim 1$ のトロイダル支配的な円盤を維持するために必要な条件は、完全には理解されていない。

手法
著者らは、AthenaKを用いた理想化された3次元全球的磁気流体力学（MHD）シミュレーションを提示している。セットアップには、固定された点質量重力ポテンシャル（ニュートン力学、一般相対論はなし）と、トロイダル磁場を持つ回転する磁化ガス雲が含まれる。ガスは初期状態では平衡状態になく、流入、円環化し、円盤を形成する。

• 状態方程式： 局所的に等温の状態方程式が使用され、ガスの熱的スケールハイト $H_{th}$ を半径に対して固定している（$H_{th} = \sqrt{2}hR$）。これにより、著者らは熱的スケールハイトに対する解像度を制御することができる。
• 初期条件： 初期ガス雲は、広範囲の半径において一様な密度と角運動量を持つ。磁場は、円環化半径（$R_{circ} = 16 r_0$）の外側では純粋なトロイダル磁場として初期化されており、広範囲においてプラズマ $\beta_\phi = 1$ である。一部のランでは、比較のために正味の垂直磁場を含めている。
• 数値解像度： シミュレーションには、最大9レベルの適応格子細分化法（AMR）を利用している。重要なパラメータは、熱的スケールハイトあたりのセル数（$h r / \Delta x$）である。著者らは、$h$ のパラメータ研究と解像度の検証を行い、収束性をテストしている。
• 境界条件： 内側の真空シンクは、$r < r_0$ に進入するガスを除去し、外側の境界は初期条件に固定されている。

主な結果

1. 解像度に依存する磁気的遷移： 最も重要な発見は、ガスの熱的スケールハイトが解像されているかどうかに基づいて、円盤の磁気状態が決定的に依存することである。

   • 未解像（$h r / \Delta x \lesssim 3$）： 円盤は高度に磁化され（$\beta \ll 1$）、高い降着ストレス（$\alpha \gg 1$）を示す。円盤は磁気圧によって垂直方向に支えられている。
   • 解像済み（$h r / \Delta x \gtrsim 4$）： 円盤は、内側から外側に向かって垂直方向に崩壊し、$\beta \sim 1$ かつ $\alpha \sim 0.1$ の中程度の磁化を持つ中間面を形成する。この状態において、中間面はガス圧支配（熱平衡に近い）であり、中間面の上部および下部の領域は強固に磁化されている（$\beta \ll 1$）。

2. 円盤の構造と進化：

   • 系は、当初、急速に降着する、強固に磁化されたケプラー回転円盤を形成する。
   • 円盤が進展するにつれ、中間面が乱流的（$\delta v^2 \sim \delta v_A^2 \sim v_A^2 \sim c_s^2$）な変動を持つ準定常状態へと遷移する。
   • 角運動量輸送は、平均的なストレスと変動的なストレスの両方によって駆動される。中間面では変動的なマックスウェル・ストレスが支配的であるが、より高い高度では平均的なストレスが支配的となる。

3. 磁場を維持するメカニズム：

   • 方位角方向の磁場は、浮力によって垂直方向に継続的に脱出するが、局所ダイナモによって補充される。
   • 解像された円盤（$R < R_{circ}$）において、ダイナモは（大規模な流入によって生成される）動径方向磁場の剪断と、乱流拡散のバランスによって駆動される。磁場の動径方向への移流は、内側円盤では副次的であるが、円環化半径の外側では重要となり、外側領域における $\beta \ll 1$ 状態の維持に寄与する。
   • 解像されたケースにおける全球的な結果は、局所的なZNVFシヤリングボックス・シミュレーション（Squire et al. 2024）で見られる「低 $\beta$ 状態」と概ね一致している。

4. ロバスト性テスト：

   • $\beta \sim 1$ の状態は、初期磁場強度、初期乱流の有無、および異なる再構成手法（PLM vs. PPM）の利用に対してロバストである。
   • 正味の垂直磁場を含めると、より強く磁化された中間面（$\beta \lesssim 1$）を持ち、平均マックスウェル・ストレスによって駆動される異なる降着モードが生じる。これは、従来のNVF研究と一致している。
   • メッシュ細分化テストにより、$\beta \sim 1$ への遷移は数値的な解像度の効果であることが確認された。すなわち、解像度が不足している円盤を細分化すると $\beta \sim 1$ 状態への崩壊が誘発され、解像度を下げる（derefine）と $\beta \ll 1$ に戻る。

意義と主張
本論文は、$\beta \sim 1$ のトロイダル支配的な円盤の存在に関する、局所シミュレーションと全球的シミュレーションの間の不一致を解決したと主張している。著者らは、以前の全球的シミュレーションがこの状態を捉えきれなかったのは、熱的スケールハイトが十分に解像されていなかったためであると考えている。

• 観測との関連性： 結果として得られる降着ストレス $\alpha \sim 0.1$ は、矮新星やX線トランジェントの噴出から推定される値と一致しており、このような円盤が、正味の垂直磁束を必要とせずにこれらの現象を自然に説明できる可能性を示唆している。
• 物理的洞察： 本研究は、熱構造が解像されている限り、全球的な設定においても強固なトロイダル磁場が浮力による損失に対して局所ダイナモによって維持され得ることを示している。
• 注意点： 著者らは、結果が厳密に収束しているわけではないことを述べている。局所的なシヤリングボックス研究によれば、さらに高い解像度（$\sim 10$ セル/スケールハイトが必要）を用いると、$\alpha$ と $\beta$ はさらに進化し続ける可能性がある。また、本シミュレーションは理想化されており（ニュートン力学、等温）、放射輸送、自己重力、および一般相対論的効果を無視している。これらは現実的な天体システムにおける定量的な詳細を変化させる可能性がある。これらの結果と他の全球的シミュレーション（例：Hopkins et al. 2024b）との違いは不明確であり、その原因は解像度や熱力学および流入特性の違いにある可能性がある。