Individual Star Sampling in Star Formation Simulations: A Semi-Deterministic Model

技術要約：星形成シミュレーションにおける個別の星のサンプリング：半決定論的モデル

問題提起
現代の銀河形成シミュレーションでは、太陽質量分解能で個々の星を解像することが増えている。しかし、恒星質量を決定する物理プロセス（雲の断片化、降着、フィードバック）は現在のシミュレーションでは解像できないスケールで動作するため、恒星質量を割り当てることは依然として根本的な課題となっている。その結果、ほとんどのシミュレーションは、普遍的で不変な初期質量関数（IMF）を確率的にサンプリング（ランダムまたは「RND」サンプリング）することに依存している。この手法は、埋没クラスターの総質量（$M_{\rm ecl}$）と、その中で最も質量の大きい星の質量（$m_{\star, \rm max}$）との間の緊密な相関関係を再現できないという、重大な観測上の課題に直面している。純粋な確率的サンプリングは、低質量クラスターにおいて $m_{\star, \rm max}$ の過度な散らつきを予測し、本来は存在すべきではない環境において大質量星を過剰に生成してしまう。さらに、確率的サンプリングは、極端な変動を抑制する環境的制約を捉えることができず、フィードバックや星形成効率（SFE）における大きな実行間の差異を生じさせる。

手法：半決定論的（SDT）スキーム
著者らは、計算の実現可能性を維持しながら、恒星質量をその親となるクラスターに結びつけるように設計された、半決定論的（SDT）サンプリングスキームを導入している。この手法は、arepo 磁気流体力学コード内の RIGEL フレームワーク内で動作する。主な構成要素は以下の通りである：

1. リザーバー粒子（RsvPs）： 未解像の崩壊するコアや原始星円盤は、「リザーバー粒子」によって表現される。ガスセルは、特定の基準（ジーンズ不安定性、ビリアル不安定性、収縮流、密度、および温度の閾値）を満たすときに RsvP へと変換される。RsvPs は流体力学ソルバーからは切り離されるが、重力には結合したままである。これらは、星を形成する能力を持つ「活動的」な状態になる前に、「不活性」なフェーズ（期間 $\sim t_{\rm dyn}$）を経る。
2. オンザフライのクラスター識別： 効率的な Friends-of-Friends (FoF) アルゴリズムが、活動的な RsvP および若い星（$< \tau_{\rm young}$）をグループ化するために周期的に（$\Delta t_{\rm FoF}$）適用され、クラスターを特定する。これにより、既存の星と利用可能なリザーバー質量の両方を含む、クラスターの総質量（$M_{\rm ecl}$）の即時的な推定値が得られる。
3. 半決定論的サンプリング・ロジック：
   • 大質量星（決定論的）： 現在の $M_{\rm ecl}$ によって許容される最も質量の大きい星に対して、本手法は「最適サンプリング」を採用する。これは、$m_{\star, \rm max} - M_{\rm ecl}$ 関係に基づき、単一の最も質量の大きい星（$m_1$）の質量を決定論的に割り当て、クラスターの進化ステップごとに、最高質量に対応する質量ビンが正確に一度だけ埋まるようにする。
   • 低質量星（確率的）： その他のすべての星は IMF から確率的にサンプリングされるが、その際の上限質量は、決定論的に割り当てられた大質量星の質量（$m_{1, \rm low}$）によって設定される。これにより、集団の大部分については IMF の形状を維持しつつ、高質量端における環境的制約を強制する。
   • 質量保存： 本手法は、サンプリングされた星の質量が局所的なリザーバー粒子の質量を超えた場合に、探索半径（$R_{\rm search}$）内の隣接する RsvP から質量を降着させることで、局所的な質量保存を確保する近傍ベース（NGB）のアプローチを利用している。

主な結果
SDT 法は、孤立した分子雲（$10^3$ から $8 \times 10^4 M_\odot$ の範囲）および矮小銀河の合体シミュレーションにおいて、標準的なランダム（RND）および近傍ベース（NGB）サンプリングと比較検証された。

