Coexisting Tayler instability-driven dynamos in radiative zones: New dynamo solution and its impacts on stellar physics

技術要約：放射層における共存するテイラー不安定駆動型ダイナモ

問題提起
Kepler、CoRoT、TESSといったミッションによる近年の星震学観測は、恒星核の内部回転速度が、標準的な回転1次元恒星進化モデルの予測よりも著しく遅い（場合によっては数桁低い）ことを明らかにしている。これらのモデルは、通常、射流（meridional circulation）や流体不安定性のみを含む流体力学的なプロセスを組み込んでいるが、観測結果と一致するために必要な効率的な角運動量（AM）輸送を説明できていない。内部重力波と大規模磁場がこの欠落した輸送メカニズムの主要な候補であるが、後者については制約が極めて不十分である。具体的には、テイラー・スプライト・ダイナモ（テイラー不安定性が剪断流によってトロイダル磁場を再生し、角運動量輸送を駆動するメカニズム）は、進化コード（例：Spruit 2002; Fuller et al. 2019）に広く実装されているが、強く成層された恒星放射層に関連する領域での堅牢な数値的確認を欠いている。

手法
著者らは、安定成層したブシネスク・磁気流体力学（MHD）流の3次元直接数値シミュレーション（DNS）を用いて、テイラー・スプライト・ダイナモを調査している。セットアップは、同心球（$r_i$ および $r_o = 4r_i$）の間にある恒星放射層をモデル化しており、ノンスリップおよび電気絶縁の境界条件を用いている。

• 強制力: 固定された境界回転を課す従来のテイラー・クエット構成とは異なり、本研究では体積的な体積力を用いて、シェラー状の微分回転プロファイルを維持している。このアプローチにより、偽の流体不安定性を最小限に抑え、MHD駆動型ダイナモの明確な特定を可能にしている。
• パラメータ: シミュレーションは、フレーム回転率とブルント・ヴァイサラ振動数の比（$\Omega/N$）によって特徴付けられる一連の成層レジームを対象としており、これは $0.0077$から$0.5$まで変化する。これは、弱い成層から極端な成層（$N \approx 130\Omega$まで）までの領域をカバーしている。磁気プラントル数は$Pm = 4$、熱プラントル数は$Pr = 0.1$ に固定されている。数値スキームには、結合MHD方程式を解くためのオープンソースの擬似スペクトル法コードであるMagICが使用されている。異なる解のブランチを探索するために、初期状態として $(\ell=1, m=0)$ または $(\ell=2, m=0)$ のポロイダル磁場を用いている。

主な貢献と結果

1. 共存するダイナモ解の発見:
   本研究は、同一のパラメータ空間内で、初期磁場幾何学および空間的位置によって区別される、2つの異なる共存するテイラー・スプライト・ダイナモ解を初めて特定したことを報告している。

• 赤道ダイナモ: 初期 $(\ell=2, m=0)$ ポロイダル場によって誘発されるこの解は、弱成層レジーム（$\Omega/N \gtrsim 0.11$）で動作する。これを駆動する不安定性は、方位角磁気回転不安定性（MRI）とテイラー不安定性の両方の特性を示し、赤道面付近で不安定モードが発達する。
• 極ダイナモ: 初期 $(\ell=1, m=0)$ ポロイダル場によって誘発されるこの解は、極端な成層を含むはるかに広い範囲（$\Omega/N \gtrsim 0.0077$）で動作する。これは標準的なテイラー不安定性によって駆動され、不安定モードは極軸付近に局在している。この解は、MRIの安定性条件が破られた場合でも持続する。

2. 磁場および輸送に関するスケーリング則:
   著者らは、DNSデータに基づいて校正された、飽和磁場成分および輸送係数の新しいスケーリング則を導出した。

