Signatures of $^{56}$Ni Mixing and Neutron-rich Ejecta in Supernovae

기술 요약: 초신성에서의 $^{56}\text{Ni}$ 혼합 및 중성자 과잉 이젝타의 징후

문제 제기
초신성(SN) 광도 곡선은 일반적으로 프로제니터(progenitor)의 특성인 이젝타 질량($M_{\text{ej}}$), 니켈 질량($M_{\text{Ni}}$), 폭발 에너지($E_k$) 등을 추론하기 위해 단일 구역 방사성 붕괴 모델(예: Arnett 1980, 1982)을 사용하여 해석된다. 이러한 모델들은 단일 유효 확산 시간 척도를 가정하며 방사성 물질의 공간적 분포를 무시한다. 그러나 유체역학 시뮬레이션에 따르면, $^{56}\text{Ni}$은 이젝타 깊은 곳에서 합성된 후 Rayleigh–Taylor 불안정성, 대류, 난류를 통해 외부로 혼합된다. 이러한 공간적 분포는 단일 구역 모델이 생략하는 결정적인 물리적 변수이다. 결과적으로, 고도로 혼합된 이젝타를 가진 초신성의 광도 곡선은 낮은 이젝타 질량으로 오해받거나, 빠른 상승 및 높은 광도를 설명하기 위해 추가적인 중심 엔진 동력(예: 마그네타 또는 폴백 강착)이 필요한 것으로 잘못 해석될 수 있다. 또한, 콜랩사(collapsar)에서 r-과정(r-process) 농축의 관측 가능성은 니켈에 의해 구동되는 배경과 중성자 과잉 이젝타 사이의 배치 및 기하학적 구조에 민el하게 의존하며, 이는 구형 모델에서 흔히 지나치게 단순화되는 요소이다.

방법론
저자들은 Redback 프레임워크의 JAX 기반 플러그인으로 구현된 준해석적 다중 껍질 모델인 snmix를 소개한다. 단일 구역 모델과 달리, snmix는 호몰로지 팽창(homologously expanding)하는 이젝타를 기하학적으로 간격이 배치된 속도를 가진 $N$개의 껍질로 이산화한다.

• 방사성 구조: 이 모델은 임의의 $^{5 de 6}\text{Ni}$ 방사형 분포(top-hat 또는 core-plus-tail 프로파일 내의 혼합 분율 $f_{\text{mix}}$로 매개변수화됨)와 r-과정 물질을 허용한다.
• 수송 물리: 각 껍질은 독립적으로 진화하며 열 에너지, 단열 손실, 복사 확산을 추적한다. 결정적으로, 모델은 전역 평균이 아닌 각 특정 껍질의 외부에 있는 누적 칼럼 밀도(column density)를 기반으로 확산 및 감마선 퇴적을 계산한다.
• 불투명도 및 SED: 모델은 회색 불투명도(gray opacities)와 온도 의존적 불투명도 법칙을 포함한다. r-과정 계산을 위해, 란타넘족 의존적 불투명도와 핵 네트워크 계산을 통해 보정된 해석적 가열률을 포함한다.
• 기하학적 구조: 주로 구형이지만, 콜랩사를 위한 2-성분 "디스크-윈드(disk-wind)" 기하학으로 확장되어, 시선 방향에 따라 가중치가 부여된 극 방향의 SN 성분과 적도 방향의 중성자 과잉 윈드를 결합한다.
• 추론: 저자들은 이 프레임워크를 사용하여 합성 볼로메트릭 광도 곡선을 생성하고, 단일 구역 피팅을 다중 껍질의 실제 값과 비교하여 매개층 편향(parameter biases)을 분석하기 위해 베이지안 추론(NUTS 샘플링 사용)을 수행한다.

주요 기여 및 결과

1. $^{56}\text{Ni}$ 혼합이 광도 곡선에 미치는 영향:

   • $^{56}\text{Ni}$의 외부 혼합은 덮여 있는 확산 칼럼 및 감마선 포획 칼럼을 감소시킨다. 이는 고정된 $M_{\text{Ni}}$, $M_{\text{ej}}$, $E_k$에 대해 더 빠르고 밝은 상승과 더 가파른 방사성 꼬리(tail)를 초래한다.
   • 빠른 상승과 밝은 상승은 낮은 이젝타 질량이나 엔진의 존재에 대한 독점적인 증거가 아니다. 고도로 혼합된 니켈 모델은 추가적인 동력원 없이도 이러한 특징들을 재현할 수 있다.

