25/05/20 進捗報告
1. 背景・問題提起
- シミュレーションデータからハロー質量ごとのSFR進化を比較したところ、最も質量の大きいハローでz=0付近でSFRが再び上昇している現象が見られた。
- 一般的な期待（高質量ハローではフィードバックにより星形成抑制が持続）とは逆の傾向。
- 「なぜ高質量ハローだけSFRが再活性化するのか？」を物理的に解明することが本研究の目的。
2. 初期考察・仮説
- 冷却流(cooling flow)問題や、冷却時間(cooling time)と自由落下時間(free-fall time)の関係が鍵ではないかと予想。
- さらには、「一度フィードバックで排出されたガス（recycled gas）」が中心に戻ることで、星形成が再度活発化している可能性も示唆された。
3. 解析手法・データ構成
使用データ
- ラジアルプロファイル（featherファイル形式）で
Density,Temperature,Metallicity,ElectronAbundance,SFRなど
- pygrackleを用いた冷却関数計算
- AGN/SN feedback関連量：
BlackHoleAccretionRate, エネルギー注入公式
冷却時間(t_cool)の計算
- \( t_{\text{cool}} = \frac{1.5 k_B T}{n_e \Lambda} \)
- 必要な物理量（T, density, metallicity, electron abundance）はプロファイルの中央値などを利用
冷却/加熱力比の評価
- \( \dot{E}_{\text{cool}} = n_e n_H \Lambda(T, Z) \cdot V \)
- \( \dot{E}_{\text{AGN}} = \epsilon_{\text{FB}} \epsilon_{\text{rad}} \dot{M}_{\text{acc}} c^2 \)
- \( \dot{E}_{\text{SN}} = \eta_{\text{SN}}(Z) \dot{M}_{\text{star}} \)
- 実際の解析では「冷却が起きそうな領域（中心またはt_cool/t_ff < 10）」の値を用いる
4. 主な解析・可視化結果
4.1 SFR進化・累積SFRプロファイル
- ハロー質量が大きいほど、z=0でのSFR回復が顕著。
- AGN feedbackの有無でSFRに大きな差が出る（特に高質量ハロー）。
4.2 冷却条件の空間分布
- \( t_{\text{cool}}/t_{\text{ff}} < 10 \) を満たす領域が高質量ハローの中心で顕著。
- 低・中質量ハローではこの条件を満たす領域が限定的。
4.3 金属量プロファイルとT-ρ-Zダイアグラム
- どのハロー質量でも中心のZは高いが、高質量ハローの中心で高Zかつ低温・高密度なガスが存在。
- これはrecycled gasの存在を示唆するが、Zのみでは確定的証拠にはなりにくい（T,ρの情報と合わせる必要あり）。
4.4 AGN/SNフィードバックの進化
- 高質量ハローではz=0に向かってAGN/SNエネルギー注入率が低下。
- 一方でBH accretion rate自体は大きく減っていない場合もあり、「feedback効率の低下」や「加熱が追いつかなくなっている」ことを示唆。
4.5 AGN有無による温度・密度・SFRプロファイル
- 低質量ハロー：AGNの有無に関わらずSFR抑制効果は小さい。
- 中質量ハロー：AGN feedbackが最も効率的にSFRと冷却を抑制。
- 高質量ハロー：AGNがあってもSFRが再燃、中心ガス密度も高く保たれている。
5. 考察・結論
- 冷却流の再発生（冷却の再活性化）が高質量ハローでSFRを再度駆動している。
- ただし冷却条件（t_cool/t_ff < 10）だけではSFR復活の説明としては不十分であり、
- (a) フィードバック（AGN/SN）が十分に弱まっている
- (b) 中心へのガス供給量が豊富
- (c) 冷却後のガスが星形成領域まで到達できている
- (d) enrichedなrecycled gasが中心で再冷却している など複合的な条件が同時に成立している必要がある。
- 中間・低質量ハローでは、冷却条件を満たしていても「供給量不足」や「フィードバックがまだ強い」などの理由でSFRが抑制されている。
6. 検証のための追加解析提案
- t_cool/t_ff < 10 かつSFRが低い領域での「冷却 vs フィードバック加熱」のエネルギー収支比較
- 冷却レートの計算では、物理的に意味のある代表値（中心・冷却流領域のmedian/proper weighted mean）を使う
- Black hole accretion region周辺の温度・密度・速度プロファイルの診断
- 時系列でrecycled gasの挙動・金属量変化をトラッキング
7. 今後の展望・方針
- エネルギー収支やガス流入/流出の時間発展を解析し、「星形成再燃の必要十分条件」を定量的に示す
- SFR回復と観測的な冷却流銀河団・AGN duty cycleの関連性を議論
- 粒子追跡やフェーズ分類（cold/warm/hot, enriched/un-enriched）など多角的な検証を進める