Forecast for the detectability of patchy hydrogen reionization in WEAVE-QSO measurements of the Lyman-$\alpha$ forest power spectrum at redshift $z \geq 4$

기술 요약: WEAVE-QSO 관측을 통한 불균일한 수소 재이온화의 검출 가능성 예측

문제 정의
성간 매질(IGM)은 수소 재이온화 시대의 "열적 유물(thermal relics)"을 간직하고 있다. 재이온화 과정 동안, 광학적으로 가열된 가스는 이온화된 거품(bubbles) 내에서 단열 냉각을 겪으며 장기간에 걸쳐 IGM 온도 분포에 공간적 불균질성을 각인시켰다. 이러한 불균질성은 적색편이 $z \sim 4$까지 지속되며, 이는 라이먼-알파($Ly\alpha$) 숲의 플럭스 파워 스펙트럼($P_{Ly\alpha}$)의 대규모(낮은 파수, $k \lesssim 10^{-3} \, \mathrm{s\,km^{-1}}$) 성분을 강화할 것으로 예측된다. 고해해상도 데이터(예: KODIAK, SQUAD, XQ-100)를 사용한 기존 연구들은 이와 같은 강화 현상에 대해 미약한 증거($\simeq 2.7\sigma$)를 발견했으나, 제한된 통계적 정밀도와 적색편이 경로 길이로 인해 높은 유의성을 가진 검출에는 도달하지 못했다. 저자들은 방대한 샘플 크기를 갖춘 차세대 WEAVE-QSO 서베이가 이 대규모 유물 흔적을 결정적으로 검출하고 특성화할 수 있는지 예측하고자 한다.

방법론
저자들은 WEAVE-QSO 서베이 구성에 맞춤화된 모의 $Ly\alpha$ 숲 스펙트럼을 생성하고 관련 불확실성을 예측하기 위해 종합적인 파이프라인을 개발하였다.

1. 유체역학 시뮬레이션: IGM의 기저 구조와 열적 역사는 Sherwood-Relics 시뮬레이션 제품군을 사용하여 모델링되었다. 이 우주론적 유체역학 시뮬레이션은 재이온화 과정과 그 이후의 IGM을 정밀하게 분해하며, 예측된 대규모 파워 강화를 생성하는 불균일한 재이온화 효과를 포함한다.
2. 모의 스펙트럼 생성 (WEAVEify): 현실적인 WEAVE-QSO 스펙트럼으로 Sherwood-Relics 출력을 순방향 모델링하기 위해 전용 소프트웨어 패키지인 WEAVEify가 개발되었다. 이 과정에는 다음이 포함된다:
   • 광학적 두께 스큐어(skewers)를 QSO 정지 좌표계(rest frames)로 매핑.
   • 기기 선형 함수(Instrumental Line Spread Function, $R \sim 5000$ 가정)를 시뮬레이션하기 위한 파장 의존적 가우시안 커널과의 컨볼루션.
   • SDSS QSO로부터 유도된 주성분 분석(PCA) 방법을 이용한 연속체(continuum) 배치.
   • 포아송 통계를 고려하여 검출기 특성, 하늘 밝기, 노출 시간($3060$ s)에 기반한 현실적인 광자 제한 노이즈(암전류 및 판독 노이즈 포함) 추가.
3. 서베이 샘플링: 저자들은 $z = 4.0, 4.2, 4.4, 4.6$을 중심으로 하는 4개의 적색편이 빈(bin)에 대해 $10,000$개의 모의 QSO 스펙트럼을 생성하였다. 샘플은 QSO 광도 함수에 기반한 몬테카를로 방법을 사용하여 $m_r < 21.9$의 $r$-밴드 등급 컷을 적용하여 구성되었다. 이는 $10,000 \, \mathrm{deg}^2$ 서베이 면적에 대해 각 빈당 각각 $2,942, 1,822, 1,385, 1,155$개의 QSO를 기대할 수 있음을 의미한다.
4. 파워 스펙트럼 추정 및 불확실성 모델링: 1D 파워 스펙트럼은 고속 푸리에 변환(FFT) 추정기를 사용하여 추정되었다. 저자들은 다음을 고려한 완전한 공분산 행렬을 구축하였다:
   • 통계적 불확실성: 모의 샘플 크기에 대한 부트스트래핑을 통해 도출.
   • 계통적 불확실성: 분해능 보정, 노이즈 제거, 상관관계가 없는 금속선 오염(C IV, Si IV, Mg II), 고밀도(HCD) 흡수체 불완전성, 연속체 배치 오류를 별도로 모델링.
   • 분석은 분해능 보정 계수가 $0.2$이상으로 유지되는 파수$k < 0.05 , \mathrm{s,km^{-1}}$ 영역으로 제한되었다.
