21cmFAST에 라잎다중산란을 반영한 반분석적 접근

초록

라잎α 광자의 감쇠날개에 의한 다중산란을 반분석적으로 모델링하고, 이를 21cmFAST에 적용해 고전적인 직선 근사와 비교하였다. Monte‑Carlo 시뮬레이션으로 얻은 거리 분포를 단일 파라미터 베타 함수로 근사화했으며, 다중산란이 21 cm 전력 스펙트럼에 미치는 영향은 고령도(스핀 온도 결합 전)에서만 약 50% 정도 차이를 만든다. 라잎α 가열 효과는 거의 변하지 않아, 현실적인 X‑레이 가열 시나리오에서는 무시할 수 있다.

상세 분석

이 논문은 21 cm 신호 시뮬레이션에서 라잎α 광자의 다중산란(MS) 효과를 효율적으로 포함시키는 새로운 방법을 제시한다. 기존의 반수치 시뮬레이터들은 라잎α 교차섹션을 δ‑함수로 근사해, 라잎β와 라잎α 사이에서 방출된 광자가 레드시프트에 의해 정확히 공명 주파수에 도달할 때까지 직선 경로를 따라 이동한다고 가정한다. 그러나 라잎α 교차섹션의 감쇠날개는 광자가 공명에 접근할수록 점차 큰 확률로 산란하게 하며, 이로 인해 산란 간 거리가 비무시적이게 된다. 이러한 현상은 특히 스핀 온도가 아직 기체 온도와 완전히 결합되지 않은 고령도(z ≈ 20–30) 구간에서 라잎α 결합 계수 xα의 공간 변동성을 크게 바꿀 수 있다.

저자들은 먼저 “Speedy Lyα Ray Tracing Algorithm”(SPαRTA)라는 공개 Monte‑Carlo 코드를 개발했다. 이 코드는 격자 기반 라디에이션 전이 시뮬레이션보다 몇 초에서 몇 분 안에 실행되며, 선형 섭동 이론을 이용해 peculiar velocity와 유한 온도 효과를 반영한다. SPαRTA를 이용해 다양한 IGM 온도와 팽창 조건에서 라잎α 광자의 흡수점까지의 거리 분포를 측정했으며, 그 결과가 베타 분포 형태를 띠고 단일 스케일 파라미터 β에 의해 완전히 기술된다는 것을 발견했다. β는 주로 광자의 초기 주파수와 IGM의 온도·밀도에 의해 결정되며, β → ∞ 일 때는 기존의 직선(δ‑함수) 근사와 일치한다.

이 베타 분포를 윈도우 함수 W(k;z,z′)의 형태로 Fourier 공간에 삽입함으로써, 21cmFAST의 기존 구조를 크게 변경하지 않고도 MS 효과를 포함시킬 수 있었다. 구체적으로, 라잎α 플럭스 Jα는 원래의 emissivity 필드와 윈도우 함수를 곱한 뒤 역변환하는 방식으로 계산된다. 이 과정에서 거리 제한 RSL(z,z′)를 사용해 라잎α 광자가 도달할 수 있는 최대 comoving 거리까지를 고려했으며, 감쇠날개에 의한 비직선 경로는 베타 분포가 제공하는 확률적 가중치로 대체되었다.

시뮬레이션 결과는 다음과 같다. (1) 고령도에서 xα ≈ 1 수준일 때, MS 효과는 라잎α 플럭스의 평균값을 약 20% 증가시키고, 전력 스펙트럼 P21(k)에서는 k ≈ 0.1 Mpc⁻¹ 부근에서 최대 50% 차이를 만든다. (2) 스핀 온도가 기체 온도와 완전히 결합된 이후(z ≲ 15)에는 xα≫1이 되므로, 플럭스 변동이 신호에 미치는 영향이 급격히 감소한다. (3) 라잎α 가열은 광자 산란 횟수와 에너지 전달 효율이 거의 변하지 않기 때문에, MS를 포함하더라도 전체 가열률은 5% 이하로 변동한다. 특히 X‑레이 가열이 지배적인 시나리오에서는 라잎α 가열 기여가 무시 가능 수준이다.

이러한 결과는 두 가지 중요한 함의를 가진다. 첫째, 파라미터 탐색을 위한 대규모 MCMC 혹은 신경망 기반 추론에 있어, 라잎α MS를 완전한 라디에이션 전이 시뮬레이션 없이도 충분히 정확히 모델링할 수 있음을 보여준다. 둘째, 현재 및 차기 21 cm 관측(예: HERA, SKA)에서 고령도 구간의 전력 스펙트럼을 해석할 때만 MS를 고려하면 되고, 낮은 레드시프트에서는 기존 직선 근사가 충분히 정확하다는 실용적인 가이드를 제공한다.