Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Francois Foucart,Patrick Chi-Kit Cheong,Matthew D. Duez,Lawrence E. Kidder,Harald P. Pfeiffer,Mark A. Scheel
2024-09-24
Abstract:Binary neutron star mergers play an important role in nuclear astrophysics: their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals, and the relative abundances of the produced r-process elements. Existing merger simulations are however limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte-Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant ($(5-10)\,{\rm ms}$). We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates, transport algorithm, as well as recent improvements by Radice et al to the two-moment algorithms can all contribute to changes at the $(10-30)\%$ level in the global properties of the merger remnant and outflows. The limitations of previous moment schemes fixed by Radice et al also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte-Carlo and moment schemes in the context of a low mass binary neutron star system.
High Energy Astrophysical Phenomena,General Relativity and Quantum Cosmology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the influence of different neutrino transport methods and neutrino - matter interaction rates on the merger results in neutron - star merger simulations. Specifically, the researchers focus on the impact of using advanced two - moment or Monte - Carlo transport schemes on the behavior of high - density regions, the collapse time, and the properties of hot tidal arms during the merger process. In addition, they also explore the application effects of different reaction rates in the simulation, and how the recent improvement of the two - moment algorithm by Radice et al. affects the overall properties of the merger remnant and the outflow. ### Background of the paper Neutron - star mergers are important events in nuclear astrophysics. The gravitational waves and electromagnetic signals they generate carry information about the equation of state of cold matter above the nuclear saturation density, and are also one of the main sources of r - process elements in the universe. Neutrino - matter interactions have a significant impact on the properties of the outflow and the relative abundances of r - process elements during and after the merger process. However, existing merger simulations have limitations in realistically modeling neutrino transport and neutrino - matter interactions. ### Research purposes This study evaluates the differences between these methods by comparing the effects of different neutrino transport methods (such as two - moment and Monte - Carlo transport schemes) on a single short - lived neutron - star merger system (with a duration of 5 - 10 milliseconds). The study also examines the impact of the choice of different reaction rates on the simulation results. ### Main findings 1. **Influence of transport methods**: - Using the most advanced two - moment or Monte - Carlo transport schemes, different transport methods have significant differences in the behavior of high - density regions, which will affect the collapse time of the system and the properties of hot tidal arms. - The best transport schemes are qualitatively basically the same in terms of the influence of neutrinos on the system, but there are 10% - 30% changes in global properties. 2. **Influence of reaction rates**: - Different choices of reaction rates, transport algorithms, and the improvement of the two - moment algorithm by Radice et al. will all lead to changes in the global properties of the merger remnant and the outflow. - These changes are sufficient to explain the significant differences in the production of heavy - lepton neutrinos previously observed in low - mass double - neutron - star systems. ### Methods - **Initial conditions**: All simulations consider the merger of two neutron stars (with masses of 1.4 solar masses and 1.3 solar masses respectively), in an initial non - rotating state, and evolve from about 6 orbits before the merger. - **GR - Hydro method**: Use the hydrodynamic method under general relativity, including the generalized harmonic form of Einstein's equations and the fluid equations in Valencia form. - **Neutrino transport**: Compare the two - moment transport scheme (including Minerbo closure) and the Monte - Carlo transport scheme, as well as the improved method proposed by Radice et al. ### Conclusions Through these studies, the authors hope to improve the understanding of neutrino transport and neutrino - matter interactions during the neutron - star merger process, so as to better explain the observed gravitational waves and electromagnetic signals, and further understand the production mechanism of r - process elements.