5. 매개변수 추론: 모의 데이터와 이론적 모델을 비교하기 위해 베이지안 마르코프 체인 몬테카를로(MCMC) 프레임워크가 채택되었다. 모델에는 표준 IGM 매개변수(유효 광학적 두께 $\tau_{\mathrm{eff}}$, 온도 $T_0$, 기울기 $\gamma$, 열 에너지 $u_0$)와 불균일한 재이온화 보정을 나타내는 추가 진폭 매개변수 $A_p$가 포함되었다. 두 가지 시나리오가 테스트되었다: $A_p$를 가변적으로 허용하는 "표준(fiducial)" 적합과 불균일한 보정을 비활성화한 "No-$A_p$" 적합이다.

주요 기여

• 최초의 상세한 WEAVE-QSO 예측: 본 연구는 $z \ge 4$에서 WEAVE-QSO 서베이에 의한 $P_{Ly\alpha}$ 측정 정밀도에 대한 최초의 구체적인 예측을 제공한다.
• 포괄적인 계통 오차 모델링: 저해상도 대규모 서베이에 특화된 계통 불확실성, 즉 HCD 댐핑 윙(damping wings)과 연속체 배치 오류가 대규모 파워에 미치는 영향을 엄격하게 다루었다.
• 검출 능력의 검증: Sherwood-Relics 물리량을 실제적인 기기 모델링과 통합함으로써, 불균일한 재이온화의 특정 신호가 표준 IGM 열적 변화와 구별됨을 입증하였다.

결과

• 불확실성 예산: 예측된 $P_{Ly\alpha}$의 총 상대 불확실성은 $k_{\mathrm{max}}$ 이하의 모든 스케일과 적색편이 빈에 대해 평균했을 때 약 $5.2\%$이다. 대규모($k \lesssim 0.01 \, \mathrm{s\,km^{-1}}$)에서는 중간 범위의 통계적 오차와 매우 큰 규모에서의 계통 오차(특히 HCD 및 연속체 배치)에 의해 지배되지만, 불확실성은 $10\%$ 미만으로 유지된다.
• 검출 유의성: 불균일한 재이온화 신호를 포함하는 모의 데이터에 MCMC 분석을 적용했을 때:
  • $A_p$를 포함하는 표준 모델은 $A_p$가 없는 모델보다 유의미하게 더 나은 적합도를 보였다 ($\chi^2/\mathrm{dof}$가 $65.28/32$에서 $36.07/28$로 개선).
  • 이 개선은 4개의 추가 매개변수에 대해 $\Delta\chi^2 = 29.21$에 해당하며, 이는 양측 가우시안 등가 유의성으로 $\simeq 4.5\sigma$에 달한다.
  • $A_p$의 사후 분포(posterior distributions)는 모든 적색편이 빈에서 입력된 진폭을 성공적으로 회복하였다.
• 스케일 민감도: 검출 유의성은 가장 큰 스케일에 크게 의존한다. 가장 낮은-$k$ 빈($k \sim 10^{-3} \, \mathrm{s\,km^{-1}}$)을 제거하면 유의성이 $\simeq 1.2\sigma$로 떨어지며, 이는 검출이 열적 유물 흔적의 영향을 가장 많이 받는 모드들에 의존함을 나타낸다.

의의 및 주장
본 논문은 WEAVE-QSO 서베이가 예측된 샘플 크기와 공분산 특성을 갖추고 있다면, $z \ge 4$에서 라이먼-알파 숲 파워 스펙트럼의 대규모 유물 흔적을 높은 통계적 유의성($\simeq 4.5\sigma$)으로 직접 검출할 수 있다고 주장한다. 이는 소규모 고해상도 데이터셋에서 나타난 기존의 미약한 검출($\simeq 2.7\sigma$)보다 상당한 진전이다. 저자들은 개발된 프레임워크가 대규모 파워 스펙트럼 측정을 통해 재이온화 시기를 평가할 수 있게 해주며, WEAVE-QSO 구성이 계통 불확실성이 예측된 수준으로 제어된다면 미약한 증거를 확실한 검출로 전환하기에 충분하다고 결론짓는다. 이 연구는 재이온화의 정확한 적색편이를 측정한다고 주장하는 것이 아니라, 균일한 UV 배경 모델 대비 불균일한 흔적의 진폭을 측정하는 것을 목표로 한다.