Formulating bulk viscosity for neutron star simulations

Abstract

In order to extract the precise physical information encoded in the gravitational and electromagnetic signals from powerful neutron-star merger events, we need to include as much of the relevant physics as possible in our numerical simulations. This presents a severe challenge, given that many of the involved parameters are poorly constrained. In this paper we focus on the role of nuclear reactions. Combining a theoretical discussion with an analysis connecting to state-of-the-art simulations, we outline multiple arguments that lead to a reactive system being described in terms of a bulk viscosity. The results demonstrate that in order to properly account for nuclear reactions, future simulations must be able to handle different regimes where rather different assumptions/approximations are appropriate. We also touch upon the link to models based on the large-eddy-strategy required to capture turbulence.

final version published in Phys. Rev. D

References in corpus (11)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Multi-messenger Observations of a Binary Neutron Star Merger
• Constraining the Maximum Mass of Neutron Stars From Multi-Messenger Observations of GW170817
• Neutron-star radius constraints from GW170817 and future detections
• Kilonovae
• On the importance of viscous dissipation and heat conduction in binary neutron-star mergers
• High resolution numerical-relativity simulations for the merger of binary magnetized neutron stars
• General-Relativistic Large-Eddy Simulations of Binary Neutron Star Mergers
• Large amplitude behavior of the bulk viscosity of dense matter
• Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model
• A covariant approach to relativistic large-eddy simulations: The fibration picture