Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

技術要約：星および円盤形成過程におけるダストの凝集、動的ドラッグ、および乱流混合のモデリング

問題提起
ダストは、ガスの冷却、電離レベル、および磁場とガスの結合に影響を与えるため、星および惑星形成の両方において極めて重要な役割を果たす。より大きな粒子（ペブル）へのダストの凝集は、微惑星形成やコア成長に不可欠であるが、これまでの星形成に関する三次元流体計算では、ガスに対するダストの動的な進化がほとんど無視されてきた。Bate (2022) による従来の手法などは、ダスト集団を各流体要素（SPH粒子）内で独立して進化させ、ダストがガスに完全に結合していると仮定していた。この近似は、初期の低密度な崩壊フェーズでは成立していたが、密度が高く粒成長が急速に進む原始惑星系円盤内のような、より長いタイムスケールにおけるダストの移動、沈降、および混合を考慮するには不十分であった。さらに、既存のモデルは、粒成長に必要な相対速度を駆動する乱流と、ダストを垂直方向に混合する乱流拡散を、いかに自己整合的に扱うかという点においても欠けていることが多かった。

手法
著者らは、以下の3つの異なる物理プロセスを組み合わせた、SPHNGという新しいSPH（粒子法）コードの実装を提示している。

1. ダストの凝集: 離散化された平均場スモルホフスキー方程式法（Hirashita & Yan 2009; Bate 2022）を利用し、コードはすべてのSPHガス粒子に対して、$N$ 個の対数的に間隔を置いたサイズビンにわたるダスト粒子の集団の進化を追跡する。この手法は、特定の凝集閾値速度以下での衝突確率を1として、二体衝突を考慮する。
2. 動的ドラッグ (MULTIGRAIN): 著者らは、MULTIGRAIN アルゴリズム（Hutchison et al. 2018）を統合している。これは、終端速度近似を用いて、ガスと複数のダスト種の混合に関する連続体流体方程式を解くものである。これにより、圧力勾配や重力に起因する、異なるサイズのダストとガスの間のバルク的な相対速度を計算することが可能となる。
3. 暗黙的積分および乱流拡散: 主要な革新は、ダストのドラッグ方程式を解くために、最近開発された暗黙的積分法（Elsender & Bate 2024）を使用している点である。これにより、停止時間が長い粒子に伴う制限的なタイムステップ制約を回避し、物理的に不自然な負のダスト分率の発生を防ぐことができる。さらに、著者らは乱流拡散をモデル化するための新しい暗黙的手法を導入している。ダストの乱流攪拌を拡散プロセス（$D_i = \alpha_{SS} c_s H / Sc$）として扱い、シュミット数（$Sc$）によって異なる粒径の結合効率を考慮する。この拡散項は、ガウス＝ザイデル・オペレータースプリッティング法を用いて解かれ、流体要素間でのダストの混合を可能にする。

本コードは、分子雲コアの崩壊から、第一水 hydrostatic コア相、そして原始惑星系円盤の形成と進化に至るまでのダストの進化を、シンク粒子を用いて恒星コア形成後まで追跡する。

主な貢献

• 初の3D結合シミュレーション: 本研究は、星およびその円盤の形成初期段階において、ダストの凝集、動的ガスドラッグ、および乱流混合を同時にモデル化した、初の三次元流体計算である。
• 暗黙的拡散スキーム: 乱流によるダストの混合を安定してシミュレートできる、堅牢な暗黙的ソルバーの開発により、明示的な拡散スキームに伴う厳しいタイムステップ制限を回避した。
• 統一された扱い: コードは、相対速度（成長を駆動）と乱流拡散（混合を駆動）の扱いを、テストケースにおいて同じ乱流パラメータ（$\alpha_{SS}$）を用いて、単一のフレームワーク内で統一している。

結果
著者らは、テスト計算（垂直カラムにおけるダストの沈降）と完全な星形成シミュレーションの両方の結果を提示している。

• テスト計算: コードは、ダストの沈降と乱流拡散の間の平衡に関する解析解を正常に再現した。テストにより、乱流拡散を含めることで、非現実的に薄いダスト層の形成が防止され、有限のダストスケールハイトが維持されることが確認された。
• 星形成シミュレーション:
  • 第一水 hydrostatic コア: 前報（Bate 2022）と一致して、第一水 hydrostatic コア内では急速なダスト成長が起こり、コア内部では局所的に単分散な集団（サイズ100 $\mu$m超）が生成される。ドラッグと拡散の導入は、この初期フェーズにおいては進化のタイムスケールが短いため、無視できる程度の影響であった。
  • 原始惑星系円盤の進化: 円盤が形成されると、乱流拡散の導入がダストの進化を著しく変化させる。
    • 成長の促進: 乱流拡散は、ダストの成長率を大幅に高める。拡散が、小さな粒子や中間サイズの粒子を、大きな粒子が成長している高密度領域へと継続的に供給することで、凝集のための「餌」の局所的な枯渇を防ぐためである。
    • 乱流への依存性: 乱流のレベル（$\alpha_{SS}$）が極めて重要である。高い乱流（$\alpha_{SS} = 10^{-3}$）では、円盤全体にわたって $\sim30-100$ $\mu$m の単分散な集団への急速な成長が起こる。より低い乱流（$\alpha_{SS} = 10^{-4}$）では、成長は遅くなり、より垂直方向に層状化したダスト集団となり、大きな粒子は中間面付近に限定される。
    • ダストの沈降: 弱い乱流を伴うシミュレーションでは、ダストの沈降により中間面のダスト・ガス比が大幅に増加（>50%）する一方で、小さな粒子は中間面付近で枯渇し、表面層に留まる。
  • ドラッグ vs 拡散: 特定のプロセスを含む有無による比較では、ダストがよく結合しているため、これらのタイムスケールにおいてダスト・ガスの動的ドラッグはダスト分布にほとんど影響を与えないことが示された。対照的に、乱流拡散の導入は、ダスト分布を滑らかにし、凝集を加速させるという極めて大きな効果をもたらした。

意義と主張
本論文は、この新しい手法が、星およびその円盤の形成過程におけるダストの成長、混合、および移動を、初めて自己整合的にモデル化できることを主張している。著者らは、ドラッグと拡散の導入が第一水 hydrostatic コアにおけるダスト成長に関する根本的な知見を変えることはないものの、その後の原始惑星系円盤の進化においては極めて重要であることを強調している。具体的には、乱流ダスト拡散が、大きな粒子が形成されている領域へ小さな粒子を供給する源として機能することで、ダスト粒子の成長率を大幅に高めると主張している。

著者らは、研究の範囲について以下のように控えめな姿勢を保っている。

• 結果は、採用された乱流モデル（例：$\alpha$-円盤 vs MHD乱流）に依存する。
• 現在のシミュレーションは凝集のみを考慮しており、分裂（フラグメンテーション）は考慮していない。異なる凝集閾値の設定によっては、分裂と小さな粒子の再生成が生じる可能性がある。
• 放射伝達の計算は、依然として標準的な星間塵のオパシティを仮定しており、進化する局所的な粒径分布を使用していない。これは、将来的な拡張が必要なステップとして特定されている。

最終的に、本研究は、エッジオン（エッジオン）円盤におけるダストの沈降や乱流攪拌に関して、ALMA、VLT-SPHERE、JWSTといった施設からの豊富な新しい観測データと、三次元流体計算を比較するためのツールを提供するものである。