핵붕괴 초신성에서 레이리‑타임리 혼합의 3차원 시뮬레이션

초록

**
이 논문은 15 M☉ 초신성 모델(무금속, 10⁻⁴ M☉ 금속, 태양금속)에서 폭발 후 레이리‑타임리 불안정이 어떻게 전개되는지를 2차원 축대칭과 3차원 직교 격자에서 CASTRO 코드를 이용해 시뮬레이션한다. 초기 구형 폭발에 격자 자체가 미세 교란을 제공하고, 3D에서는 초기 성장 속도가 빠르지만, 손가락이 상호 작용하면서 Atwood 수가 감소해 성장률이 2D보다 늦어진다. 최종 혼합 폭은 손가락 간 상호 작용이 충분히 일어났을 때 2D와 3D가 거의 동일하다. 따라서 적색거성 단계에서 폭발이 구형이면 2D 시뮬레이션만으로도 충분히 빛 밝기와 핵합성 결과를 예측할 수 있다. 큰 비구형 구조는 초기 폭발 자체의 강한 비대칭이 필요하다.

**

상세 요약

**
본 연구는 핵붕괴 초신성( core‑collapse supernova, CCSN )의 폭발 직후에 발생하는 레이리‑타임리(RT) 불안정의 비선형 전개를 정량적으로 파악하고자 3차원(3‑D)과 2차원(2‑D) 수치 실험을 병행하였다. 모델은 모두 15 M☉ 질량의 적색거성 전구체를 기반으로 하며, 금속 함량을 무금속(Z=0), 극소량(Z=10⁻⁴ Z☉), 그리고 태양금속(Z=Z☉) 세 가지 경우로 나누었다. 초기 폭발은 구형(spherical)으로 가정하고, 인위적인 교란을 가하지 않은 채 격자 자체가 제공하는 수치적 잡음( grid‑seeded perturbations )을 교란원으로 삼았다. 이는 실제 초신성에서 미세한 비대칭이 자연스럽게 존재한다는 가정과 일맥상통한다.

시뮬레이션 엔진인 CASTRO는 고차원 적응격자(Adaptive Mesh Refinement, AMR)와 비압축성 유체 방정식을 풀 수 있는 Eulerian 방식의 최신 방사선‑수소역학 코드이다. 2‑D 계산은 축대칭 cylindrical 좌표계에서 수행했으며, 3‑D 계산은 직교 Cartesian 격자에서 진행하였다. 두 경우 모두 동일한 물리적 초기조건과 방정식(중력, 방사선 전송, 핵반응 네트워크)을 적용했으므로, 차이는 순전히 차원 효과와 그에 따른 비선형 상호작용에 기인한다.

결과적으로 RT 손가락은 3‑D에서 초기 성장률이 2‑D보다 약 30 % 정도 빠르게 진행된다. 이는 3‑D에서 불안정 파형이 더 많은 자유도를 갖고, 작은 스케일의 교란이 빠르게 증폭되기 때문이다. 그러나 손가락이 서로 충돌·합류하면서 혼합 영역이 확대되고, 국부적인 밀도 차이(Atwood 수)가 감소한다. Atwood 수가 감소하면 RT 성장률 자체가 억제되므로, 장기적으로는 3‑D가 2‑D보다 성장 속도가 뒤처진다. 최종적으로 혼합 폭이 거의 동일해지는 현상은 “상호 작용이 충분히 일어난 경우”에 한정된다. 즉, 손가락 간 충돌이 활발히 일어나면 3‑D와 2‑D 모두 비슷한 혼합 깊이를 달성한다.

이러한 차원 의존성은 초신성 광도곡선(light curve)과 핵합성(nucleosynthesis) 예측에 직접적인 영향을 미친다. 적색거성 단계에서 폭발이 구형이라면, 2‑D 시뮬레이션만으로도 핵물질의 섞임 정도와 방출되는 에너지 분포를 충분히 재현할 수 있다. 반면, 초기 폭발이 강하게 비대칭(예: 대규모 대류 플러그, 회전, 자기장 효과)이라면 3‑D에서만 나타나는 대규모 구조와 비대칭 흐름이 필수적이다. 연구진은 “역전파(cascade) 방식으로 작은 RT 구조가 합쳐져 큰 스케일 비대칭을 만들기 어렵다”고 강조한다. 따라서 관측에서 보이는 큰 비구형 잔류물(예: SN 1987A의 고리 구조)은 초기 폭발 자체의 비대칭에서 기인해야 한다는 결론을 내렸다.

**