지역 은하 내 우주선 이온화율 3차원 지도: H₃⁺ 흡수와 3D‑PDR 모델링의 결합

초록

H₃⁺ 흡수선을 이용해 27개 확산 분자 구름 시선 중 12개에서 우주선에 의한 H₂ 이온화율 ζ(H₂)를 측정하고, 3D‑PDR 모델과 Gaia‑Edenhofer 차등 소광 지도에 기반한 밀도 프로파일을 결합해 1 kpc 이내의 3차원 ζ(H₂) 지도를 구축하였다. 평균 ζ(H₂)=5.3×10⁻¹⁷ s⁻¹(σ=2.5×10⁻¹⁷ s⁻¹)이며, 수십 파섹 규모에서는 거의 균일하지만 100 pc 정도에서는 5배 차이가 나타난다. 이는 우주선 가속원 근접성에 따라 이온화율이 변한다는 시나리오를 지지한다.

상세 분석

본 연구는 확산 분자 구름에서 H₃⁺가 우주선에 의한 H₂ 이온화의 가장 직접적인 지표임을 재확인하고, 최신 적외선 고해상도 분광기(iSHELL, CRIRES)를 활용해 3.7–3.9 µm 파장에서 네 개의 전이(R(1,1)ᵤ, R(1,0), R(1,1)ₗ, Q(1,1))를 관측하였다. 데이터 감소 과정에서는 플랫필드, 대기선 보정, LSR 변환, 그리고 다중 관측 결합을 체계적으로 수행했으며, CH 4300 Å 흡수를 템플릿으로 사용해 H₃⁺ 라인의 중심과 폭을 고정하고 깊이만 자유롭게 하여 약한 신호에서도 신뢰할 수 있는 등가폭을 도출하였다.

화학적 해석에서는 H₂⁺ → H₃⁺ 생성과 전자 재결합 파괴 과정을 가정한 정적 평형 모델을 이용해 ζ(H₂)=kₑ xₑ n_H N(H₃⁺)/N(H₂) 식을 도출하였다. 여기서 kₑ는 전자 재결합 계수, xₑ는 전자분율, n_H는 수소 핵 밀도이며, N은 관측된 컬럼밀도이다. 이 식은 1D 화학 모델과 3D‑PDR 시뮬레이션 결과와 2배 이내의 일치를 보이며, 복잡한 선형 밀도 구조를 가진 경우 ζ(H₂)/n_H 형태로 추정하는 것이 합리적임을 제시한다.

3D‑PDR 모델링은 Edenhofer 등(2024)의 차등 소광 지도에서 파생된 3D 가스 밀도 분포를 입력으로 사용하였다. 소광‑밀도 변환 계수 n_H=1710 (dE_GRZ/ds) cm⁻³를 채택하고, 12개의 후방 샘플 평균을 이용해 각 시선에 대한 1D 밀도 프로파일을 생성했다. 이러한 밀도 정보와 관측된 H₃⁺, H₂ 컬럼을 결합해 베이시안 최적화를 수행, 각 구름에 대한 ζ(H₂) 값을 도출하였다. 결과적으로 전체 표본의 평균 ζ(H₂)는 5.3×10⁻¹⁷ s⁻¹이며, 표준편차는 2.5×10⁻¹⁷ s⁻¹이다.

공간적 분석에서는 ζ(H₂)가 수십 파섹 규모에서는 거의 일정하지만, 약 100 pc 거리에서는 최대 5배 차이가 나타나는 지역을 확인했다. 이는 우주선 가속원(예: 초신성 잔해, 펄서 풍선)에서의 거리 의존적 확산 및 에너지 손실 메커니즘이 지역적인 이온화율 변동을 주도한다는 기존 이론을 실증적으로 뒷받침한다. 또한, 이전 H₃⁺ 기반 조사에서 제시된 ζ(H₂)≈3.5×10⁻¹⁶ s⁻¹와 비교해 약 1 dex 낮은 값을 얻음으로써, 최신 밀도 재평가와 3D 모델링이 이온화율 추정에 미치는 영향을 강조한다.

이 연구는 (1) 고해상도 적외선 흡수선 관측, (2) CH 기반 라인 프로파일 템플릿, (3) 3D‑PDR 시뮬레이션, (4) Gaia‑Edenhofer 차등 소광 지도라는 네 가지 최신 기법을 통합함으로써, 지역 은하 내 우주선 이온화율의 3차원 분포를 최초로 정량화한 점에서 의의가 크다. 향후 더 많은 시선과 고밀도 구름을 포함한 확장 연구가 우주선 전파 모델과 ISM 화학 모델을 정밀하게 연결하는 데 핵심 자료가 될 것이다.