Pressure-regulated feedback-modulated star formation as a subgrid model for galaxy formation simulations

技術要約：サブグリッドモデルとしての圧力制御型フィードバック変調星形成

問題提起
銀河形成の大規模な宇宙論的シミュレーション（IllustrisTNG、EAGLE、SIMBAなど）は、銀河集団の統計的性質を再現することには成功しているが、巨大分子雲（GMC）や多相中間物質（ISM）の物理的スケールを解像することはできない。その結果、これらのシミュレーションは、星形成率（SFR）と有効な状態方程式（EoS）に対して、経験的に較正されたサブグリッド処方箋に依存している。従来のモデルは通常、比星形成率がガス密度に $\dot{m}_*/m_g \propto \rho^{1/2}$ （逆自由落下時間）として比例すると仮定しており、その正規化係数は局所宇宙の観測に基づいて設定されている。

著者らは、このアプローチにおける2つの主要な限界を特定している：

1. 物理的結合の欠如： 経験的なスケーリングは、星形成を制御する物理プロセス、具体的には重力とステラ・フィードバックのバランスとの直接的な関連性を欠いている。
2. 解像度および環境への敏感さ： これらの処方箋は、ガス密度や圧力が局所宇宙よりも著しく高い領域（例：高赤方偏移銀河）において、しばしば失敗する。また、局所的な体積密度に依存するモデルは、数値解像度に敏感である。ガス円盤のスケールハイトが解像されていない場合、測定される体積密度は不正確になり、異なる解像度間で一貫性のない星形成予測を招くことになる。

手法
本論文では、圧力制御・フィードバック変調（PRFM）理論（Ostriker & Kim 2022）に基づく新しいサブグリッド・フレームワークを導入している。このフレームワークでは、SFRはISMの熱的および垂直方向の動的平衡を維持するように調整される。ここで、全中間面圧力（$P_{tot}$）はISMの重み（$W$）と均衡し、フィードバック収率（$\Upsilon_{tot}$）を介して単位面積あたりのSFR（$\Sigma_{SFR}$）に比例する。

著者らは、TIGRESSスイートによる高解像度磁気流体力学シミュレーションを用いて較正された、2つの補完的なバージョンのモデルを移動格子コード Arepo に実装した。

1. PRFM-vol (Volumetric/体積型): ガス円盤の垂直方向の厚さが数値的に解像されている（通常、質量解像度 $\sim 10^5 M_\odot$）シミュレーション向け。

   • TIGRESSで較正されたEoS（$P_{tot} \propto n_H^{1.8}$）を用いて有効圧力（$P_{eff}$）を各セルで計算する。
   • 動的時間（$t_{dyn}$）は、ガス、星、およびダークマターの局所的な体積密度から導出される。
   • 枯渇時間（$t_{dep}$）は、体積法則の垂直方向の積分を考慮するための再正規化因子（$R_f \approx 2$）に基づき設定される。

2. PRFM-int (Integrated/積分型): 円盤のスケールハイトが解像されていない粗い解像度（通常、質量解像度 $\gtrsim 10^6 M_\odot$）向け。

   • 局所的な体積密度に依存する代わりに、表面密度（$\Sigma_g, \Sigma_*$）とISMの重み（$W$）を用いて、垂直平衡の検討から中間面密度と圧力を再構成する。
   • 均衡ガススケールハイト（$H_g$）を解き、これらの積分量から直接 $t_{dep}$ を導出することで、解像度に敏感な体積密度の使用を回避し、SFRの処方箋を決定する。

主な貢献

• 物理に基づいた較正： 本モデルは、局所銀河への経験的なフィットではなく、TIGRESSシミュレーションに対して較正されたPRFM平衡条件（$W = P_{tot} = \Upsilon_{tot}\Sigma_{SFR}$）からSFRとEoSを導出している。
• 解像度に対する堅牢性： PRFM-int アプローチの開発は、体積モデルが低解像度で失敗するという問題を具体的に解決するものである。これは、表面密度と平衡制約を利用することで、垂直方向の円盤構造を解像できないことに対する感度を回避している。
• 可変的な効率： 自由落下時間あたりの効率が一定であると仮定する従来のモデルとは異なり、PRFMモデルは、星形成効率（$\epsilon_{dyn}$）が圧力とともに増加する（$\epsilon_{dyn} \propto P_{tot}^{0.43}$）ことを予測する。

結果
著者らは、質量解像度が $10^5, 10^6, 10^7 M_\odot$ の孤立した銀河系類似ディスクシミュレーションにおいて、これら2つの実装を標準的なIllustrisTNGの処方箋と比較検証した。

• 高解像度における一致： $10^5 M_\odot$ において、PRFM-volとPRFM-intは互いに密接に一致し、SFRおよび枯渇時間に関してTIGRESSの結果とも一致した。両者は、TIGRESSで較正されたOstriker-Kim関係（$\Sigma_{SFR} \propto P_{tot}^{1.2}$）を正常に再現した。
• 体積モデルの解像度感受性： PRFM-volとTNGモデルは、顕著な解像度依存性を示した。解像度が $10^6$ および $10^7 M_\odot$ に粗くなると、円盤スケールハイトを解像できないことが原因で、不正確な体積密度測定が生じる。これにより、PRFM-volモデルは枯渇時間を過大評価し、高解像度のSFR履歴に一致しなくなる。
• PRFM-intの堅牢性： PRFM-intの実装は、テストされたすべての解像度において堅牢であった。これは、高解像度（$10^5 M_\odot$）のPRFM-volシミュレーションのグローバルな星形成履歴と表面密度分布を、$10^7 M_\odot$ においても正確に追跡することに成功した。
• TNGとの比較： TNGモデルと比較して、PRFMモデルは、特に高圧・高密度領域において、より短いガス枯渇時間を予測する。PRFMモデルは、TNGモデルが自由落下時間あたりの一定の効率を維持するのに対し、枯渇時間の非常に広いダイナミックレンジと、圧力に対する急峻な依存性を示す。

意義
本論文は、PRFM由来のサブグリッド処方箋が、銀河形成シミュレーションのダイナミックレンジ全体にわたって、物理に基づいた数値的に安定したフレームワークを提供すると主張している。

• 高解像度シミュレーションに対して： PRFM-volモデルは、従来の経験的な法則よりも現実的なISM物理の取り扱いを提供し、SFRを銀河の動的バランスに直接結びつける。
• 大規模な宇宙論的シミュレーションに対して： PRFM-intモデルは、決定的な「解像度依存性」の問題を解決する。これにより、円盤スケールハイトを解像できない大規模ボックスの宇宙論的シミュレーションであっても、高解像度の物理と整合性を保ちながら、物理的に動機付けられた星形成則を利用することが可能になる。
• 高赤方偏移への影響： 著者らは、PRFMモデルの圧力依存的な効率が、高赤方偏移銀河（圧力がより高い銀河）における星形成が現在のTNGのようなモデルの予測よりも著しく効率的である可能性を示唆していると述べており、これはJWSTによって観測されている明るい銀河の解釈に影響を与える可能性がある。

本研究は、「Learning the Universe Collaboration」の目標、すなわち、高解像度のISMシミュレーション（TIGRESSのような）が、より広いパラメータ空間（例：金属量の変化）をカバーするように拡張されるにつれて更新可能な、経験的にチューニングされたものではなく「予測的」なサブグリッドモデルを開発するという目標への一歩として提示されている。