Robust Ly$\alpha$ forest constraints on reionization from semi-analytic galaxy formation models

기술 요약: 반-해석적 은하 형성 모델을 통한 Ly$\alpha$ 포레스트의 강건한 재이온화 제약

문제 정의
재이온화 시대(EoR)는 우주의 바리온 물질이 겪는 마지막 주요 상전이를 나타낸다. 이전 연구(구체적으로 Qin et al. 2025, 이하 Q25)는 고적색편이 Ly$\alpha$ 포레스트가 재이온화의 시기와 형태에 대한 엄격한 제약을 제공한다는 것을 입증했으나, 이러한 추론의 강건성(robustness)에 대해서는 중요한 의문이 남아 있다. 기존의 EoR 추론 프레임워크는 은하의 특성과 헤일로 질량 사이의 단순화된 스케일링 관계에 의존하는 경우가 많으며, 이는 은하 형성의 복잡성을 결정론적인 매핑으로 압축한다. 고적색편이 은하는 별 질량, 별 형성사, 그리고 탈출 분율(escape fraction) 측면에서 상당한 확률론적 변동성(stochasticity)을 보일 것으로 예상된다. 본 논문은 확률론적 변동성, 활동 은하 핵(AGN), 그리고 은하 주변 매질(CGM) 감쇄를 포함하는 더 물리적으로 유연한 반-해석적 모델을 사용하여 테스트했을 때, Ly$\alpha$ 포레스트 데이터로부터 도출된 타이트한 후기 재이온화 역사가 여전히 유효한지 조사한다.

방법론
저자들은 $(210 h^{-1} \text{cMpc})^3$의 우주론적 $N$-body 시뮬레이션 부피에 적용된 은하 형성과 재이온화의 반-해석적 모델인 Meraxes를 활용한다. Meraxes는 가스 유입, 냉각, 별 형성, 초신성 피드백, 블랙홀 성장을 포함한 은하의 조립 역사를 추적한다. 재이온화 예산(ionizing budget)을 지배할 것으로 예상되는 저질량 은하 집단을 해상하기 위해, 몬테카를로 알고리즘을 사용하여 원자 냉각 임계값(atomic-cooling threshold)까지 머거 트리(merger trees)를 확장하였다.

본 연구는 은하 이온화 탈출 분율($f_{\text{esc}}^{\text{gal}}$)이 서로 다르게 매개변수화된 일련의 소스 모델들을 구축한다:

1. 표준 모델(Fiducial Model): $f_{\text{esc}}^{\text{gal}}$이 헤일로 질량($M_{\text{vir}}$)의 거듭제곱 법칙(power law)에 따라 변화함.
2. AGN 보조 모델(AGN-Assisted Model): 관측된 희미한 AGN 집단에 맞춰 보정된, 강착 블랙홀로부터 나오는 이온화 광자를 추가함.
3. CGM 감쇄 모델(CGM-Attenuated Model): 이온화 광자가 성간 매질(IGM)에 도달하기 전 밀집된 가스 구조에 의해 흡수되는 현상을 현상학적으로 모델링하여, 국소 클럼핑 인자($C$)와 CGM 칼럼 밀도에 기반한 추가적인 탈출 광자 억제를 도입함.
4. 은하 특성 의존 모델(Galaxy-Property-Dependent Models): $f_{\text{esc}}^{\text{gal}}$을 별 질량($M_*$), 별 형성률(SFR), 비정상 별 형성률(sSFR), 그리고 별 형성 디스크 표면 밀도($\Sigma_{\text{SF}}$)의 함수로 매개변수화함.

모든 모델은 XQR-30+ 데이터셋(Bosman et al. 2022)으로부터 관측된 Ly$\alpha$ 유효 광학적 두께($\tau_{\text{eff}}$)의 누적 분포에 맞춰 보정된다. 이온화장은 미해결된 재결합과 자기 차폐(self-shielding)를 고려하기 위해 하위 격자 밀도 분포를 포함하는 커스텀 버전의 21cmFAST를 사용하여 계산된다. Ly$\alpha$ 포레스트 투과율은 고해상도 유체역학 시뮬레이션에 의해 보정된 확률론적 매핑을 통해 변동하는 건-피터슨 근사(fluctuating Gunn-Peterson approximation)를 사용하여 순방향 모델링(forward-modeled)된다.

