The Entangling of Supernova Feedback Impacts with Coarsening Simulation Resolution

技術要約：超新星フィードバックの影響と粗いシミュレーション解像度の絡み合い

問題提起
超新星（SN）フィードバックは、星形成（SF）を制御し、銀河流出を駆動する主要なメカニズムとして広く認識されている。しかし、SNがいかなる具体的なチャネルを通じてこれらのマクロな効果を達成するのかは、依然として不明である。高解像度のシミュレーションは、超新星残骸の熱膨張と放射冷却を明示的に解像できるが、ほとんどの銀河形成シミュレーションはより低い解像度で動作している。これらの低解像度モデルは、単一のガスセル内の星間物質（ISM）の状態は均質であり、シミュレーションで利用可能なデータから導出可能であると仮定する「サブグリッド」近似に依存している。ここで扱われる中心的な問題は、粗い解像度のシミュレーションが、異なる環境（例：希薄な領域 vs 密な領域）で発生するSNの明確に異なる物理的影響を自己整合的に区別できるのか、あるいは、粗いグリッドによる平均化がこれらの重要な物理プロセスを覆い隠し、冷却半径や星形成率の誤推定を招くのか、ということである。

手法
著者らは、孤立した矮小銀河を用いた2つの独立した超高解像度シミュレーション、LYRA（ガス質量解像度 $m_{gas} \sim 4 M_\odot$）および RIGEL（$m_{gas} \sim 1 M_\odot$）を分析している。両シミュレーションは、移動メッシュ流体計算コード AREPO を使用しており、同様の初期条件（低バリオン分率、$M_{vir} = 2 \times 10^{10} M_\odot$）を共有している。極めて重要な点として、両シミュレーションは、SNによるエネルギー注入と放射冷却の影響のみを分離するために、若い星からの放射フィードバックを除外している。

本研究は、個々のSNイベントの直近の環境（プロジェニター星の周囲の最も近い $M_{env}$ のガスによって定義される）に焦点を当てている。著者らは以下の手順を行った：

1. 局所環境の特性評価： 両シミュレーションにおける局所的なガス密度（$\rho$）と比エネルギーの変化（$\Delta e$）を測定した。
2. 解像度の影響の分析： 高解像度データを、局所環境の質量スケール（$M_{env}$）を $20 M_\odot$（微細）から $200 M_\odot$ および $2000 M_\odot$（粗い）まで変化させることで、系統的に「粗化」した。これは、低解像度シミュレーションにおける平均化プロセスを模したものである。
3. グローバルな効果との相関： SNの特性を、銀河流出率（ディスク上方の特定高度におけるエネルギーフラックス）およびSNイベント前後における局所的な星形成率（SFR）と比較した。
4. 解析モデルによる検証： シミュレーションデータを、SN衝撃波の進化に関する解析解（例：Kim & Ostriker 2015のセドフ・テイラー相、冷却半径の推定）と比較した。

主要な貢献と結果

• SN環境の二峰性分布： 分析により、SNが発生する局所的な密度における明確な二峰性分布が明らかになった。

  • 低密度モード： 希薄な領域（$n \sim 10^{-4} - 10^{-1} \text{ cm}^{-3}$）で発生するSN。これらのイベントは、大きな比エネルギー変化（$\Delta e$）を示し、しばしばディスクの脱出速度を超える。これらは銀河流出と時間的に相関しており、放射損失が支配的になる前にキロパーセックスケールまで膨張することが可能である。
  • 高密度モード： 密な星形成領域（$n \sim 10^4 \text{ cm}^{-3}$）で発生するSN。これらのイベントは、効率的な放射冷却により、より小さな $\Delta e$ を示す。これらはパーセクスケールで局所的な星形成塊を効果的に破壊するが、有意な流出を駆動することはない。

• 明確に異なるフィードバックチャネル： 本論文は、これら2つのモードが別々の創発的挙動を駆動することを確立している：

  • 低密度SN は、銀河流出（銀河周辺物質へのエネルギー注入）の主要な駆動力である。
  • 高密度SN は、直接的な星形成抑制（高密度ガスの破壊）の主要な駆動力である。
  • 決定的なことに、低密度SNは、周囲のガスが希薄すぎるため、その直近の星形成を有意に抑制することはない。また、高密度SNは、シェルがディスクを脱出する前にエネルギーが放射されてしまうため、効率的に流出を駆動することはない。

• 解像度の閾値： この2つのモードの区別は、非常に高い解像度（$M_{env} \lesssim 100 M_\odot$）においてのみ保持される。

  • $M_{env}$ が $200 M_\odot$ に増加すると、密度の二峰性が失われ始める。
  • $M_{env} = 2000 M_\odot$（多くの宇宙論的シミュレーションに典型的）では、2つのモードは区別不可能になる。平均化プロセスによって密度のコントラストが消失し、粗いデータのみから、あるSNが流出を駆動するのか、あるいは星形成を抑制するのかを判断することは不可能になる。

• 冷却半径の誤推定： 粗い解像度の密度推定（$\rho_{2000}$）は、SNの真の冷却半径（$r_{cool}$）を予測できない。

  • 与えられた粗い密度に対して、真の冷却半径（微細な密度から導出される）は、数桁単位で変動し得る。
  • 粗い推定値は、低密度環境において $r_{cool}$ を過小評価する傾向があり、また高密度環境における効果の二峰性を捉えることができない。これは、熱エネルギーを注入するか運動量を注入するかを決定するために局所密度に依存するサブグリッドモデルを根底から揺るがすものである。

• 星形成閾値： 粗い質量スケール（$2000 M_\odot$）で平均化すると、平均密度がわずか $10^{-2} \text{ cm}^{-3}$ 程度の領域でも非ゼロの星形成率を示すことがある。これは、高解像度シミュレーションにおいて分子雲を模すために用いられる高密度閾値（$\rho_{th} \sim 10^3 - 10^4 \text{ cm}^{-3}$）が、混合相を含む粗い解像度のセルに適用した場合、人工的に制限的すぎる可能性を示唆している。

意義と主張
著者らは、SNフィードバックのマクロな効果は一様ではなく、ISMの局所的な小スケールの構造に強く依存すると結論づけている。本論文の主張は以下の通りである：

1. 自己整合性には高解像度が必要である： SNフィードバックの二重の役割（流出とSF抑制）を自己整合的にモデル化しようとするシミュレーションは、セルあたりの質量スケールを $\lesssim 100 M_\odot$ まで解像しなければならない。この閾値を下回ると、フィードバックの明確に異なる物理的チャネルは絡み合い、消失してしまう。
2. サブグリッドモデルの限界： 一般的に使用される粗い解像度のサブグリッドモデルのパラメータ（SN冷却半径やSF密度閾値など）は、再評価が必要である可能性がある。これらのモデルは、セル内のISMが均質であることを前提としているが、本論文はその仮定がSNの結果を予測する上で重大な誤差を招くことを示している。
3. 代替的なモデリング戦略： 必要な解像度を達成できないシミュレーションの場合、高解像度の結果を明示的に解像するのではなく、その「創発的な効果」を模倣する有効なモデルが必要である。著者らは、「デカップルド・ウィンド（decoupled wind）」モデル（流出を低密度ガスに限定するもの）やマルチフルイド（多流体）アプローチが、高密度塊のジーンズ長を明示的に解像することなく、これらの効果を捉えるための有望な手段であると強調している。

本論文は、低解像度のシミュレーションが無用であると主張しているのではなく、それらが粗い解像度の情報のみに基づいて、特定のSNイベントの具体的なフィードバックチャネル（流出か抑制か）を自己整合的に決定することは不可能であると主張している。