[논문 리뷰] The Impact of the WHIM on the IGM Thermal State Determined from the Low-$z$ Lyman-$α$ Forest

Teng Hu, Vikram Khaire외 5명|arXiv (Cornell University)|2023. 08. 28.

Galaxies: Formation, Evolution, PhenomenaPhysics and Astronomy인용 수 3

한 줄 요약

이 연구는 저 redshift 간성간 매질(IGM)의 열적 상태, 즉 $T_0$와 $\gamma$로 매개변수화된 상태가, 따뜻한-핫 간성간 매질(WHIM)에 의한 상당한 충격 예열에도 불구하고 라이아프 숲(Ly$\alpha$ forest)을 통해 신뢰성 있게 측정될 수 있음을 보여준다. Nyx, IllustrisTNG, Illustris 시뮬레이션을 대상으로 기계학습 기반 시뮬레이션 기반 추론을 수행한 결과, $b$--$N_{\text{H\,{\tiny i}}}$ 분포는 WHIM 비율의 차이에 대해 민감하지 않으며, 이로 인해 $|\Delta\log(T_0/\text{K})|\lesssim 0.05$ dex 및 $|\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.

연구 동기 및 목표

• 저 redshift ($z \lesssim 1$)에서 따뜻한-핫 간성간 매질(WHIM)의 존재가 라이아프 숲을 통한 IGM 열적 상태 $[T_0, \gamma]$ 측정에 왜곡을 유도하는지 평가하는 것.
• 다양한 피드백 모델 하에서 IGM 밀도-온도 관계의 추적 수단으로서의 $b$--$N_{\text{H\,{\tiny i}}}$ 분포의 강건성 평가.
• 우주론적 시뮬레이션 간 WHIM 비율의 차이가 관측된 라이아프 숲 통계로부터 추론되는 열적 매개변수에 영향을 미치는지 확인하는 것.
• 다른 피드백 설계를 가진 시뮬레이션에서 라이아프 숲 가짜 스펙트럼으로부터 $[T_0, \gamma]$를 복원하는 기계학습 기반 추론 프레임워크의 검증.

제안 방법

• 고해상도 Nyx 유체역학 시뮬레이션을 기반으로 훈련된 기계학습 기반 시뮬레이션 기반 추론 방법을 사용하여, $z = 0.1$에서 $b$--$N_{\text{H\,{\tiny i}}}$ 분포로부터 $[T_0, \gamma]$를 추론하였다.
• 이와 동일한 추론 프레임워크를 피드백 강도와 그로 인한 WHIM 비율이 상당히 다른 IllustrisTNG 및 Illustris 시뮬레이션에서 유도된 라이아프 숲 가짜 스펙트럼에 적용하였다.
• 관측된 $b$--$N_{\text{H\,{\tiny i}}}$ 분포를 $[T_0, \gamma]$에 대한 사후 분포로 매핑하기 위해 밀도 추정 기반의 가능도 없는 추론(DELFI)을 사용하여, 복잡한 비선형성에도 불구하고 매개변수 복원이 가능하도록 하였다.
• 피드백의 영향과 이온화력 역사에 의한 영향을 분리하기 위해 시뮬레이션을 동일한 광이온화율 $\Gamma_{\text{H\,{\tiny i}}}$로 후처리하였다.
• 다양한 피드백 모델 하에서의 시뮬레이션된 라이아프 숲 흡수체의 $T$--$\Delta$ 분포를 비교하여 피드백과 WHIM 함량에 대한 민감도를 평가하였다.
• 각 시뮬레이션의 진짜 입력 값과 복원된 $[T_0, \gamma]$를 비교하여 편향과 정밀도를 정량화함으로써 결과를 검증하였다.

실험 결과

연구 질문

• RQ1 저 redshift ($z \lesssim 1$)에서 WHIM의 존재가 라이아프 숲을 통한 IGM 열적 상태 $[T_0, \gamma]$ 측정에 유의미한 편향을 유도하는가?
• RQ2 IGM에 상당한 WHIM($T > 10^5$ K)이 포함되어 있어도, 라이아프 숲 내의 $b$--$N_{\text{H\,{\tiny i}}}$ 분포가 $[T_0, \gamma]$를 신뢰성 있게 복원할 수 있는가?
• RQ3 다른 피드백 모델(예: Illustris vs. IllustrisTNG)은 라이아프 숲 내에서의 IGM 열적 상태 매개변수 $[T_0, \gamma]$와 WHIM 탐지 가능성에 어떤 영향을 미치는가?
• RQ4 피드백 메커니즘과 무관하게, 감지 가능한 라이아프 숲 흡수체의 물리적 특성과 광이온화율 간의 상관관계는 어느 정도인가?
• RQ5 광이온화율을 일정하게 유지할 경우, $b$--$N_{\text{H\,{\tiny i}}}$ 분포는 체적 평균 WHIM 비율에 민감한가?

주요 결과

• 기계학습 기반 추론 프레임워크는 Nyx, IllustrisTNG, Illustris의 세 시뮬레이션 전반에서 진짜 $[T_0, \gamma]$ 매개변수를 성공적으로 복원하였으며, 편향이 $|\Delta\log(T_0/\text{K})| \lesssim 0.05$ dex 이하, $|\Delta\gamma| \lesssim 0.1$ 이하로 나타났다. 이는 일반적인 측정 불확도보다 작다.
• 체적 평균 WHIM 비율의 큰 차이—Nyx에서는 4%, IllustrisTNG에서는 10%, Illustris에서는 38%—에도 불구하고, 현실적인 조건 하에서 $b$--$N_{\text{H\,{\tiny i}}}$ 분포는 이러한 변화에 거의 민감하지 않다.
• 광이온화율 $\Gamma_{\text{H\,{\tiny i}}}$를 일정하게 유지할 경우, 다양한 시뮬레이션 간에 시뮬레이션된 라이아프 숲 흡수체의 $T$--$\Delta$ 분포는 거의 구분되지 않으며, 이는 피드백 효과가 관측된 숲 통계에서 가려진다는 것을 시사한다.
• IllustrisTNG와 Illustris 간의 WHIM 비율이 크게 다름에도 불구하고, 라이아프 숲 흡수체 내의 평균 WHIM 비율은 거의 동일한 11.6%로 나타났다. 이는 흡수체가 피드백 강도와 관계없이 동일한 물리적 조건을 선별적으로 추적한다는 것을 의미한다.
• 실제 조건 하에서 WHIM의 영향을 고려할 경우, $b$--$N_{\text{H\,{\tiny i}}}$ 분포는 민감하지 않으며, 이는 충격 예열이 상당히 존재하더라도 라이아프 숲을 통한 표준 열적 상태 측정이 여전히 강건하다는 것을 의미한다.
• 이 연구는 IGM 열적 상태의 $T_0$--$\gamma$ 관계가 충격 예열과 피드백으로 인한 상당한 산란이 존재하는 저 redshift($z \lesssim 1$)에서도 여전히 타당하고 측정 가능한 기술로 남아 있음을 확인한다.