SN 1987A의 3차원 분자 지도와 폭발 시뮬레이션 대조 연구

초록

ALMA 고해상도 관측으로 SN 1987A의 CO와 SiO 분자를 3차원으로 복원하고, Prometheus‑HotB 코드로 만든 10가지 전구체 모델과 비교하였다. 이산화탄소‑산소 풍부한 핵심 영역은 이진 합병 청색거성 모델과 가장 잘 맞으며, SiO는 외부 혼합이 부족해 관측과 차이를 보인다. 또한, 가장 무거운 C/O 물질의 질량중심이 중성자별 킥 방향과 일치한다는 이론적 예측을 제시한다.

상세 분석

본 논문은 SN 1987A의 분자 방출을 ALMA로 0.03–0.06″(≈10–20 AU) 해상도로 두 시점(≈10 000 일, 11 800 일)에서 관측한 뒤, 복셀 단위의 3차원(좌표 x‑y‑z) 분포를 재구성한 것이 핵심이다. 복원 과정에서 방출 강도와 속도 채널을 결합해 자유 팽창(homologous expansion) 가정하에 각 픽셀을 ‘속도 × 시간’으로 거리로 변환하였다. 결과적으로 CO는 비교적 평탄하고 원반형에 가까운 구조를, SiO는 고속 외피에 집중된 비대칭 클럼프(clump) 형태를 보였다.

이러한 관측 결과를 10가지 전구체 모델(레드·블루 초거성, 이진 합병 청색거성 등)과 비교하기 위해, Prometheus‑HotB 코드로 1 초에서 수천 초까지의 초기 폭발 과정을 시뮬레이션한 뒤, 장기(수천 년) 진화를 추적해 물질 혼합과 속도 분포를 얻었다. 모델별로 CO와 SiO의 질량‑속도 분포를 정량화하고, 3차원 형태를 시각적으로 비교하였다.

가장 좋은 일치는 ‘이진 합병 청색거성’ 시나리오였다. 이 모델은 핵심에서 C와 O가 풍부한 물질이 강하게 집중되고, 그 질량중심이 폭발 비대칭성에 따라 이동한다는 점에서 CO 관측과 일치한다. 특히, 두 개 이상의 모델이 CO의 원반형 구조와 고속 외피를 재현했으며, 관측된 ‘중심 공백’(velocity≈0 km s⁻¹에서의 낮은 방출)도 자연스럽게 나타난다.

반면 SiO는 모델이 충분히 외부까지 Si를 섞어 넣지 못해, 관측에서 보이는 고속(>3000 km s⁻¹) 클럼프가 부족하다. 이는 전구체 단계에서 Si가 이미 외부에 충분히 분포되지 않았거나, 시뮬레이션의 혼합 효율(예: Rayleigh‑Taylor 불안정성)과 네트워크가 제한적임을 시사한다. 저자들은 이 문제를 해결하기 위해 초기 전구체의 화학 구성을 재조정하거나, 3D 네트워크 해상도를 높여야 한다고 제언한다.

또한, 시뮬레이션은 가장 무거운 C/O‑풍부 물질의 질량중심이 중성자별(NS) 킥 방향과 강하게 상관관계가 있음을 보여준다. 관측된 CO 방출이 가장 강한 영역과 질량중심이 일치한다면, NS는 관측자 쪽으로 약 45° 북쪽 방향으로 움직이고 있을 것으로 예측한다. 이는 기존 X‑ray·광학에서 추정된 NS 위치와도 일치한다(관측된 비대칭성, 라디오 파동의 비대칭 팽창 등).

결론적으로, 고해상도 ALMA 분자 지도는 폭발 초기 물리(핵심 혼합, 비대칭성, NS 킥)와 직접 연결되는 귀중한 데이터이며, 현재의 3D 수치 모델은 CO 재현에 성공했지만 SiO와 같은 고온·고속 물질의 혼합을 아직 충분히 설명하지 못한다는 점을 강조한다. 향후 더 정교한 전구체 모델링과 장기 3D 수치 시뮬레이션이 필요하다.