폭발 직전 별의 맥동이 구조를 바꾸다

초록

이 논문은 질량 10.5–15 M⊙의 적색 초거성(RSG)이 핵연료 소진 말기에 κγ‑기전으로 유도되는 방사형 맥동을 겪으며, 맥동 단계에 따라 별의 반지름·광도·유효온도가 크게 변동한다는 것을 보여준다. 이러한 변동은 사전‑초신성 관측값을 단일 L–T 측정만으로는 해석할 수 없게 만들고, 폭발 시 초기 반지름에 따라 초기 광도 급증(early‑excess)과 광속도 분포가 달라진다. 저자들은 SN 2023ixf와 SN 2024ggi의 관측과 비교해, 특히 SN 2023ixf는 맥동 모델과 일치하지만 SN 2024ggi는 제시된 주기가 스피처 궤도와 일치해 실제 맥동 증거가 부족함을 지적한다. 결과적으로 15 M⊙ 이상 고질량 RSG에서는 수소층 구조와 초신성 특성이 크게 달라지며, 정적(수력학적) 모델만으로는 충분하지 않다는 결론을 내린다.

상세 분석

본 연구는 MESA와 SNEC를 이용해 비회전, 단일성 적색 초거성의 최후 단계(핵탄소 연소 말기)부터 핵붕괴 직전까지의 수리적 수치 시뮬레이션을 수행하였다. 10.5 M⊙, 12.5 M⊙, 15 M⊙ 세 가지 초기 질량에 대해 인공 점성(artificial viscosity)을 포함한 암시적 수압학 솔버를 적용해 envelope가 κγ‑기전(주로 Fe‑group 원소의 불투명도와 방사압 변화에 의존)으로 자극받아 방사형 맥동을 시작한다는 점을 확인했다. 맥동 주기는 동역학적 시간척도 τ_dyn ≈ √(R³/GM_env)와 거의 비례하며, 질량이 클수록 R_max가 크게 늘어나 τ_dyn가 길어져 주기가 수백에서 천일 수준으로 확대된다. 이는 관측된 RSG의 주기‑광도 관계와도 일치한다(부록 B 참고).

맥동에 따른 별의 구조 변화는 HRD에서 루프 형태로 나타나며, 반지름 변동폭은 15 M⊙ 모델에서 30 % 이상, 광도 변동은 0.5 dex에 달한다. 특히 고질량 모델은 비선형적인 서브구조를 보이며, 이는 envelope 내부의 다중 층 비동조(pulsation mode coupling)와 열 재구성에 기인한다. 이러한 비정상적인 구조는 사전‑초신성 관측에서 단일 L–T 값만으로 핵질량이나 초기 질량을 추정하는 것을 불가능하게 만든다.

폭발 단계에서는 동일한 핵심 매개변수(중성자별 질량, 폭발 에너지, ⁵⁶Ni 질량)를 적용했음에도, 맥동 위상 φ_exp에 따라 외피의 밀도·반지름 프로파일이 크게 달라진다. SNEC 시뮬레이션 결과, 큰 반지름(맥동 최대 확장 시)에서는 충격파가 저밀도 외피를 통과하면서 초기 광도 급증(early‑excess)이 나타나며, 이는 전통적으로 CSM 상호작용으로 해석되는 현상과 유사하다. 또한 초기 광속도(v_ph)도 반지름이 클수록 낮아져, 관측된 Fe II λ 6159 흡수선 최소값과 비교했을 때 고질량 모델은 평균보다 높은 속도를 보이지만, 맥동 위상에 따라 1σ 범위 내에 들어간다.

SN 2023ixf에 대한 적용에서는, 사전 관측된 1100 일 주기의 변동이 모델의 15 M⊙ 맥동 주기와 일치하고, 폭발 직전의 반지름·광도 조합이 초신성 초기 광도와 광속도 곡선 모두를 재현한다. 반면 SN 2024ggi의 경우, 보고된 주기가 스피처 궤도(≈100 일)와 일치해 실제 별의 맥동 증거가 부족함을 저자들은 통계적으로 입증한다.

결론적으로, 15 M⊙ 이상 고질량 RSG에서는 맥동이 envelope 구조와 폭발 후 광도·광속도 특성을 결정적인 수준으로 좌우한다. 기존의 정적(수력학적) 전구 모델은 이러한 동적 변화를 반영하지 못하므로, “missing RSG problem”(관측된 고질량 RSG가 초신성 전구에 부족함) 해소와 초신성 유형 구분에 있어 수리적 수압학 모델이 필수적이다.