[論文レビュー] The Impact of the WHIM on the IGM Thermal State Determined from the Low-$z$ Lyman-$α$ Forest

Teng Hu, Vikram Khaire+5 more|arXiv (Cornell University)|Aug 28, 2023

Galaxies: Formation, Evolution, PhenomenaPhysics and Astronomy被引用数 3

ひとこと要約

本研究は、低赤方偏移における銀河間物質（IGM）の熱的状態が、$T_0$ および $\gamma$ でパラメータ化されているが、温熱銀河間物質（WHIM）による大きな衝撃加熱があっても、Ly$\alpha$森から信頼性高く測定可能であることを示している。Nyx、IllustrisTNG、Illustrisのシミュレーションを用いた機械学習を活用したシミュレーションベースの推論により、WHIMの割合の違いに対して $b$--$N_{\text{H\,{\tiny i}}}$ 分布がほとんど感度を示さないことが明らかになった。その結果、$|\Delta\log(T_0/\text{K})|\lesssim 0.05$ デキシタル、$|\Delta\gamma|\lesssim 0.1$ のバイアスが得られ、これは通常の観測誤差よりも小さい。

ABSTRACT

At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, $T=T_0 (ρ/ {\barρ})^{γ-1}$. This study explores the impact of the WHIM on measurements of the low-$z$ IGM thermal state, $[T_0,γ]$, based on the $b$-$N_{H I}$ distribution of the Lyman-$α$ forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [$T_0$, $γ$] can still be reliably measured from the $b$-$N_{H I}$ distribution at $z=0.1$, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at $z=0.1$. The results suggest that the underlying $[T_0,γ]$ of both simulations can be recovered with biases as low as $|Δ\log(T_0/ ext{K})| \lesssim 0.05$ dex, $|Δγ| \lesssim 0.1$, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the $b$-$N_{H I}$ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-$α$ absorbers, and discover that although their $T$ and $Δ$ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

研究の動機と目的

• 低赤方偏移 ($z \lesssim 1$) における温熱銀河間物質（WHIM）の存在が、Ly$\alpha$森からのIGM熱的状態 $[T_0, \gamma]$ 測定に歪みを及えるかどうかを評価すること。
• さまざまなフィードバックモデル下でのIGM密度-温度関係のプローブとしての $b$--$N_{\text{H\,{\tiny i}}}$ 分布の頑健性を評価すること。
• 宇宙論的シミュレーション間のWHIM割合の違いが、観測されたLy$\alpha$森統計から推定される熱的パラメータに影響を及ぼすかどうかを特定すること。
• 異なるフィードバック処置を有するシミュレーションにおいて、$[T_0, \gamma]$ をモックLy$\alpha$森スペクトルから回復する機械学習ベースの推論フレームワークの妥当性を検証すること。

提案手法

• 高解像度のNyx流体力学的シミュレーションを用いてトレーニングされた機械学習を活用したシミュレーションベースの推論手法を採用し、$z = 0.1$ における $b$--$N_{\text{H\,{\tiny i}}}$ 分布から $[T_0, \gamma]$ を推定した。
• 同様の推論フレームワークを、フィードバック強度およびそれに伴うWHIM割合に顕著な差があるIllustrisTNGおよびIllustrisシミュレーションから導出されたモックLy$\alpha$森スペクトルに適用した。
• 密度推定に基づく尤度フリー推論（DELFI）を用い、観測された $b$--$N_{\text{H\,{\tiny i}}}$ 分布を $[T_0, \gamma]$ に関する後確率分布にマッピングすることで、複雑な非線形性に対してもパラメータ回復が可能となった。
• シミュレーションを写し合わせ、同じ光電離率 $\Gamma_{\text{H\,{\tiny i}}}$ を共有することで、フィードバックの影響とイオン化歴の影響を分離した。
• シミュレーション間でのLy$\alpha$吸収体の $T$--$\Delta$ 分布を比較し、フィードバックおよびWHIM含有量への感受性を評価した。
• 各シミュレーションにおける真の入力値と回復された $[T_0, \gamma]$ を比較することで、バイアスと精度を定量的に評価し、結果を検証した。

実験結果

リサーチクエスチョン

• RQ1 低赤方偏移 ($z \lesssim 1$) におけるWHIMの存在が、Ly$\alpha$森からのIGM熱的状態 $[T_0, \gamma]$ 測定に顕著なバイアスを及えるか？
• RQ2 IGMに顕著なWHIM（$T > 10^5$ K）が存在する場合でも、Ly$\alpha$森における $b$--$N_{\text{H\,{\tiny i}}}$ 分布が $[T_0, \gamma]$ を信頼性高く回復できるか？
• RQ3 異なるフィードバックモデル（例：Illustris と IllustrisTNG）が、熱的状態パラメータ $[T_0, \gamma]$ およびLy$\alpha$森におけるWHIMの検出可能性にどのように影響を与えるか？
• RQ4 観測可能なLy$\alpha$吸収体の物理的性質が、フィードバック機構に依存しない形で光電離率とどの程度相関するか？
• RQ5 光電離率を一定に保った状態で、$b$--$N_{\text{H\,{\tiny i}}}$ 分布が体積平均WHIM割合に感度を示すか？

主な発見

• 機械学習ベースの推論フレームワークは、Nyx、IllustrisTNG、Illustrisの3つのシミュレーションすべてにおいて真の $[T_0, \gamma]$ パラメータを的確に回復した。バイアスは $|\Delta\log(T_0/\text{K})| \lesssim 0.05$ デキシタル、$|\Delta\gamma| \lesssim 0.1$ であり、これは通常の測定誤差よりも小さい。
• Nyxでは4%、IllustrisTNGでは10%、Illustrisでは38%と顕著な体積平均WHIM割合の差があるにもかかわらず、$b$--$N_{\text{H\,{\tiny i}}}$ 分布は現実的な条件下でこれらの変動に対してほとんど感度を示さなかった。
• 光電離率 $\Gamma_{\text{H\,{\tiny i}}}$ を一定に保った場合、シミュレーション間でLy$\alpha$吸収体の $T$--$\Delta$ 分布はほとんど区別がつかず、フィードバック効果が観測された森の統計に隠蔽されていることが示された。
• IllustrisTNGとIllustrisの両方で、Ly$\alpha$吸収体内のWHIM割合は平均11.6%とほぼ同一であり、IGM全体のWHIM割合に顕著な差があるにもかかわらず、これは吸収体がフィードバック強度にかかわらず同じ物理的状態を優先的にトレースしていることを示唆している。
• $b$--$N_{\text{H\,{\tiny i}}}$ 分布は現実的な条件下でWHIMに感度を示さないため、衝撃加熱が顕著であっても、標準的なIGM熱的状態測定がLy$\alpha$森から堅牢に可能であることが示された。
• 本研究は、IGMが衝撃加熱およびフィードバックによる大きな分散を示しても、$T_0$--$\gamma$ 関係が $z \lesssim 1$ において依然として有効かつ測定可能であることを確認した。