Resolving Star Cluster Formation in Galaxy Simulations with Cosmic Ray Feedback

技術要約：宇宙線フィードバックによる銀河シミュレーションにおける星団形成の解明

問題提起
星団は、大規模な星形成の主要な場であり、それが間接的に星形成を制御し、星間物質（ISM）を形作るフィードバック過程を駆動する。銀河のフィードバックを包括的に理解するには、パーセクスケールの星団形成と多相的なISMを解像する必要がある。宇宙線（CR）はISMの圧力と加熱に寄与することが知られているが、その星団形成率（SFR）、星団質量関数、および星団の束縛状態に対する具体的な影響については依然として不明である。先行研究では、銀河流出の文脈で宇宙線フィードバックを調査したり、星団に対してサブグリッドモデルを使用したりしてきたが、個々の星粒子、動的に結合した宇宙線、およびフィードバックによって制御された多相的ISMを同時に解像した研究はほとんどない。本研究は、自己閉じ込め理論に基づく宇宙線輸送モデルが、星団の形成とその特性にどのように影響するかという理解の空白を埋めるものである。

手法
著者らは、移動メッシュコード Arepo を用いた高解像度の「tallbox」シミュレーションを分析している。シミュレーションには、水素と金属の冷却・加熱プロセスを追跡する Crisp 非平衡熱化学・フィードバックモデルを採用している。本研究では、宇宙線輸送の影響を分離するために、3つの異なる物理ケースを利用している：

1. MHD: 宇宙線を一様な加熱/電離項として扱う磁気流体力学。
2. CR-NL: 2モーメント宇宙線磁気流体力学（CRMHD）を用い、非線形ランダウ減衰（NLLD）のみを減衰機構とする。これにより、宇宙線とガスの結合が一様になる。
3. CR-NL-IN (Fiducial): CRMHDを用い、NLLDとイオン中性減衰（IND）の両方を適用する。このケースでは、宇宙線は冷たい中性ガスからデカップル（分離）し、星形成への局所的な影響が減少する。

シミュレーションは、質量分解能 $10 M_\odot$、空間スケール $\sim 0.5$ pc までガスを解像している。星形成は、ガスセルが密度閾値（$n_H = 10^3 \text{ cm}^{-3}$）を超えた場合に、自由落下時間あたり100%の局所効率で衝突のない星粒子へと変換される、シュミット型の処方箋を通じて実装されている。恒星フィードバックは、宇宙線と超新星（SNe）の影響を分離するため、超新星のみを通じてモデル化されている。放射フィードバックとステラウィンド（星風）は省略されている。SNeは、年齢3〜38 Myrの星に対して確率的に注入され、CRMHDケースでは、SNeエネルギーの5%が追加で宇宙線として注入される。

星団は、時間座標を組み込んだ修正 Friends-of-Friends (FoF) アルゴリズムを用いて特定される。これにより、特定されたグループが、2 pcの空間連結長内にある同齢の星の集団であることを保証している。解析は、大きな潮汐破壊を受ける前の初期特性を特徴付けるため、「若い」星団（年齢 $\sim 10$ Myr）に焦点を当てている。

主な結果

1. 星団質量関数: すべての3つのシミュレーションケースにおいて、星団質量関数は質量 $> 10^3 M_\odot$ においてべき乗則（$dN/dM \propto M^{\alpha_{PL}}$）に従い、傾きは $-2.0$から$-2.4$の間となる。これは、傾きが星形成率表面密度（$\Sigma_{SFR}$）と正の相関を持つという観測的傾向と一致している。宇宙線の導入は、これらの観測された関係に沿ってデータポイントをシフトさせるが、非物理的な結果を生むことはない。
2. 星団形成効率 (CFE): 大規模な星団（$M > 5000 M_\odot$）について計算した場合、CFE（$\Gamma$）は $\Sigma_{SFR}$ と正の相関を示し、これは理論的予測（例：Kruijssen 2012）および観測と一致している。しかし、より低い質量閾値（$M > 1000 M_\odot$）については、CRMHDケースは高い一定のCFEを示す一方で、MHDケースはより低い値を示す。著者らは、CRMHDケースにおける高いCFEは、使用した密度閾値と効率の高さによるアーティファクト（人工的な結果）である可能性を示唆している。
3. 超新星のクラスタリング: すべてのケースにおいて、約50%のSNeは希薄で前処理されたガス（$n_H < 10^{-2} \text{ cm}^{-3}$）の中で発生している。この割合は宇宙線モデルにほとんど依存しておらず、宇宙線フィードバックが、最初の数個のSNeが爆発する局所的な環境を大きく変えることはないことを示している。これは、宇宙線が希薄なISMの全容積比率を変化させているにもかかわらずである。
4. 星団の特性とビリアルパラメータ:
   • 半径と速度分散: MHDケースの星団は、CRMHDケースと比較して、より大きな半質量半径と高い1次元速度分散を示す。
   • 束縛状態: MHDの星団は平均して非束縛状態（$\alpha_{vir} > 1$）であるが、CRMHDの星団は平均して束縛状態（$\alpha_{vir} < 1$）である。
   • メカニズム: 著者らは、この違いをISMの乱流エネルギー収支に帰している。MHDケースはSFRが最も高く、その結果、SNeによるエネルギー注入率が高くなり、より乱流的なISMとなる。この高い周囲の乱流が形成される分子雲に引き継がれ、その結果、速度分散が高く、重力的に束縛されにくい星団が生じる。宇宙線フィードバックはSFRを減少させ、それによって乱流エネルギー収支を低下させ、より強く束縛された星団を生成する。

意義と主張
本論文は、動的に結合した宇宙線、多相的ISM、および大質量星を代表する個々の星粒子を同時に含む、初のシミュレーションを提示していると主張している。主な意義は、以下のことを示した点にある：

• Crispフィードバックモデルは、放射フィードバックの事前注入がない状態でも、微調整なしに、観測された質量関数勾配や環境依存性を持つ星団形成を自己整合的に再現できること。
• 宇宙線フィードバックの主要な影響は、全般的なSFRの減少であること。この減少はISMの乱流エネルギー収支を低下させ、その結果、効果的な宇宙線フィードバックがない場合と比較して、より重力的に束縛された（速度分散が低い）星団を生成すること。
• 宇宙線は、個々の高密度分子雲の崩壊を動的に阻止したり、SNeの局所的な前処理を著しく変化させたりはせず、その影響はISMの状態の全般的な制御を通じて媒介されること。

著者らは、理想化された設定により放射フィードバックとステラウィンドが省略されているものの、宇宙線の影響に関する各ケース間の相対的な差異は堅牢であると結論付けている。これらの結果は、宇宙線輸送モデルに対して新たな制約を与えるとともに、宇宙線フィードバックが、未知の領域を作り出すのではなく、観測された星団の関係性に沿ってシミュレーション結果を移動させる役割を果たすことを示唆している。