Emission line models for the lowest mass core-collapse supernovae -- II. 3D NLTE radiative transfer modelling of a $9.0\,M_\odot$ neutrino-driven explosion

기술 요약: 최저 질량 핵붕괴 초신성을 위한 방출선 모델 – II. 9.0 $M_\odot$ 중성미자 구동 폭발에 대한 3D NLTE 복사 전달 모델링

문제 및 동기
폭발 후 수개월에서 수년이 경과한 핵붕괴 초신성(CCSNe)의 성운기(nebular phase)는 방출선의 프로파일에 폭발의 비대칭성을 기록한다. 정확한 예측을 위해서는 비국소 열역학적 평형(NLTE) 복사 전달이 필수적이지만, 저질량 전구체($M_{\text{ZAMS}} \approx 8-12 M_\odot$)에 대한 기존 연구들은 주로 1D 모델에 의존해 왔다. 이러한 1D 모델은 인위적인 껍질 구조를 생성하며, 중성미자 구동 폭발에 내재된 복잡한 3D 혼합과 비대칭성을 포착할 수 없다. 특히, 철 핵(iron-core) CCSNe와 전자 포획 초신성(ECSNe)을 구분하는 진단적 특징들이 3D 기하 구조에서도 유지되는지, 그리고 저에너지 폭발에서 관측 각도 변화가 선 광도와 프로파일에 어떤 영향을 미치는지 여부는 여전히 불확별하다. 본 연구는 9.0 $M_\odot$ 전구체의 수력학 시뮬레이션에 업그레이드된 3D NLTE 복사 전달 코드를 적용함으로써 이러한 공백을 메우고자 한다.

방법론
저자들은 최근 1D/2D 프레임워크에서 완전한 3D NLTE 복사 전달 솔버로 업그레이드된 ExTraSS(EXplosive TRAnsient Spectral Simulator) 코드를 활용한다. 주요 기술적 업그레이드는 다음과 같다:

• 복사 전달: 이전의 "on-the-spot" 근사법을 대체하여, 복사장에서 유도된 완전한 광이온화 및 광들뜸 비율을 구현하였다. 코드는 메모리 부하를 관리하기 위해 도메인 분할(Brunner & Brantley 2009에서 영감을 얻음)을 사용하는 광선 추적(ray-tracing) 기법을 채택하며, 방위각 슬라이스에 따라 계산을 병렬화하기 위해 $\phi$ 기반 분할을 사용한다.
• 미세물리(Microphysics): 비열적 들뜸을 위한 에너지 침적을 계산하고, H I, O I, Mg I, Fe I, Fe II에 대한 특정 재결합 비율을 사용하며, 업데이트된 광이온화 단면적(Verner et al. 1996)을 포함한다.
• 폭발 모델: 본 연구는 9.0 $M_\odot$ 태양 금속량 전구체의 3D 중성미자 구동 폭발인 s9.0 모델(Stockinger et al. 2020)을 사용한다. 시뮬레이션은 탄성 충격파 발생(bounce) 후 19.74일까지 진화되었으며, 이 시점의 이젝타(ejecta)를 성운기(250일 및 400일) 스펙트럼 합성을 위해 동종 팽창 격자($38 \times 15 \times 30$ 셀)로 매핑하였다.
• 비교: 합성 스펙트럼은 400개의 관측 각도($20 \times 20$ 그리드)에 대해 생성되었으며, SN 1997D 및 SN 2016bkv의 관측값과 비교되었다. 결과는 동일한 전구체에 대한 Jerkstrand et al. (2018, J18)의 1D 모델과 벤치마킹되었다.

주요 결과

