Papaer II
出典：Kim et al. 2016, ApJ, 833, 202, doi:10.3847/1538-4357/833/2/202）。
📘 AGORA Project Paper II の目的と方法
▶ 対象：
- Milky Way質量の孤立円盤銀河（非宇宙論的）
- ガス＋ダークマター＋星成分を持つ安定した初期条件
▶ 参加コード（9種）：
- SPH系：CHANGA, GASOLINE, GEAR, GADGET-3
- AMR系：ART-I, ENZO, RAMSES
- メッシュレス系：GIZMO-MFM（Godunov型）
- その他：GEAR（化学進化付き）
▶ 共通設定：
- 同一の初期条件ファイル
- 物理モデル共通化：
- 冷却・加熱：GRACKLE
- 星形成：統一スレッショルド・効率（
n_thresh = 10 cm⁻³,ε_SF = 1%） - フィードバック：Type II超新星ごとに \(10^{51}\) ergを注入
- 星形成にはJeans pressure floorあり（分解能に応じた最低温度を設定）
🔬 主要な結果（具体的に）
1. 大局的構造（ガス/星の分布、回転曲線）
- 非常によく一致
- 図5・6：ガスと星の表面密度プロファイル → 各コードとも指数的円盤を形成
- 図8：回転曲線 → 8コードで重ねて見えるほど一致
- 回転速度のピーク位置・値ともに一致
2. 速度分散
- 図9：垂直方向のガス速度分散（σ_z）
- 半径全域で±20%程度の差（極めて良好な一致）
3. 星形成率とクランプ
- 図10：SFR vs 時間 → 全体傾向は共通（時間とともに緩やかに増加）、ピーク値の時刻にややズレ
- 図11：新形成星のクランプ質量関数
- 最大で約3倍の差（特に高質量側）
- メッシュコードはクランプ形成がやや抑制される傾向あり（拡散性の違い）
4. 密度・温度のPDF（Probability Distribution Function）
- 図14：Gas density PDF（質量加重）
- 低密度ガス（n_H ≲ 0.1 cm⁻³）ではメッシュ vs 粒子でばらつき
- 高密度ガス（n_H ≳ 100 cm⁻³）では、スパイク的構造の違いが見られ、粒子系がより濃密なクランプを形成
- 温度PDFも類似傾向
⚠️ システマティックなコード差があった箇所
| 特徴 | 原因と見られる要因 |
|---|---|
| 高密度ガスの“しっぽ”のばらつき | 数値拡散の差（AMRはdiffusive） |
| 低密度ガスでの挙動の違い | 解法の本質的性質（粒子 vs メッシュ） |
| 星形成クランプ数・質量分布 | セル内 vs 粒子間分解能の差や時間積分アルゴリズムの差 |
🧠 論文の結論
- 共通の初期条件・物理設定・解析手法を導入すれば、最新のコード間では銀河スケールの基本物理に関して非常に良い一致が得られる
- 違いが出る部分は、コードの「数値解法」そのものではなく、サブグリッド物理や拡散スキームの違いによる
- 高密度・低密度・速度場など、異なるスケールと物理過程でコード特性が現れる
✅ 要点まとめ（箇条書き）
- 同じ設定下では、主要な銀河構造はどのコードでもよく再現
- 数値解法の違いは、星形成クランプやガスPDFの端部で主に現れる
- フィードバックの扱い方や数値的拡散が、最もコード差の原因となる
- 結果的に、「コードを問わず信頼できる」銀河形成シミュレーションの可能性を示した