바나드 5에서 보는 탈수소 비율과 CO 고정의 새로운 지도

초록

IRAM 30 m 텔레스코프를 이용해 Barnard 5의 별이 없는 핵과 원시성 핵을 동시에 관측했다. N₂H⁺/N₂D⁺와 HCO⁺/DCO⁺의 탈수소 비율과 C¹⁸O(2–1) 데이터를 결합해 CO 고정도(f_d)를 지도화했으며, 별이 없는 핵에서 탈수소 비율이 0.43(N₂H⁺)·0.09(HCO⁺)까지 상승하고 CO 고정도는 5.0에 이른다. 반면 원시성 핵에서는 각각 0.15, 0.05, f_d ≈ 4.1이다. 결과는 프로토스타의 열적 피드백이 고밀도 핵 내부의 탈수소를 억제하지만, 저밀도 외피에서는 미세하게만 영향을 미친다는 점을 확인한다.

상세 분석

이 연구는 Barnard 5(B5)라는 퍼시스스 분자 구름 내에서 별이 없는 핵과 Class 0/I 원시성 핵(IRAS 03445+3242)을 동시에 조사함으로써 탈수소 화학과 CO 고정 사이의 상관관계를 공간적으로 해석한 최초 사례이다. 관측은 IRAM 30 m에서 N₂H⁺(1–0), N₂D⁺(1–0), H¹³CO⁺(1–0), DCO⁺(2–1) 및 C¹⁸O(2–1) 전이를 이용했으며, 동일한 공간 해상도(33.6″)와 픽셀 크기(12″)로 재처리하였다. 하이퍼파인 구조를 갖는 전이들의 광학 깊이와 여기 온도(T_ex)를 pyspeckit의 HFS 피팅으로 추정하고, 마코프 체인 몬테카를로(MCMC) 기법을 통해 파라미터 공간을 탐색해 불확실성을 최소화하였다.

컬럼 밀도는 광학 깊이와 T_ex에 민감하게 반응함을 확인했으며, 특히 N₂H⁺와 N₂D⁺는 T_ex가 1 K 변할 때 20–27 % 정도 변동한다. 따라서 전체 지도에 동일한 T_ex 값을 적용했으며, 이는 별이 없는 핵과 원시성 핵 사이의 온도 차이가 1 K 이하라는 관측 결과와 일치한다.

탈수소 비율 R_D는 N₂H⁺와 HCO⁺ 각각에 대해 정의했으며, 별이 없는 핵에서는 R_D(N₂H⁺)=0.43±0.10, R_D(HCO⁺)=0.09±0.02로, 원시성 핵에서는 0.15±0.03, 0.05±0.01에 불과했다. CO 고정도 f_d는 C¹⁸O와 ¹³CO(1–0) 전이로부터 전체 CO 풍부도를 추정해 계산했으며, 별이 없는 핵에서 5.0±0.1, 원시성 핵에서 4.1±0.1을 보였다.

이러한 결과는 두 가지 핵심 메커니즘을 시사한다. 첫째, CO가 기체 상태에 존재하면 H₃⁺와의 반응을 통해 H₂D⁺ 생성이 억제되어 탈수소 비율이 낮아진다. 별이 없는 핵에서는 CO가 고정(얼음)되어 있어 H₃⁺가 H₂D⁺로 전환될 여지가 커지고, 따라서 높은 R_D가 관측된다. 둘째, 원시성 핵 주변에서는 열적 피드백(≈15 K)으로 인해 CO가 다시 기체 상태로 방출되고, 동시에 온도 상승으로 H₂D⁺↔H₃⁺ 평형이 역전되어 탈수소 비율이 감소한다.

흥미롭게도 HCO⁺를 추적하는 저밀도 외피에서는 R_D(HCO⁺) 감소가 상대적으로 작았다. 이는 HCO⁺가 주로 n(H₂)≈10⁴–10⁵ cm⁻³ 영역을 추적하고, 이 영역은 원시성 핵의 열적 영향이 제한적이기 때문이다. 따라서 탈수소 화학은 밀도와 온도에 따라 층화된 구조를 보이며, 고밀도 핵 내부와 저밀도 외피가 서로 다른 화학적 진화를 겪는다.

이 논문은 기존 단일점 관측에 비해 공간 해상도를 크게 향상시켜, 탈수소 비율과 CO 고정이 동시에 변하는 모습을 직접 시각화했다. 또한, Barnard 5에서 보고된 “높은 탈수소 비율·낮은 CO 고정”이라는 이례적 현상이 실제로는 별이 없는 핵과 원시성 핵이 혼재된 복합 구조 때문임을 입증한다. 결과는 화학 모델링에 온도·밀도·CO 고정의 상호작용을 보다 정밀하게 구현해야 함을 강조한다.