• 磁場: 方軸対称トロイダル磁場（$B_\phi^{m=0}$）は、以前の解析的予測（例：Fuller et al. 2019）が示唆した $( \Omega/N_{eff} )^{1/3}$ やそれ以上の依存性とは異なり、有効ブルント・ヴァイサラ振動数（$N_{eff}$）に依存しないことが判明した。径方向磁場（$B_r^{m=0}$）は、Fuller et al. (2019) によって導出されたスケーリングに従う。
• 最小剪断（$q_{min}$）: $B_\phi^{m=0}$ に関する新しいスケーリングに基づき、著者らはテイラー不安定性を引き起こすために必要な最小剪断率の新しい処方箋を推論した。この新しい閾値は、特に強成層レジームにおいて、以前の解析的推定値よりも著しく低くなっており、ダイナモが以前は抑制されていると考えられていた領域でも動作できることを示唆している。
• 角運動量輸送: 輸送はマックスウェル応力（粘性 $\nu_M$）によって支配されており、これは $\nu_M \propto (\Omega/N_{eff})^{9/4}$ とスケーリングする。これは Fuller et al. (2019) による $(\Omega/N_{eff})^2$ という予測よりもわずかに効率が低いが、Spruit (2002) による $(\Omega/N_{eff})^4$ よりも遥かに効率的である。レイノルズ応力は、マックスウェル応力の約 $3 \times 10^{-4}$ であり、無視できるほど小さいことが示された。

3. 星震学への影響:
   本研究は、生成された磁場の観測的シグネチャーを評価している。

• 赤色巨星: 生成された径方向磁場は、赤色巨星の水素燃焼殻で検出されている $10^4 - 10^5$ G の磁場を説明するには弱すぎる。さらに、トロイダル磁場の卓越は、振動周波数の磁気シフト（$\delta\nu_{mag}$）が $\nu^{-1}$ に比例することを意味するが、これは現在の赤色巨星における観測結果と矛盾する。
• 主系列星: $\gamma$ Doradus 型の星において、このダイナモは強いトロイダル磁場（$\sim 10^6$ G）と径方向磁場（$\sim 10^4$ G）を生成し得る。著者らは、これらの磁場が磁気重力慣性（MGI）モードの診断を通じて検出可能である可能性を示唆しており、ダイナモの存在を制約する手段を提供できる可能性がある。

意義と主張
本論文は、初期磁場幾何学が極ブランチを支持する場合、テイラー・スプライト・ダイナモは進化した恒星に見られる極端な成層レジームにおいても動作可能な、堅牢なメカニズムであると主張している。極ダイナモが $\Omega/N \approx 0.0077$ で機能できるという発見は、強い成層がこのメカニズムを抑制するという従来の仮定に疑問を投げかけるものである。

著者らは、特に最小剪断（$q_{min}$）の低い閾値に関する新しいスケーリング則が、テイラー・スプライト・ダイナモが以前のモデルよりも広い恒星内部の体積で動作できることを意味していると強調している。この拡張された動作範囲は、このメカニズムが赤色巨星の内部回転を説明するのに十分なほど効率的である可能性を示唆しており、これはこれまで恒星進化モデリングにおける課題となっていた問題である。しかし、著者らは亜巨星段階については慎重な姿勢を保っており、輸送が効率的すぎて、さらなる校正なしには亜巨星の回転率と一致しない可能性があると述べている。

最後に、本研究は2つのダイナモ解のバイスタビリティ（双安定性）がもたらす複雑さを強調している。異なる不安定性の混合（MRI/テイラー vs 純粋なテイラー）および異なる初期条件によって駆動される2つの異なるブランチの存在は、恒星磁場の進化が初期状態に非常に敏感である可能性を示唆しており、これは現在の1次元進化コードでは捉えきれない要因である。著者らは、$\nu_M$ と $q_{min}$ の新しい処方箋を将来の恒星進化グリッドに組み込み、観測された回転プロファイルや表面組成を再現するダイナモの能力をより良くテストすべきであると結論付けている。