ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

기술 요약: "ECHO21: 암흑기부터 재이온화 시기까지의 전역 21-cm 신호 모델링을 위한 도구"

문제 및 동기
암흑기와 우주 새벽(Cosmic Dawn)은 직접적인 관측 프로브의 부재로 인해 우주 연대기에서 이해가 부족한 시기이다. 중성 수소의 21-cm 초미세 전이선은 최초의 별, 은하간 매질(IGM), 그리고 우주론적 매개변수의 물리학을 밝혀낼 수 있는 잠재적인 탐사 도구를 제공한다. 21cmFAST나 BEORN과 같은 수치적 또는 준수치적(semi-numerical) 코드들은 정확하고 물리적으로 풍부한 예측을 제공하지만, 계산 속도가 느려 대규모 데이터셋으로부터 효율적인 천체물리학적 또는 우주론적 추론을 수행하는 데 한계가 있다. 반대로, ARES나 Zeus21과 같은 기존의 해석적(analytical) 코드들은 속도는 빠르지만, 특히 Lyman-$\alpha$ (Ly$\alpha$) 광자에 의한 가열 효과와 같은 특정 물리적 구성 요소를 결여하고 있다. 저자들은 Ly$\alpha$ 가열의 누락이 신호 모델링에 강한 편향(bias)을 초래할 수 있음을 확인하였다. 빠르고 충분히 현실적이며, 대규모 매개변수 공간을 다룰 수 있는 유연한 도구에 대한 필요성을 해결하기 위해, 저자들은 ECHO21(Exploring the Cosmos with Hydrogen Observation)을 소개한다.

방법론
ECHO21은 적색편이 $z \approx 1500$부터 현재까지의 전역 21-cm 신호를 모델링하기 위해 해석적 프레임워크를 활용하는 Python 패키지이다. 이 코드는 IGM을 벌크(bulk) IGM과 완전히 이온화된 H II 영역의 두 영역 모델로 구분하여 작동한다.

• 물리적 프레임워크: 코드는 가스 운동 온도($T_k$), 벌크 IGM 전자 분율($x_e$), 그리고 H II 영역의 부피 충전율($Q$)에 대한 결합된 미분 방정식을 푼다.

  • 열적 역사(Thermal History): $T_k$의 진화는 단열 냉각, 컴프턴 가열(Compton heating), X-선 가열, 그리고 결정적으로 Ly$\alpha$ 가열을 고려한다. 저자들은 이전의 해석적 코드들이 Ly$\alpha$ 가열을 무시했던 반면, ECHO21은 이를 포함하며, 별 형성율 밀도(SFRD)로부터 유도된 Ly$\alpha$ 비특이 강도($J_{Ly}$)를 기반으로 가열률을 계산한다고 명시하였다.
  • 이온화 역사(Ionization History): 코드는 X-선 및 UV 광자에 의해 구동되는 이온화와 벌크 IGM의 재결합을 추적한다. 또한 복사 재결합에 대한 이온화 광자의 방출률을 비교함으로써 H II 영역의 부피 충전율을 모델링한다.
  • 별 형성율 밀도(SFRD): ECHO21의 참신한 특징은 세 가지 뚜렷한 SFRD 모델을 포함한다는 점이다:
    1. 물리적 동기 모델: 일정한 별 형성 효율을 사용하여 암흑 물질 헤일로로의 바리온 붕괴(Press-Schechter, Sheth-Tormen 또는 Tinker 헤일로 질량 함수 사용)를 기반으로 한다.
    2. 준경험적(Semi-empirical) 모델: 특성적인 별 형성 시간 척도와 최소 비릴 온도에서의 급격한 컷오프를 포함하는 물리적 모델의 변형이다.
    3. 경험적 동기 모델: UV 광도 함수 데이터(저 $z$에서는 Madau-Dickinson 형태, 고 $z$에서는 지수적 감소)에 의해 제약되는 모델이다.

• 매개변수: 모델은 6개의 자유 천체물리학적 매개변수인 Ly$\alpha$ 정규화($f_{Ly}$), Lyman 계열 스펙트럼 지수($s$), X-선 정규화($f_X$), X-선 스펙트럼 지수($w$), 이온화 광자의 탈출 분율($f_{esc}$), 그리고 최소 비릴 온도($min(T_{vir})$)에 의해 제어된다. 또한, 표준 우주론적 매개변수($H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$)의 변화를 허용한다.

• 성능: 해석적 접근 방식 덕 인해 ECHO21은 단일 CPU에서 약 $O(1)$초 만에 하나의 전역 21-cm 신호 인스턴스를 생성할 수 있다. 이 코드는 MPI 병렬화를 위해 설계되어, 몇 분 만에 수천 개의 모델을 생성할 수 있다.

주요 기여