주요 결과
이온화 광자가 헤일로 질량, 은하 광도, 그리고 소스 유형(별 vs. AGN) 사이에 어떻게 배분되는지에 대한 상당한 변화에도 불구하고, 보정된 모든 유효 모델은 유사한 후기 재이온화 역사로 수렴한다:

• 재이온화 타이밍: 재이온화 중간 지점($z|_{x_{\text{HI}}=0.5}$)은 $z \sim 6.2$와 $6.9$사이**에서 발생한다. 최종 단계($z|{x{\text{HI}}=0.05}$)는 **$z \sim 5.36$과$5.57$ 사이에서 발생한다.
• 소스 재분배:
  • 표준 모델과 AGN 보조 모델은 탈출 광자 예산을 저질량 헤일로 쪽으로 이동시킨다. AGN 모델의 경우, AGN은 $z \sim 5$에서 이온화 방출량의 최대 $\sim 15\%$를 기여하며, 이는 후기에 IGM을 과소 이온화할 수 있는 가파른 별 탈출 분율을 보완한다.
  • CGM 감쇄 모델은 질량이 큰 헤일로로부터의 더 큰 고유 기여를 허용하지만, 관측된 불투명도와 일치시키기 위해 강력한 후기 억제( $C^3$에 비례)를 필요로 하며, 이는 재이온화 말기에 탈출하는 별의 방출량을 $\sim 15\%$ 효과적으로 감소시킨다.
  • 은하 특성 모델: $M_*$ 및 SFR에 연결된 모델들은 헤일로 질량과의 강한 상관관계로 인해 표준 모델과 유사하게 행동한다. sSFR 모델은 상당한 산포를 도입하지만 유사한 질적 계층 구조를 유지한다. SF-디스크 모델은 가장 질량이 큰 헤일로에 가중치가 놓인 시나리오이지만, 보정 후에도 동일한 후기 역사로 수렴하며, 이를 위해 추가적인 CGM 감쇄가 필요하다.
• 클럼핑 인자: 표준 모델은 재이온화 종료 시점에 $\sim 10$까지 상승하는 클럼핑 인자($C$)를 예측하며, 이는 IGM을 과잉 이온화하지 않으면서 높은 이온화 방출량을 균형 있게 맞춘다.
• CMB 광학적 두께: 보정된 모델들의 톰슨 산란 광학적 두께($\tau_e$)는 0.0535에서 0.0652 사이이며, 이는 Planck 제약 조건과 일치한다.

의의 및 주장
본 논문은 Ly$\alpha$ 포레스트 데이터로부터 추론된 늦고 확장된 재이온화 역사가 은하-헤일로 연결, 완만한 AGN 기여, 그리고 현상학적 CGM 감쇄와 같은 물리적으로 동기화된 변동에 대해 **강건함(robust)**을 주장한다.

저자들은 Ly$\alpha$ 포레스트가 재이온화의 마지막 단계에 대한 타이밍과 진행 과정을 강력하게 제한하지만, 이온화 광자 예산의 물리적 기원(즉, 저질량 은하, 질량이 큰 은하, AGN, 그리고 흡수 사이의 구체적인 균형)을 유일하게 결정하지는 않는다고 결 결론짓는다. 보정된 Meraxes 모델들 사이의 편차는 소스 모델링에서 발생하는 EoR 추론의 계통 오차(systematic uncertainty)에 대한 추정치를 제공한다. 이 연구는 단순한 헤일로 기반 처리를 넘어 더 상세하고 확률론적인 반-해석적 프레임워크로 이동하더라도, 이전 연구(Q25)에서 도출된 후기 EoR에 대한 엄격한 제약이 유효함을 재확인한다.