2. 단일 구역 추론에서의 편향:

   • 혼합된 광도 곡선을 단일 구역 모델로 피팅하면 시각적으로는 좋은 피팅을 얻지만, 물리적 매개변수에는 편향이 발생한다.
   • 고정된 불투명도에서, 모델은 누락된 방사형 구조를 흡수하기 위해 현저히 낮은 $M_{\text{ej}}$와 더 높은 니켈 분율($f_{\text{Ni}}$)을 추론한다. 반면 $M_{\text{Ni}}$는 상대적으로 안정적으로 유지된다(고정된 불투명도에서 약 3% 이내로 회복됨).
   • 만약 불투명도가 가변적으로 허용된다면, 편향은 훨씬 낮은 유효 불투명도를 추론하는 방향으로 이동한다($\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$ for fully mixed cases).
   • 이러한 편향은 폭발 에너지 및 프로제니터 매핑에 대한 잘못된 추론으로 이어진다.

3. SN 1998bw 사례 연구:

   • SN 1998bw에 snmix를 적용한 결과 중간 정도의 혼합 분율($f_{\text{mix}} \approx 0.59$)을 얻었다. 단일 구역 최대 가능도 추정치(maximum likelihood estimate)가 사후 분포(posterior)와 겹치기는 하지만, $M_{\text{ej}}$에서는 약 13%, $E_k$에서는 약 5% 차이를 보였다. 이는 볼로메트릭 데이터만으로는 질량, 에너지, 불투명도, 혼합을 유일하게 분해할 수 없음을 보여준다.

4. 콜랩사 r-과정 징후:

   • r-과정 농축의 징후는 보편적으로 후기 근적외선(NIR) 과잉(excess)인 것은 아니다.
   • 배치(Placement): 묻혀 있는 r-과정 물질은 높은 불투명도로 인해 초기에는 광학적 빛을 억제하지만, 농축된 층이 투명해진 후에야 NIR 과잉이 나타난다.
   • 기하학적 구조: 디스크-윈드 모델에서, 만약 중성자 과잉 물질이 적도 평면에 국한되어 있다면, 축 방향 GRB-SNe(극 방향 시선)는 극 방향의 SN 성분이 지배적이므로 r-과정 징후가 거의 나타나지 않는다. 신호는 오프-액시스(off-axis) 관측자에게 가장 강하거나, 물질이 높은 위도($\gtrsim 30^\circ$)까지 혼합되었을 때 가장 뚜렷하다.
   • 진단: 색 진화($\Delta(r-H)$)는 가치 있는 진단 도구이지만, 광학적 억제와 H-밴드 밝아짐을 분리하는 과정이 필요하다. 쌍을 이룬 광학 및 NIR 모니터링(10~60일)이 후기 H-밴드 점 데이터만 사용하는 것보다 더 정보력이 높다.

의의
본 논문은 "방사성 구조"(가열 및 불투명도의 방사형 분포)가 초신성 해석에서 중요하지만 종종 간과되는 변수임을 주장한다.

• 니켈 구동 초신성의 경우: 빠르게 상승하는 Ic-BL 및 GRB-SNe를 설명하기 위해 관습적으로 중심 엔진을 호출하는 것에 이의를 제기하며, 외부 혼합이 유효한 대안이 될 수 있음을 시사한다. 또한, 단일 구역 피팅이 "좋은" 통계적 피팅을 제공하면서도 물리적으로는 잘못된 매개변수를 제공할 수 있다고 경고한다.
• 콜랩사의 경우: 단순한 "NIR 과잉" 휴리스틱을 넘어 r-과정 농축 탐색 전략을 정교화한다. 이젝타가 적도에 국한되어 있다면, 축 방향 GRB-SNe에서 r-과정 징후가 검출되지 않는 것이 r-과정 생성의 부재를 의미하지는 않는다는 점을 강조한다.
• 관측 전략: 저자들은 이러한 퇴보성(degeneracy)을 깨기 위해서는 관측값들을 결합해야 한다고 결론짓는다. 즉, 낮은 이젝타 질량, 높은 혼합, 또는 엔진 동력을 구분하기 위해서는 광구 속도, 색 진화, 스펙트럼, 그리고 후기 볼로메트릭 감소를 모두 결합하는 것이 필수적이다.

이 연구는 snmix가 단순화된 준해석적 도구임에도 불구하고, 방사형 구조의 영향을 효과적으로 격리하여 표준적인 단일 구역 해석에 필요한 교정을 제공한다는 점을 강조한다.