• $m_{\star, \rm max} - M_{\rm ecl}$ 関係の再現： SDT 法は、観測されたクラスター質量と最大星質量の間の緊密な関係を自然に再現し、最適サンプリングの予測に密接に従う。対照的に、RND および NGB 手法は顕著な確率的散らつきを示し、低質量クラスター（$M_{\rm ecl} \lesssim 10^3 M_\odot$）において大質量星を系統的に過剰生成する。
• 実行間の差異の低減： SDT シミュレーションは、異なるランダムシード間での大質量星（$m > 8 M_\odot$ および $m > 30 M_\odot$）の数において無視できる程度の変動を示すが、RND シミュレーションは、特に低質量雲において発散的な傾向を示す。
• 形成時期と質量分離： SDT（および NGB）における環境的制約は、ポアソン期待値と比較して大質量星の形成を $\sim 0.15$ Myr 遅らせる。この遅延は、恒星フィードバックの開始におけるコヒーレントなタイムラグをもたらす。さらに、本手法は初期の質量分離を生み出し、大質量星が密度ピークで優先的に形成されることを示すが、これは RND シミュレーションには見られない特徴である。
• フィードバックと星形成効率（SFE）： 大質量星の形成が遅延し規制されているため、SDT シミュレーションはフィードバック光度（EUV、風、超新星爆発）の実行間の変動が最小限である。その結果、統合された SFE は SDT ランにおいてより高く、かつ安定している。これは、ランダムにサンプリングされた大質量星によってガスが早期に破壊されることが回避されるためであり、特に小さな雲において顕著である。
• 銀河規模の IMF（gwIMF）： 矮小銀河の合体シミュレーションにおいて、SDT 法は星形成率（SFR）に対する gwIMF スロープの強い依存性を予測する。低 SFR において、gwIMF は著しく急峻になり（スロープ $\sim -3$）、「トップライト（top-light）」となる。これは、統合銀河 IMF（IGIMF）理論と一致している。RND モデルは、この変動を再現できず、SFR に関わらず標準的なスロープを維持する。
• H$\alpha$ 診断： 本研究は、低 SFR 環境において、大質量星の不足により SDT モデルの比 H$\alpha$ 光度が低くなることを明らかにしている。これは、不変の IMF を仮定する標準的な H$\alpha$-to-SFR 変換係数が、静穏な低質量銀河における固有の SFR を、因子 2–3 の割合で系統的に過小評価することを示唆している。

意義と主張
論文は、SDT スキームが、高解像度シミュレーションにおける恒星質量サンプリングの問題に対する、物理的に動機付けられた計算効率の高い解決策を提供すると主張している。最も質量の大きい星を観測的制約に基づいて決定論的に扱う一方で、集団の大部分については確率的サンプリングを維持することにより、SDT は以下のことを達成している：

1. 確率的サンプリングと観測された $m_{\star, \rm max} - M_{\rm ecl}$ 関係との間の不一致を解消する。
2. フィードバックおよび SFE のランダムシードに対する人工的な感受性を、特に矮小銀河のような低質量環境において低減させる。
3. 銀河規模の IMF が SFR と共に変化することを自己整合的に再現する。これは、技術的な困難のためにシミュレーションではしばしば見落とされる現象である。
4. 初期質量分離やフィードバックのタイミングを含む、星団の初期力学的進化を研究するための、より現実的な枠組みを提供する。

著者らは、現在のモデルは星を衝突のない粒子として扱い、固定された IMF スロープを仮定しており、金属量やガス密度への依存性を無視しているという限界があることを認めている。また、RsvP 内での 100% の星形成効率の仮定は、原始星のフィードバック（ジェットや風）を無視しており、それがさらに SFE を制御する可能性があることも述べている。これらの簡略化にもかかわらず、著者らは SDT 法が、銀河進化シミュレーションにおける星形成の自己整合的なモデリングに向けた重要な一歩であると主張している。