1. 물리적 조건: 3D 모델은 1D J18 모델에 비해 더 균일한 온도 구조를 보인다. 3D에서의 인위적인 껍질 층상 구조 결여는 이젝타 전체에 걸쳐 더 분산된 $\gamma$-선 에너지 침적을 유도하며, 이는 더 높은 코어 밀도와 더 효율적인 재결합으로 인해 더 낮은 전체 트래핑 효율(400일째에 1D 침적의 82%)과 더 낮은 전자 분율($x_e$)을 초래한다.
2. 스펙트럼 진화 및 원소 분포: 3D 이젝타는 강한 혼합을 보이며, 주요 $^{56}$Ni 플룸(plume)은 최대 1400 km s$^{-1}$의 속도까지 확장된다. 이러한 혼합은 1D 대응 모델에 비해 3D 모델에서 Fe I 방출을 크게 강화하며, 특히 6200–6700 Å 및 7900–8400 Å 영역에서 그러하다.
3. 관측과의 비교 (SN 1997D 및 SN 2016bkv):
   • 모델은 SN 1997D에 대해 He, C, O, Mg 선의 프로파일을 상당히 잘 재현한다.
   • 그러나 모델은 SN 1997D에 비해 H$\alpha$ 및 Fe I 선의 세기를 과대 예측한다.
   • SN 2016bkv의 경우, 모델은 7900–8400 Å 영역에서 강한 Fe I 선이 나타나지 않는 현상을 재현하지 못하며, 이는 해당 천체에 대해서도 모델의 Fe I 방출이 너무 강함을 시사한다.
4. 관측 각도 효과: 선 프로파일과 광도는 관측 각도에 따라 크게 변한다.
   • H$\alpha$: 가장 밝은 프로파일은 청색 편이(blueshifted)되어 있으며, 가장 어두운 프로파일은 적색 편이(redshifted)되고 더 넓다. 가장 밝은 관측과 가장 어두운 관측 사이의 광도 차이는 약 16%이다.
   • Mg I]: 짧은 파장에서의 높은 불투명도로 인해 훨씬 더 큰 광도 변화(2배 이상의 인자)를 보인다.
   • 중심 이동(Centroid Shifts): 수력학적 킥(kick)이 중심 이동을 결정하는 고(high)-킥 모델과 달리, 이 모델의 낮은 킥 속도(중성미자에 의한 비대칭성에 의해 지배됨)는 중성자별 킥 벡터와 엄격하게 일치하지 않는 복잡한 중심 이동 패턴을 결과한다.
5. 진단적 특징: 철 CCSNe와 ECSNe를 구분하는 핵심 진단선들(예: He I $\lambda$7065, [C I] $\lambda$8727, [C I] $\lambda\lambda$9824, 9850)은 3D 모델에서도 여전히 존재하며, 이는 3D 혼합이 존재하는 상황에서도 1D 기반 진단법의 유효성을 뒷받침한다.

의의 및 주장
본 논문은 저질량 CCSNe의 성운기를 정확하게 모델링하기 위해서는 3D NLTE 복사 전달이 필수적임을 입증한다. 1D 모델은 방사성 니켈의 혼합과 그로 인한 선 프로파일의 영향을 포착하지 못하기 때문이다. 저자들은 다음을 보여준다:

• 3D 진단법의 유효성: 저질량 철 CCSNe를 ECSNe와 구분하는 주요 스펙트럼 특징들은 3D 혼합 효과에도 불구하고 견고하게 유지된다.
• 비대칭성 탐지: 선 프로파일의 관측 각도에 대한 강한 의존성은 폭발의 비대칭성을 추론하기 위한 상세한 비교를 가능하게 한다.
• $^{56}$Ni 추적: 가장 빠른 $^{56}$Ni 이젝타(최대 1400 km s$^{-1}$)조차 성운기 선 프로파일에 추적 가능한 흔적을 남긴다.
• 모델의 한계: 현재의 s9.0 모델은 SN 1997D 및 SN 2016bkv의 관측치에 비해 Fe I 및 H$\alpha$ 선을 너무 강하게 생성하는 경학년을 보이는데, 이는 3D 프레임워크가 상당한 개선임에도 불구하고 특정 저광도 사건들에 맞추기 위해서는 특정 폭발 파라미터나 미세물리적 정교화가 여전히 필요함을 나타낸다.

이 연구는 저질량 전구체를 위한 1D에서 3D 모델링으로의 경로를 구축하며, 저광도 Type IIP 초신성의 향후 성운기 관측을 해석하기 위한 더 현실적인 프레임워크를 제공한다.