우주 환경에서 페닐리움 반응성 재조명

초록

본 연구는 최근 Kocheril 등(2025)이 제시한 페닐리움(C₆H₅⁺)의 비반응성을 반박한다. VUV 광이온화를 이용해 순수한 페닐리움을 생성하고, 아세틸렌(C₂H₂)과의 반응을 조사한 결과, 장벽 없는 입체적 부가 반응이 저온 우주 환경에서도 효율적으로 진행됨을 확인하였다. DFT 계산(M06‑2X/AVTZ) 역시 반응 경로에 에너지 장벽이 없음을 보여주며, 실험적 교차섹션과 반응 속도 상수는 이론값과 일치한다. 압력 차이에 따른 생성물 분포 차이는 설명 가능하지만, 반응 자체가 억제된 것은 아니다. 따라서 페닐리움은 첫 번째 방향족 고리 형성에 중요한 전구체로 남아 있으며, 기존 astrochemical 모델에 반영할 필요가 있다.

상세 분석

이 논문은 두 가지 핵심 질문에 답한다. 첫째, 페닐리움이 아세틸렌과 반응할 수 있는가? 둘째, Kocheril 등에서 제시된 “비반응성”이 실험적·이론적 조건 차이에 기인한 오류인지 여부이다. 실험적으로 저에너지 VUV(≈11 eV) 광자를 이용해 니트로벤젠을 광이온화함으로써 거의 순수한 C₆H₅⁺ 이온을 생성한다. 이는 기존 연구에서 사용된 고에너지 전자 충격이나 충돌 이온화와 달리, 비반응성 이성질체(예: l‑C₆H₅⁺)가 동시에 생성될 가능성을 최소화한다. 생성된 이온을 8.5 × 10⁻⁵ mbar의 아세틸렌 가스와 단일 충돌 조건에서 반응시켰으며, 질량 분석을 통해 m/z 103 (C₈H₇⁺) 부가체와 m/z 102 (C₈H₆⁺) H 손실 생성물을 정량화했다. 교차섹션은 내부 에너지 0–3 eV 범위에서 급격히 감소하지만, 13 eV 이하에서는 부가체 형성이 주된 경로임을 확인했다. 이는 이전 실험(예: Knight et al., 1987)에서 보고된 4 × 10⁻¹⁰ cm³ s⁻¹와 거의 일치한다.

이론적으로는 M06‑2X/AVTZ 수준의 DFT 계산을 수행해 반응 경로의 잠재 에너지 면을 탐색하였다. Figure 1a에서 보여지듯, C₂H₂가 C₆H₅⁺의 비수소 탄소에 접근할 때 장벽이 전혀 존재하지 않는다. 최적화된 전이 상태와 중간체(1‑3)는 모두 0 kJ mol⁻¹ 이하의 장벽을 보이며, 제로점 에너지 보정 후에도 반응은 발열적이다. 특히, 5‑membered 고리 형성(iii)만이 약간의 발열을 보이지만, 전체 반응은 에너지적으로 유리하다. 이러한 “장벽 없는” 특성은 저온(10–100 K) 우주 환경에서도 반응이 일어날 수 있음을 의미한다.

압력 효과에 대한 논의도 중요하다. Kocheril 등은 10⁻⁸ mbar 수준의 초저압에서 실험했으며, 이는 부가체가 충돌에 의해 안정화될 기회가 거의 없음을 뜻한다. 반면 본 연구는 10⁻⁵ mbar 수준으로, 충돌에 의한 에너지 전달이 충분히 일어나 부가체가 안정화된다. 따라서 두 실험의 차이는 반응성 자체가 아니라 생성물의 분포와 검출 가능성에 기인한다. 또한, C₄H₃⁺+C₂H₂ 경로가 실제로는 비환원성 l‑C₆H₅⁺를 주로 생성한다는 선행 AIMD 연구(Peverati et al., 2016)를 인용해, Kocheril 연구에서 관찰된 “비반응성”은 잘못된 전구체 할당에서 비롯된 것으로 결론짓는다.

천체화학 모델링 측면에서는 Nautilus 3‑phase 코드를 이용해 업데이트된 네트워크를 시뮬레이션했다. C₆H₅⁺+H₂ 라디에티브 결합 반응은 온도와 압력에 거의 민감하지 않아, ISM(Interstellar Medium) 조건에서도 10⁻⁹ cm³ s⁻¹ 수준의 속도 상수를 유지한다. 모델 결과는 C₆H₅⁺가 실제로 C₈H₇⁺(=C₆H₅CCH₂⁺)와 C₈H₆⁺를 통해 C₆H₇⁺(벤젠 전구체)로 전환되는 경로가 주요함을 보여준다. 그러나 C₄H₃⁺+C₂H₂ 경로를 제외해도 벤젠 및 기타 방향족 종의 총량은 1 dex 정도만 감소한다. 이는 페닐리움 자체가 “핵심 전구체”라기보다 “다양한 경로 중 하나”임을 시사한다.

결론적으로, 이 논문은 (1) 페닐리움은 아세틸렌과 장벽 없는 반응을 수행한다, (2) Kocheril 등에서 보고된 비반응성은 실험 조건과 전구체 할당 오류에 기인한다, (3) 페닐리움은 여전히 첫 번째 방향족 고리 형성에 기여하지만, 모델에 반영할 때는 다른 효율적인 경로와 함께 고려해야 함을 강조한다.