A Measurement of the Thermal and Ionization State of the IGM at $z < 0.5$

기술 요약: $z < 0.5$에서의 IGM의 열적 및 이온화 상태 측정

문제 제기
저적색($z < 0.5$)에서의 은하간 매질(IGM)의 열적 및 이온화 상태는 관측 데이터와 이론적 예측 사이에 상당한 불일치를 보인다. 고적색 IGM 모델은 Ly$\alpha$ 포레스트 관측과 잘 일치하지만, 허블 우주 망원경(HST)의 코스믹 오리진 분광기(COS)를 통한 저적색 데이터는 예측보다 훨씬 넓은 Ly$\alpha$ 흡수선을 보여준다. He II 재이온화($z \approx 3$) 이후의 단열 냉각을 고려하는 표준 우주론 모델은 $z=0$에서 평균 IGM 온도($T_0$)가 $\sim 4000$ K이고 온도-밀도 기울기($\gamma$)가 $\sim 1.6$일 것으로 예측한다. 그러나 이전 연구들은 저적색 IGM이 이러한 정형적인 값들보다 더 뜨겁거나 더 난류적일 수 있음을 시사해 왔으나, 제한된 데이터 크기와 수소 광이온화율($\Gamma_{\text{HI}}$}$) 사이의 퇴행성으로 인해 정밀한 측정이 저해되어 왔다.

방법론
본 연구는 네 개의 적색 범위($z=0.1, 0.2, 0.3, 0.4$)에 걸쳐 IGM의 열적 상태($T_0, \gamma$)와 이온화 상태($\Gamma_{\text{HI}}$)를 측정하기 위해 머신러닝 기반의 추론 프레임워크를 사용한다.

1. 데이터 처리: 분석에는 높은 신호 대 잡음비(S/N)를 가진 82개의 아카이브 HST/COS 퀘이사 스펙트럼(Danforth et al. 2016)을 활용한다. 저자들은 결합 도플러 파라미터($b$)와 중성 수소 칼럼 밀도($N_{\text{HI}}$) 분포를 추출하기 위해 VPFIT 코드의 커스텀 래퍼를 사용하여 자동화된 보이트 프로파일 피팅(Voigt profile fitting)을 수행한다. 금속 흡수선과 근접 영역(proximity zones)은 Danforth et al. (2016)의 식별에 따라 마스킹되며, 최종 샘플은 657개의 Ly$\alpha$ 흡수선이다.
2. 시뮬레이션 및 순방향 모델링: 저자들은 다양한 열적 역사를 가진 51개의 Nyx 우주론적 유체역학 시뮬레이션 세트를 사용한다. 광이온 가열만으로는 저적색에서 높은 온도의 모델을 생성하기 어렵다는 점을 해결하기 위해, $T_0$가 $\sim 30,000$ K에 달하는 모델을 생성하도록 사후 처리 온도 재조정(post-processing temperature rescaling) 기법을 적용한다. 또한 $\Gamma_{\text{HI}}$를 13개의 값에 걸쳐 변화시킨다. 결정적으로, 순방향 모델링에는 COS 장비의 특정한 비가우시안 선 퍼짐 함수(LSF)가 포함되며, 관측적 인위적 현상에 맞추기 위해 노이즈 플로어를 추가한다.
3. 추론 엔진: 본 연구는 신경망 밀도 추정기를 사용하는 밀도 추정 무결론적 추론(DELIFI)을 활용한다. 신경망은 시뮬레이션된 모의 스펙트럼으로부터 조건부 확률 분포 $P(b, N_{\text{HI}} | T_0, \gamma, \Gamma_{\text{HI}})$를 학습하도록 훈련된다. 이 에뮬레이터는 흡수선 밀도($dN/dz$)에 대한 가우시안 프로세스 에뮬레이터와 결와되어, 관측된$b-N_{\text{HI}}$ 분포과 선 개수로부터 IGM 파라미터를 직접 베이지안 추론함으로써 열적 상태와 이온화율 사이의 퇴행성을 효과적으로 해소한다.

