[논문 리뷰] Self-consistent picture of the mass ejection from a one second-long binary neutron star merger leaving a short-lived remnant in general-relativistic neutrino-radiation magnetohydrodynamic simulation

Kenta Kiuchi, Sho Fujibayashi외 4명|arXiv (Cornell University)|2022. 11. 14.

Gamma-ray bursts and supernovae인용 수 6

한 줄 요약

이 연구는 부드러운 상태방정식(SFHo)를 가진 1초 동안의 이중 중성자별 병합을 일반 상대성 이론, 중성미자 복사 전도를 고려한 자기유체역학(simulation)을 통해 처음으로 수행하였으며, 동역학적 및 난류적 후 병합 질량 분출을 모두 해결하였다. 이는 빠르고 약간 상대성 이론적 속도를 가진 동역학적 분출과 느리고 난류 점성에 의해 유도된 후 병합 분출이라는 두 가지 구분된 분출 성분을 확인하였으며, 후자에서는 r-과정의 임계 값에 가까운 전자 분포를 보여, 희토류 원소의 효율적 생성 가능성을 시사한다.

ABSTRACT

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on Japanese supercomputer Fugaku using about $72$ million CPU hours with $20,736$ CPUs. We consider an asymmetric binary neutron star merger with masses of $1.2$ and $1.5M_\odot$ and a `soft' equation of state SFHo. It results in a short-lived remnant with the lifetime of $\approx 0.017$\,s, and subsequent massive torus formation with the mass of $\approx 0.05M_\odot$ after the remnant collapses to a black hole. For the first time, we confirm that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the post-merger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

연구 동기 및 목표

• 자기 일관성 있는 중성미자 복사 및 자기유체역학을 포함한 1초 동안의 전체 이중 중성자별 병합 시뮬레이션 수행.
• 동역학적 단계와 후 병합 토러스에서의 질량 분출 과정을 해석하고, 특히 자기회전 불안정성(MRI)이 분출을 어떻게 유도하는지 규명.
• 분출의 핵합성 영향을 조사하며, 특히 전자 분포의 분포와 r-과정 원소 생성과의 관련성 분석.
• GW170817과 GRB 170817A 사이의 관측된 1.7초 지연을 고려할 때, 상대성 속도를 가진 제트가 어떻게 발생할 수 있는지 조사.
• 시뮬레이션 지속 시간, 해상도, BH 스핀이 Poynting-플럭스 지배적 외부 유동의 발생에 미치는 영향 평가.

제안 방법

• Fugaku 슈퍼컴퓨터를 사용하여 20,736개의 CPU와 약 8500만 CPU 시간을 사용해 일반 상대성 이론, 중성미자 복사 자기유체역학(GRMHD) 시뮬레이션 수행.
• 중성자별 물질에 대해 부드러운 상태방정식(SFHo)을 사용하고, 질량이 1.2 및 1.5 M☉인 비대칭 이중 병합을 시뮬레이션.
• 분출 성분의 전자 분포, 에ント로피, 최종 속도 분포를 계산하기 위해 트레이서 입자를 사용해 질량 분출 추적.
• 난류 점성과 MRI 성장의 해상도를 확보하기 위해 유한체적, 고해상도 유한체적 방법과 적응 메쉬 구역 강화 기법 사용.
• 에너지 및 에ント로피 이동을 모델링하기 위해 단순화된 중성미자 수송 기법을 사용해 중성미자 빛기 및 에너지 플럭스 계산.
• 분출 시점에 따라 분출 성분 정의: 동역학적 분출(~0.1초 이내), 후 병합 분출(0.1초에서 BH 형성 시점인 약 0.017초까지).

실험 결과

연구 질문

• RQ1 이중 중성자별 병합 중 및 병합 후 질량 분출을 유도하는 독립적인 물리적 메커니즘은 무엇이며, 한 번의 장기 시뮬레이션으로 이를 해결할 수 있는가?
• RQ2 실제적인 GRMHD 조건에서 MRI에 의해 유도된 난류 점성이 후 병합 분출의 형성과 특성에 어떻게 기여하는가?
• RQ3 후 병합 분출의 전자 분포는 어떻게 되며, 이는 효율적인 r-과정 핵합성에 기여하는가?
• RQ4 왜 시뮬레이션에서 Poynting-플럭스 지배적 외부 유동이 관측되지 않았는가? BH 스핀, 격자 해상도, 재분포 역학이 제트 형성에 미치는 영향은 무엇인가?
• RQ5 모의 분출 특성(질량, 속도, 조성)이 AT 2017gfo의 관측 제약 조건을 얼마나 잘 충족하는가?

주요 결과

• 시뮬레이션은 약 0.017초 동안 지속되는 초질량 중성자별 잔여물이 붕괴하여 BH로 전환되는 것을 해결하였다.
• 두 가지 구분된 분출 성분이 확인되었으며, 질량은 각각 약 6×10⁻³ M☉인 동역학적 분출과 MRI에 의해 유도된 난류 점성에 의해 유도된 후 병합 분출(약 8×10⁻³ M☉)이다.
• 후 병합 분출의 전자 분포 Yₑ는 피크 시점에서 약 0.26–0.27로, 란타니드-프리 r-과정 핵합성의 임계 값(약 0.25)에 매우 가까운 것으로 나타났다.
• BH 붕괴 후 형성된 토러스 질량은 약 0.002 M☉이며, 병합 후 약 1.1초에 총 분출 질량은 약 0.014 M☉로 추정되었다.
• 중력파가 발생한 후 Poynting-플럭스 지배적 외부 유동은 관측되지 않았으며, 이는 중간 정도의 BH 스핀, 재분포 압력, 짧은 시뮬레이션 지속 시간, 해상도 제한 등의 영향일 가능성이 높다.
• 이 시뮬레이션 결과는 SFHo 상태방정식가 관측된 AT 2017gfo의 단순한 킬로노바 광선 곡선 모델링에서 요구하는 질량(Mₑ𝗂𝗂 ≈ 0.05 M☉)보다 분출 질량이 부족하므로, SFHo 상태방정식가 불리할 수 있음을 시사한다.