[논문 리뷰] Dynamical ejecta of neutron star mergers with nucleonic weak processes I: Nucleosynthesis

I Kullmann, S. Goriely외 4명|arXiv (Cornell University)|2021. 09. 06.

Pulsars and Gravitational Waves Research참고 문헌 101인용 수 52

한 줄 요약

이 연구는 중성자별 병합 시의 동역학적 분출물에서 핵자성 약한 상호작용—특히 중성미자 상호작용—이 r-과정 핵합성과 방사성 가열에 미치는 영향을 조사한다. 자체 일관된 중성미자 수송(ILEAS 프레임워크)을 적용한 고급 유체역학 시뮬레이션을 통해, 중성미자에 의한 전자 분율(Ye) 변화가 A ≈ 90까지 태양 유사 r-과정 생성분포를 만들어내며, 방사성 우라늄 및 플루토늄 생성은 감소하고 스트론티움(Sr) 농도는 증가함을 입증한다. 이 결과는 상태방정식 또는 질량 비대칭성과 무관하게 유지된다.

ABSTRACT

We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r- process elements with mass numbers $A \gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the dynamically ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations, at least if the dynamical component dominates the overall ejecta. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models.

연구 동기 및 목표

• 중성자별 병합의 동역학적 분출물에서 핵자성 약한 상호작용(중성미자 흡수 및 방출)이 r-과정 핵합성에 미치는 영향을 평가하기 위해.
• 약한 상호작용을 무시한 시뮬레이션과 비교하여 중성미자 상호작용이 전자 분율(Ye)을 어떻게 변화시키고, 그로 인해 분출물의 최종 조성에 어떤 영향을 미치는지 규명하기 위해.
• 다양한 상태방정식(DD2, SFHo)과 질량 비대칭성(1.35–1.35 M⊙ 대비 1.25–1.45 M⊙)에 대해 핵합성 수확량과 가열률의 강인성을 평가하기 위해.
• 킬로노바 광선곡선과 관측 데이터로부터 분출물 성질의 탐지 가능성에 미치는 영향을 검토하기 위해.
• 후기 가열과 복사에 기여하는 자유 중성자와 분열 사슬의 역할을 평가하기 위해.

제안 방법

• 다양한 질량 비율과 상태방정식(DD2 및 SFHo)을 가진 중성자별-중성자별 병합을 대상으로 한 3D, 일반상대론적 유체역학 시뮬레이션 4건을 활용한다.
• 강성한 물질을 고려한 ILEAS(암시적 라그랑주-오일러 방법) 프레임워크를 적용하여 중성미자-물질 상호작용을 모델링하며, 밀도가 높은 영역에서는 중성미자 평형을, 밀도가 낮은 영역에서는 재흡수 효과를 고려한다.
• 약한 상호작용을 통한 전자 분율(Ye) 시간적 변화를 추적한다: νe + n ⇌ p + e⁻ 및 ḡνe + p ⇌ n + e⁺.
• 유체역학에서 확보한 최종 열역학 경로를 기반으로 후처리 핵합성 계산을 수행하여 동위원소 농도를 산출한다.
• r-과정 뉴클리드의 방사성 붕괴 가열률을 계산하며, β 붕괴, 자발적 분열(예: ²⁵⁴Cf), 및 α 붕괴 사슬 기여도를 포함한다.
• 초기 조건, 특히 질량 비대칭성과 상태방정식에 따른 핵합성 및 가열의 의존도를 분석한다.

실험 결과

연구 질문

• RQ1 핵자성 약한 상호작용—특히 중성미자 흡수 및 방출—은 중성자별 병합에서 동적으로 분출된 물질의 전자 분율(Ye)에 어떤 영향을 미치는가?
• RQ2 그 결과로 나타나는 r-과정 핵합성 수확량은 무엇인가? 특히 태양 유사 농도 패턴, 스트론티움(Sr), 랜타니드 및 악티니다이드의 생성 여부는?
• RQ3 이상의 최종 농도와 가열률은 이분화 상태방정식 또는 질량 비대칭성의 변화에 대해 어느 정도 민감하지 않게 유지되는가?
• RQ4 약한 상호작용을 포함함으로써 분열 및 α 붕괴 사슬이 후기 가열에 기여하는 정도는 어떻게 변화하는가?
• RQ5 잔류 자유 중성자와 약한 상호작용에 의해 수정된 조성으로 인해 나타나는 관측 가능한 서명(예: 광선곡선 특징, X선 방출)은 무엇인가?

주요 결과

• 핵자성 약한 상호작용을 포함함으로써 동역학적 분출물의 전자 분율(Ye)이 크게 증가하여, 극 지역에서 <0.1에서 ∼0.3–0.4 수준으로 상승한다. 이는 주로 중성미자 흡수에 기인한다.
• 그 결과로 나타나는 r-과정 핵합성은 질량수 A ≥ 90인 원소에 대해 태양 유사 농도 패턴을 생성하며, 특히 스트론티움(Sr) 농도가著 증가한다. 이는 약한 상호작용을 무시한 시뮬레이션과는 대조된다.
• 중성미자 상호작용이 없는 모델에 비해 악티니다이드(예: uraninite 및 californium) 생성이 크게 감소한다. 이는 중성자 농도 감소로 인한 것이다.
• 초기 조건의 차이에도 불구하고, 네 모델 전부에서 최종 농도 패턴과 가열률이 놀랄 정도로 유사하여, 상태방정식과 질량 비대칭성에 대해 약한 탈감도를 보임을 시사한다.
• 가열률은 병합 후 약 10분, 7시간, 30일 경과 시점에서 초기 Ye 값에 가장 민감하며, 주로 잔류 r-뉴클리드의 β 붕괴에 기인한다. 약한 상호작용을 포함할 경우 자발적 분열 및 α 붕괴 사슬의 기여도는 감소한다.
• 약한 상호작용이 포함된 상태에서도 고속으로 확장되는 분출층은 여전히 무시할 수 없는 비율의 자유 중성자를 유지하며, 이는 후기 X선 방출에 기여하고, 킬로노바 광선곡선에 탐지 가능한 서명을 남길 수 있다.