[論文レビュー] Cosmological baryon spread and impact on matter clustering in CAMELS

Matthew Gebhardt, Daniel Anglés‐Alcázar+10 more|arXiv (Cornell University)|Jul 21, 2023

Cosmology and Gravitation Theories被引用数 4

ひとこと要約

本稿では、CAMELSプロジェクトの広大なシミュレーション・スイートを用いて、超新星およびAGNからのフィードバックが、初期のダークマター分布に対してバリオンをどのように再配分するかを定量化するための宇宙論的バリオン拡散度を導入する。この度は、バリオン拡散が物質パワースペクトルの抑制と強く相関しており、特にSIMBAシミュレーションでは効率的なAGNフィードバックのおかげで約40%のバリオンが1 Mpc以上離れて分布している一方、IllustrisTNGやASTRIDではわずか約10%にとどまる。また、シンボリック回帰により、多様なフィードバックモデルにわたるパワー抑制を予測する解析的関数が同定された。

ABSTRACT

We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighboring dark matter owing to feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial extsc{SIMBA} model spreading $\sim$40\% of baryons $>$1\,Mpc away compared to $\sim$10\% for the IllustrisTNG and extsc{ASTRID} models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired $N$-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number ($k$) and baryonic spread up to $k \sim 10\,h$\,Mpc$^{-1}$ while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation.

研究の動機と目的

• 宇宙論的ハイドロダイナミクス・シミュレーションにおける初期ダークマター近傍粒子と比較したバリオンの宇宙論的拡散を定量化すること。
• 特にAGNおよび星間フィードバックが大規模なバリオン再配分に与える影響と、物質クラスタリングへの影響を評価すること。
• 多様な銀河形成モデルや宇宙論的パラメータにわたって、物質パワースペクトルのグローバルな抑制を予測するための、バリオン拡散が強力で普遍的な指標であるかどうかを評価すること。
• 訓練パラメータ空間を超えて外挿可能な、バリオン拡散と物質パワー抑制の関係を捉えるシンボリック回帰モデルの構築と検証すること。

提案手法

• 本研究では、CAMELSプロジェクトの3,000以上の高解像度ハイドロダイナミクス・シミュレーションを用い、同じ初期条件を持つフルハイドロダイナミクスとペアリングされたN体シミュレーションを含む。
• 新たな「バリオン拡散」度は、各バリオンとその初期に最も近いダークマター粒子との間の径方向距離を、ハローのラグランジュ領域全体で平均化することで定義される。
• ハイドロダイナミクスとN体シミュレーションの両方で物質パワースペクトルを計算し、フィードバックに起因する大規模構造クラスタリングの抑制を定量化する。
• 波数（k）、バリオン拡散、物質パワー抑制の関係を記述する閉形式の解析的関数を学習するため、シンボリック回帰を適用し、4つのフィードバックパラメータを変化させたSIMBAシミュレーションで訓練する。
• 訓練済みモデルの予測能力を、異なるフィードバック実装とパラメータセットを有する他の銀河形成モデル（IllustrisTNG、ASTRID）のシミュレーションでテストする。
• 非線形領域における物質パワースペクトルの予測を目的として、バリオン拡散と宇宙論的パラメータを用いたランダムフォレスト回帰により、サーヴェイモデルを訓練する。

実験結果

リサーチクエスチョン

• RQ1 宇宙論的ハイドロダイナミクス・シミュレーションにおける、異なる銀河形成モデルにおける初期ダークマター近傍粒子と比較したバリオン拡散は、どのように変動するか？
• RQ2 AGNフィードバックは星間フィードバックと比較して、どれほど大規模なバリオン再配分を促進するか？
• RQ3 バリオン拡散は、多様なフィードバック処方および宇宙論的パラメータにわたって、物質パワースペクトルの抑制を普遍的に予測する指標として機能するか？
• RQ4 SIMBAモデルで訓練されたシンボリック回帰モデルは、異なるフィードバック実装（IllustrisTNG、ASTRID）を有するシミュレーションにおけるパワー抑制をどれほど正確に予測できるか？

主な発見

• 重力的力学的不規則性のため、ダークマターは拡散を示し、質量の大きなハローでより大きな拡散を示す。また、Ωmやσ8が高くなるほど拡散が増加する。
• フィードバックの影響により、ガスはダークマターと比較してはるかに遠くまで拡散する。特にSIMBAシミュレーションでは、約40%のバリオンが1 Mpc以上離れて分布しているが、IllustrisTNGやASTRIDではわずか約10%にとどまる。
• AGNフィードバックの効率を高めるとバリオン拡散が増加するが、強力な星間フィードバックは非線形結合のためブラックホール成長とAGN活動を抑制するため、逆に拡散を減少させる可能性がある。
• バリオン拡散度は物質パワースペクトルの抑制と強く相関しており、フィードバックが大規模構造に与える影響を示す強力なグローバル指標である。
• シンボリック回帰により、SIMBAシミュレーションの20以上の宇宙物理学的パラメータにわたるkおよびバリオン拡散の関数としてパワー抑制を捉えるシンプルな解析的関数が得られたが、拡散度が低いモデル（例：IllustrisTNG や ASTRID）には一般化できない。
• 本研究では、銀河形成物理学の実装の違いですら、要約統計量（例：バリオン拡散）を用いても、フィードバックモデルの堅牢性を確保することが難しいことが浮き彫りにされた。