주요 결과
추론 결과는 $z=0.1$에서의 IGM에 대해 다음과 같은 측정값을 제공한다:

• 온도: $\log(T_0/\text{K}) = 4.45^{+0.08}_{-0.12}$이며, 이는 $T_0 = 28,183^{+5,700}_{-6,804}$ K에 해당한다. 이는 정형적인 예측치인 $\sim 4000$ K보다 약 7배 더 뜨겁다.
• 열적 기울기: $\gamma = 1.06^{+0.13}_{-0.09}$로, 저적색 IGM이 예상되는 $\gamma \sim 1.6$와는 극명하게 대조되는 거의 등온(isothermal) 상태임을 나타낸다.
• 이온화율: $\log(\Gamma_{\text{HI}}/\text{s}^{-1}) = -13.70^{+0.10}_{-0.08}$이다. 이 값은 현재의 UV 배경 합성 모델(예: Khaire & Srianand 2019)의 예측보다 약 $4\sigma$ 낮다.

이 결과들은 네 개의 적색 범위($z=0.1$에서 $0.4$) 전체에서 일관되게 나타나며, 표준 모델이 예측하는 것보다 지속적으로 뜨겁고 거의 등온인 IGM과 표준 모델보다 낮은 광이온화율을 보여준다.

대안적 설명 및 테스트
넓은 Ly$\alpha$ 흡수선의 기원을 조사하기 위해, 저자들은 두 가지 대안적 가설을 테스트했다:

1. 소규모 난류: 표준 열 모델($T_0 \sim 4000$ K)에 난류 속도 성분($v_{\text{tur}}$)을 도입함으로써, 저자들은 $z=0.1$에서 $v_{\text{tur}} \approx 15$ km s$^{-1}$의 난류 분산이 관측된 선 폭을 재현할 수 있음을 발견했다. 요구되는 난류는 저적색으로 갈수록 증가하며, 이는 에너지 주입의 지속적인 원천이 존재함을 시사한다.
2. 장비 해상도: 저자들은 HST COS 해상도가 과대평가되었는지 확인하기 위해 동일한 천체의 더 높은 해상도인 HST/STIS E140M 스펙트럼과 COS 스펙트럼을 비교하였으며, 이를 통해 보고된 COS LSF가 정확함을 확인했다. 그들은 표준 열 모델로 관측된 $b$-분포를 재현하려면 유효 해상도가 $R \sim 7000$ (FWHM $\sim 43.6$ km s$^{-1}$)이어야 함을 입증하였는데, 이는 알려진 장비 성능과 모순된다.

의의 및 주장
본 논문은 저적색 IGM이 광이온 가열 및 단열 냉각에만 의존하는 표준 우주론 모델의 예측보다 훨씬 더 뜨겁고 더 등온적이라고 주장한다. 저자들은 이 불일치가 장비 해상도 오류로 인해 발생할 수 없다고 단언한다. 결과적으로, 이 발견은 다음 중 하나를 시사한다:

• 새로운 가열 메커니즘: $z \sim 1$ 부근에서 중요해지고 현재까지 지속되는 암흑 물질 소멸, 암흑 광자 가열, 감마선 블레이자 또는 먼지 가열과 같은 과정.
• 미해결 난류: 은하 형성 피드백이나 현재의 시뮬레이션에 포착되지 않은 다른 과정에 의해 유도되는 소규모 난류 운동.

이 연구는 높은 온도와 낮은 광이온화율의 조합이 지난 100억 년 동안의 IGM의 열적 역사 전반에 대한 현재의 이해에 도전한다는 점을 강조한다. 저자들은 이 불일치를 해결하기 위해서는 새로운 물리적 가열원을 식별하거나 IGM 모델에 미해결 난류를 포함시켜야 한다고 결론짓는다.