3차원 시뮬레이션으로 밝힌 Ia형 초신성 잔해의 X‑레이 다양성

초록

이 연구는 순수 결핍, 지연 폭발, 이중 폭발 모델 각각 두 가지씩 총 6가지 Ia형 초신성 폭발 시나리오를 3차원 유체역학과 비평형 이온화(NEI)를 결합해 1000 yr까지 진화시킨 뒤, XRISM 수준의 1 eV 해상도로 X‑레이 스펙트럼을 합성한다. Fe‑Kα 중심 에너지와 라인 세기, 비대칭적인 적·청색 이동 프로파일 등 모델마다 뚜렷한 차이를 보이며, 이는 관측된 Ia형 초신성 잔해의 다양성을 3‑D 모델링으로 재현하고 폭발 메커니즘 및 전구체 시스템을 정성적으로 구분할 수 있음을 시사한다.

상세 분석

본 논문은 Ia형 초신성(SN) 폭발 메커니즘의 다양성을 SNR 단계까지 일관되게 추적하는 최초의 3‑D 시뮬레이션을 제시한다. 핵심 기술은 VH‑1 기반 PPMLR 스키마에 비평형 이온화(NEI) 모듈을 결합하고, SNR이 팽창함에 따라 격자 자체를 동적으로 확대하는 ‘코모빙’ 그리드 방식을 도입한 점이다. 이를 통해 초기 SN 단계(수년)부터 1000 yr에 이르는 SNR 규모(수십 파섹)까지 해상도 손실 없이 계산할 수 있었다. 트레이서 입자를 10만 개(70 k ejecta, 30 k ISM) 배치해 각 입자별 전자·이온 온도, 밀도, 원소 조성, 이온화 상태를 실시간으로 업데이트함으로써, 후처리가 아닌 수치 흐름과 동기화된 NEI 계산을 수행했다. 이온화·재결합율은 ATOMDB/APED 데이터를 사용했으며, 충격 전후 전자와 양성자 온도 비율을 m_e/m_p 로 설정해 초기 전자 가열 효과를 최소화했지만, 금속 풍부한 ejecta의 경우 전자는 충돌 이온화에 의해 빠르게 냉각되므로 결과에 큰 영향을 주지 않는다.

폭발 모델은 N100ddt·N5ddt(지연 폭발), N100def·N5def(순수 결핍), Pakmor et al. (2022)에서 제시한 OneExp·TwoExp(이중 폭발) 등 여섯 가지이며, 각각 질량·에너지·56Ni·IGE·IME 비율이 표에 정리돼 있다. 초기 조건은 3‑D SN 모델을 구면 평균해 1‑D 라디얼 프로파일로 시각화했으며, 자유 팽창 단계(≤3 yr) 이후 균일 밀도(ρ≈2×10⁻²⁴ g cm⁻³) ISM에 진입하도록 설정했다.

시뮬레이션 결과는 크게 두 축으로 분석된다. 첫째, 플라즈마 물리량(전자 온도, 밀도, 이온화 시간)과 원소별 부피 분포가 시간에 따라 어떻게 진화하는가이다. 지연 폭발 모델은 높은 56Ni 함량과 강한 내부 압력으로 빠른 역충격을 일으켜, Fe‑Kα 라인의 중심 에너지가 6400–6500 eV 사이에 머무는 반면, 결핍 모델은 낮은 에너지(≈6420 eV)와 넓은 라인 폭을 보였다. 이중 폭발 모델은 외부 He‑detonation에 의해 얇은 고온 껍질이 형성돼, Fe‑Kα 라인이 상대적으로 높은 에너지(≈6540 eV)와 강한 청색 이동을 나타냈다. 둘째, 합성된 X‑레이 스펙트럼은 SOXS 툴을 이용해 XRISM/Resolve의 응답을 적용, 1 eV 해상도로 재현했다. 라인 프로파일은 비대칭적이며, 특정 방향에서는 2000 km s⁻¹ 수준의 적·청색 이동이 관측된다. 이는 3‑D ejecta 구조가 비구형이며, 불균일한 속도장과 혼합이 일어나기 때문이다.

이러한 결과는 관측된 Ia형 SNR(예: Tycho, SN 1006, 0509‑67.5)의 Fe‑Kα 중심 에너지 분포와 라인 세기 변동성을 자연스럽게 설명한다. 또한, XRISM와 같은 고해상도 미션이 제공하는 도플러 시프트 측정을 통해, 특정 라인 프로파일의 비대칭성을 직접 비교함으로써 전구체 시나리오(단일·이중성)와 폭발 메커니즘(지연·결핍·이중) 사이의 차이를 정량화할 수 있음을 시사한다. 다만, 현재 모델은 균일 ISM만을 고려했으며, CSM 상호작용, 입자 가속에 의한 비열 플라즈마, 자기장 효과 등은 제외했다. 향후 이러한 요소를 포함한 파라미터 스터디와 관측 데이터와의 정량적 매칭이 필요하다.