1. Ly$\alpha$ 가열의 포함: ECHO21은 IGM의 열적 진화에 Ly$\alpha$ 가열을 명시적으로 포함한 최초의 해석적 코드이다. 저자들은 이 항을 무시하는 것이 특히 X-선 방출이 낮고 Ly$\alpha$ 방출이 강한 시나리오에서 상당한 편향을 초래할 수 있음을 입증하였다.
2. SFRD 모델의 다양성: 코드는 SFRD를 모델링하는 세 가지 다른 접근 방식(물리적, 준경험적, 경험적)을 독특하게 제공하여, 사용자가 별 형성 규정(prescription)의 선택이 결과적인 21-cm 신호에 어떤 영향을 미치는지 평가할 수 있게 한다.
3. 우주론적 유연성: 많은 전문화된 코드들과 달리, ECHO21은 천체물리학적 매개변수와 동일한 용이성으로 우주론적 매개변수를 다루어, 21-cm 신호에 미치는 우주론적 불확실성(예: 허블 텐션)에 대한 연구를 용이하게 한다.
4. 속도 및 확장성: ECHO21은 빠른 추론을 위해 최적화되어 있으며, 50개의 CPU를 사용하여 6차원 천체물리학적 매개변수 공간에서 15,625개의 모델을 약 17분 만에 생성할 수 있다.

결과

• 표준 모델(Fiducial Model): 물리적 동기 SFRD와 Press-Schechter 헤일로 질량 함수를 사용하는 표준 모델은 $z \approx 5.3$까지 완료되는 재이온화 역사를 예측하며, 50% 이온화 적색편이는 $z_{re} = 7.42$이다. 전자 산란 광학 깊이는 $\tau_e \approx 0.0592$로 계산되었으며, 이는 Planck 2018 제약 조건과 일치한다.
• 신호 특성: 표준 전역 21-cm 신호는 $z=87.5$(암흑기)에서 $-41.1$ mK의 첫 번째 흡수 피크와, $z=19.2$(우주 새벽)에서 $-142.4$ mK의 더 깊은 두 번째 흡수 피크를 보인 후, $z=11.2$에서 $22.2$ mK의 방출 피크를 나타낸다.
• Ly$\alpha$ 가열의 영향: Ly$\alpha$ 가열을 제외한 모델과의 비교를 통해, 이 항의 누락이 더 강한(더 깊은) 21-cm 흡수 신호를 유발함을 보여준다. 저자들은 Ly$\alpha$ 가열이 냉각 효과에 대항하고 X-선 매개변수에 따라 비단조적(non-monotonically)으로 신호 강도를 변화시킬 수 있으므로, 정확한 모델링을 위한 필수적인 구성 요소라고 강조한다.
• SFRD 모델의 영향: SFRD 모델의 선택은 신호를 크게 변화시킨다. 물리적 모델에 비해 별 형성의 시작을 늦게 예측하는 경험적 모델은 신호를 약 8.8의 적색편이만큼 지연시키고 약 39.4 mK만큼 얕게 만든다. 마찬가지로, 헤일로 질량 함수(HMF)의 선택 또한 신호 특징의 타이밍과 깊이에 영향을 미친다.
• 매개변수 민감도: 신호는 $f_{Ly}, f_X, w, min(T_{vir})$에 가장 민감하다. 현재의 관측 불확실성 범위 내에서 우주론적 매개변수를 변화시켜도 신호 진폭(예: 피크 흡수 깊이의 약 $\sim 18$ mK 범위)에는 상대적으로 작은 변화만을 준다.

의의 및 주장
저자들은 ECHO21을 다가오는 21-cm 데이터셋(예: REACH, SARAS 또는 미래의 우주 기반 실험)으로부터 "빠르지만 충분히 현실적인 천체물리학적 추론"을 수행하기 위해 설계된 유연하고 확장 가능한 오픈 소스 도구로 정의한다. 그들은 이 코드가 다음과 같은 작업에 특히 적합하다고 주장한다:

• 대규모 매개변수 탐색을 통한 견고한 천체물리학적 및 우주론적 추론 수행.
• 순수하게 경험적인 피팅이 아닌, 물리적으로 동기화된 신호를 사용하여 데이터 분석 파이프라인(특히 REACH 협업을 위한)을 테스트.
• 21-cm 관측과 BBN 예측, CMB 이방성, 허블 텐션과 같은 다른 우주론적 탐사 도구 간의 시너지 구축.

논문은 ECHO21이 (다중 Ly$\alpha$ 산란 등을 무시하는 해석적 근사로 인해) 준수치적 모델보다는 덜 정확할 수 있음을 겸허히 인정하면서도, 그 속도와 모듈성이 초기 추론을 수행하고 향후 버전에서 새로운 물리학을 통합하는 데 이상적인 도구임을 강조한다. 코드는 PyPI와 GitHub를 통해 공개적으로 이용 가능하다.