[논문 리뷰] A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers

David Radice, Sebastiano Bernuzzi외 2명|arXiv (Cornell University)|2021. 11. 29.

Pulsars and Gravitational Waves Research참고 문헌 136인용 수 135

한 줄 요약

이 논문은 중성자별 병합 시뮬레이션에서 처음으로 전체 비선형 도플러 효과와 중성자와의 상호작용을 포함한 새로운 순간 기반 일반 상대론적 중성미자 운반 코드인 THC_M1을 제시한다. 이 코드는 70ms까지 안정적인 장기 3차원 일반 상대론적 시뮬레이션을 가능하게 하며, 디스크에서의 중성미자 조사가 등방성 방사형 조건에서 전자 비율을 증가시키지만, 이는 이전의 방법들에 비해 핵합성 수율에 미치는 영향은 제한적임을 드러낸다.

ABSTRACT

We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in $\upsilon/c$, retaining all nonlinear neutrino-matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino-matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim}70\,$ms. The latter is the longest dynamical spacetime, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the postmerger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries and equations of state.

연구 동기 및 목표

• 일반 상대론에서 비선형 도플러 이동과 전체 중성자-물질 상호작용을 정확히 포착할 수 있는 순간 기반 중성미자 운반 기법을 개발하는 것.
• 최소 70ms까지의 장기 후기병합 시뮬레이션을 처음으로 전체 중성미자 운반을 포함해 수행하는 것.
• 중성미자 운반의 영향을 후기병합 역학, 중력파 방출 및 핵합성 수율에 대해 평가하는 것.
• 새로운 M1 기법과 이전의 M0+누출 및 M1 기법의 결과를 비교하여 전자 비율과 핵합성 수율의 차이를 평가하는 것.

제안 방법

• 코드는 일반 상대론에서 에너지 통합 볼츠만 방정식을 해결하기 위해 순간 기반 접근법을 사용하며, 𝑣/𝑐의 모든 비선형 항을 유지한다.
• 전체 중성자-물질 상호작용 항을 포함하여, 이질적으로 회전하는 병합 잔여물과 같은 상대론적 운동을 하는 매체에서의 중성자 포획을 정확하게 모델링할 수 있다.
• 순서 시스템에 대한 클로처 관계를 사용하며, THC(중성미자 포함 상대론적 유체역학) 코드 프레임워크에 구현되어 있다.
• 안정성과 수렴성을 검증하기 위해 두 가지 클로처(M1 및 M0+누출)와 다양한 해상도로 시뮬레이션을 수행한다.
• 충격파 튜브 및 구형 대칭 문제와 같은 엄격한 테스트 세트를 통해 코드를 검증한다.
• 이론적 기대치와 이전 시뮬레이션과 비교하여 중성미자 빛기와 평균 에너지를 계산한다.

실험 결과

연구 질문

• RQ1 표준 누출 또는 M0 기법과 비교해 전체 중성자-물질 상호작용과 도플러 효과를 포함함으로써 중성자별 병합 잔여물의 역학이 유의미하게 변화하는가?
• RQ2 부피 점성도와 같은 비평형 효과가 후기병합 진동과 중력파 방출에 얼마나 큰 영향을 미치는가?
• RQ3 중앙 잔여물과 순환 디스크에서 오는 중성미자 조사가 역학적 및 천천히 방출되는 물질의 전자 비율과 조성에 어떤 영향을 미치는가?
• RQ4 중성미자 운반은 특히 r-과정 원소의 핵합성 수율에 어떤 영향을 미치는가?
• RQ5 새로운 M1 기법의 결과는 M0+누출 기법과 SXS 협력 결과와 정량적으로 어떻게 비교되는가?

주요 결과

• 새로운 코드는 현재까지 가장 장기적인 3D 일반 상대론적 전체 중성미자 운반 시뮬레이션으로, 병합 후 70ms까지 성공적으로 수행하였다.
• 이질적으로 회전하는 병합 잔여물에서 중성자 포획을 정확히 모델링하기 위해 전체 중성자-물질 상호작용이 필수적임을 입증하였으며, 이는 표준 누출 기법으로는 포착되지 않는다.
• 비선형 효과를 포함한 M1 기법에서도 후기병합 진동의 감쇠를 유도하는 부피 점성도나 기타 비평형 효과에 대한 유의미한 증거는 발견되지 않았다.
• 중성미자 빛기와 평균 에너지는 이론적 기대치와 이전 시뮬레이션과 양호한 일치를 보였으며, M1 결과는 SXS 데이터와 일치하고 M0+누출은 전자 플레어 중성미자 빛기를 최대 2배까지 과대평가하였다.
• M1 기법은 M0+누출보다 방사형에서 더 높은 전자 비율을 예측한다(Ye ≈ 0.4–0.55), 그러나 재처리로 인해 핵합성 수율은 제한적으로 영향을 받는다.
• 디스크에서 오는 중성미자 조사는 등방성 방사형에서 전자 비율을 증가시키며, 에너지 의존성 흡수 단면적 때문에 관측 각도에 따라 달라지는 킬로노바 빛세기 변화를 유도한다.