[논문 리뷰] Non-Spherical Core-Collapse Supernovae II. Late-Time Evolution of Globally Anisotropic Neutrino-Driven Explosions and Implications for SN 1987A

K. Kifonidis, T. Plewa외 3명|ArXiv.org|2005. 11. 11.

Gamma-ray bursts and supernovae참고 문헌 1인용 수 136

한 줄 요약

이 연구는 15 M☉ 비회전성 블루 슈퍼지엄의 핵붕괴 초신성에 대한 2차원 유체역학 시뮬레이션을 제시하며, 전반적으로 비대칭적이며 저차수 모드가 주도하는 폭발(주로 l=2 및 l=1 모드)이 SN 1987A의 주요 특징을 재현함을 보여준다. 이 모델은 높은 속도의 철족 원소 분사(최대 3300 km/s), 수소/헬륨 인터페이스에서의 강한 혼합(수소가 500 km/s까지 혼합됨), 그리고 프로라트(ejecta 비대칭성, 주축/소축 ≈1.6)를 조기에 발생하는 충격파 변형과 반대 충격파 상호작용 억제를 통해 달성하며, 특수 물리학을 요구하지 않는 중성미자 구동 폭발 메커니즘을 지지한다.

ABSTRACT

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l=2 and l=1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**51 erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event. (abbreviated)

연구 동기 및 목표

• 전반적으로 비대칭적이며 저차수 모드가 주도하는 중성미자 구동 폭발이 SN 1987A의 관측된 특징을 재현할 수 있는지 조사하기 위해.
• 이전 모델—특히 고차수 모드 시뮬레이션—이 SN 1987A에서 관측된 대규모 혼합과 고속 분사 물질을 재현하지 못한 실패 원인을 해결하기 위해.
• 초기 충격파 변형과 초기 속도 구조가 고속의 금속 농축 클러스터가 반대 충격파 파괴로부터 보존되는 데 미치는 역할을 조사하기 위해.
• SN 1987A에서 관측된 프로라트 비대칭성과 수소의 내부 혼합이 조기에 발생하는 유체역학적 불안정성에 의해 설명될 수 있는지 평가하기 위해.
• 최소한의 가정으로 중성미자 구동 폭발 이론이 SN 1987A의 관측 데이터를 완전히 설명할 수 있는지 평가하기 위해.

제안 방법

• 핵 붕괴 후 20ms부터 5시간 이상 후까지 15 M☉ 블루 슈퍼지엄 원형성의 2차원 고해상도 유체역학 시뮬레이션을 수행하였다.
• 최근의 볼츠만 운반 시뮬레이션과 일치하는 낮은 초기 중성미자 플럭스와 느린 감쇠를 가진 수정된 중성미자 가열 조건을 적용하여 폭발 스케일을 연장하였다.
• 세포 경계에서의 수치적 플럭스 계산에 대해 AUSM+ 방법을 사용하였으며, 상류 재구성과 임계 음속 조건을 통해 경계면 음속을 계산하여 날카운 충격파 해상도를 확보하였다.
• 충격파 변형과 분사 비대칭성의 진화를 정량화하기 위해 레지오르드 모드 분해(Legendre 모드 분해, l=1 및 l=2)를 사용하였다.
• 유체역학적 불안정성—특히 리히트마이어-메쉬코프 불안정성—의 성장을 분석하기 위해 t ≈ 100 s 경계에서 토크 밀도의 변화를 분석하였다.
• 철족 원소의 속도 및 혼합 패턴을 추적하기 위해 고도의 핵합성 계산을 통합하였다.

실험 결과

연구 질문

• RQ1 중성미자 구동 폭발에서 저차수의 전반적 비대칭 모드(l=1 및 l=2)가 SN 1987A에서 관측된 고속 철족 분사(최대 약 3300 km/s)를 재현할 수 있는가?
• RQ2 초기 충격파 변형은 반대 충격파 상호작용을 지연시켜 고속 금속 농축 클러스터의 생존성과 속도에 어떤 영향을 미치는가?
• RQ3 He/H 경계면에서의 토크 밀도의 변화는 충격파가 구형이 된 후에도 내부 분사의 늦은 시기 비대칭성 생성에 어떤 역할을 하는가?
• RQ4 왜 고차수 모드 모델은 SN 1987A에서 관측된 혼합과 비대칭성을 재현하지 못하며, 이러한 차이를 설명하는 물리적 메커니즘은 무엇인가?
• RQ5 중성미자 구동 폭발 메커니즘만으로도 SN 1987A의 주요 관측 특징을 모두 설명할 수 있으며, 더 논란의 여지가 있는 다른 물리적 이론을 도입할 필요가 없는가?

주요 결과

• 저차수의 전반적 비대칭 폭발은 최대 약 3300 km/s의 철족 분사 속도를 달성하며, SN 1987A의 관측 결과와 일치한다.
• 초기 충격파의 전반적 비대칭성으로 인해 t ≈ 100 s 경계에서 He/H 경계면에 조기에 토크 밀도가 축적되어 강력한 리히트마이어-메쉬코프 불안정성을 유도하며, 이는 늦은 시기의 비대칭성을 유도한다.
• 내부 분사의 프로라트 비대칭성은 t ≈ 10,000 s에 주축/소축 비율 약 1.6에 도달하여 관측된 비대칭성과 일치한다.
• 불안정성에 의해 유도된 혼합으로 인해 수소가 속도 공간에서 최소 약 500 km/s까지 아래로 혼합되며, 이는 SN 1987A의 관측된 조성 구조와 일치한다.
• 금속 농축 클러스터의 초기 최대 속도는 고차수 모델 대비 약 40% 높아, 전체 경로 동안 음속 이하를 유지하며 반대 충격파 상호작용으로 인한 에너지 소모를 피할 수 있었다.
• 시뮬레이션은 중성미자 구동 폭발 메커니즘을 최소한의 가정으로 적용할 경우, SN 1987A의 주요 관측 특징을 완전히 설명할 수 있음을 입증하며, 더 논란의 여지가 있는 물리적 이론을 도입할 필요가 없음을 시사한다.