우주 먼지 위에서 원자 첨가로 형성되는 HNSO: 산소·질소 확산이 이끄는 새로운 황 화학 경로

초록

본 연구는 양자 화학 계산과 천체 화학 모델링을 결합해, 암석질 물(ASW) 표면에서 O·N 원자 확산에 의해 NS와 SO 라디칼과 반응하여 중간체 NSO를 형성하고, 이후 수소 첨가로 HNSO(시스형)가 효율적으로 생성된다는 것을 밝힌다. 모델은 고체 HNSO가 OCS와 비슷한 어는 상태 풍부도를 가질 수 있음을 예측한다.

상세 분석

이 논문은 HNSO(시스‑HNSO)라는 네 원소(H, N, S, O)로 구성된 소분자의 우주적 기원을 규명하기 위해 세 단계의 접근을 취한다. 첫째, ωB97M‑D4‑gCP/def2‑SVP 수준에서 구조 최적화와 진동수 계산을 수행하고, ωB97M‑D4/def2‑TZVPPD로 단일점 에너지를 보정한 DFT 방식을 사용하였다. 이는 다중구성 특성을 보이는 반응 경로에 대해 실용적인 정확도를 제공한다. 두 번째로, 33개의 물 분자로 구성된 비정질 고체 물(ASW) 클러스터 모델에 SO와 NS 라디칼을 세 가지 결합 부위(Pocket, dH, Pentamer)에 흡착시켜 각각의 결합 에너지(BE)를 구했으며, 평균 BE는 SO≈4.4 kcal mol⁻¹(≈2220 K), NS≈5.4 kcal mol⁻¹(≈2720 K)로 나타났다. 이는 원자 O·N이 ASW 표면에서 비교적 낮은 장벽으로 확산할 수 있음을 시사한다. 세 번째로, O·N 원자 첨가 반응을 12가지 경우(3 부위 × 2 라디칼 × 2 첨가 위치)로 조사했으며, 대부분이 무장벽(barrierless) 혹은 매우 낮은 활성화 에너지(≤0.3 kcal mol⁻¹)로 진행됨을 확인했다. 특히 SO+N→NSO와 NS+O→NSO 반응은 각각 –70.8 kcal mol⁻¹, –80.4 kcal mol⁻¹의 반응 엔탈피를 가지며 거의 무장벽이다. 반면, NO+S 경로는 높은 재배열 장벽으로 비현실적이라고 판단한다.

이후 NSO의 수소화 경로를 탐색했으며, H + NSO→HNSO는 시스형이 가장 안정적이고, 트랜스형은 약 3 kcal mol⁻¹ 상위에 머무른다. 트랜스→시스 이성질 전이는 CCSD(T)-F12/VTZ‑F12 수준에서 자유 에너지 장벽이 30 kcal mol⁻¹ 이상으로, 터널링 전이 속도 상수 κ는 천문학적 시간(>10⁶ yr)보다 훨씬 느리다. 따라서 우주 환경에서는 주로 시스‑HNSO가 관측될 것으로 기대된다.

천체 화학 모델링에서는 다중 결합(multibinding) 스키마를 도입해, 표면에서 O·N 원자의 확산 확률을 결합 에너지 분포에 따라 가중 평균하였다. 초기(≤10⁵ yr) 단계에서 O·N 확산이 활성화되면 NSO와 HNSO의 고체 상대 풍부도가 급격히 상승하며, 최종적으로 HNSO/OCS 비율이 0.5–1.0 수준에 도달한다. 가스상에서는 HNSO가 OCS보다 1–2 dex 낮게 유지되며, 이는 관측된 가스‑상 HNSO 비감소와 일치한다. 모델은 또한 NS와 SO 라디칼의 초기 풍부도가 HNSO 생산에 결정적임을 보여, G+0.693와 같은 충격‑유도 해방 환경에서 이러한 라디칼이 풍부히 존재함을 강조한다.

결과적으로, 이 연구는 (1) O·N 원자 확산이 비수소화 경로를 열어 NSO를 형성하고, (2) NSO의 수소화가 시스‑HNSO를 생성하며, (3) 다중 결합 스키마가 초기 확산 효율을 정확히 반영함으로써 관측 풍부도와 모델을 일치시킨다는 세 가지 핵심 통찰을 제공한다. 이는 황‑산소‑질소 복합 화학이 복잡 유기 분자와 별개로도 풍부히 진행될 수 있음을 시사하고, 향후 HNSO와 유사한 NSO‑계열 분자의 탐색을 촉진한다.