Simulating Binary Neutron Star Mergers with Finite-temperature Equations of State: The influences of the slope of the symmetry energy and artificial heating

We present a new set of numerical-relativity simulations of merging binary neutron stars, aiming to identify possible observable signatures of the slope of the symmetry energy Lsym. To achieve this goal, we employ a set of equations of state based on a parameterization of the covariant density functional theory of nuclear matter that allows controlled variations of Lsym and the skewness Qₛₐₜ, holding the latter fixed. For a set of our simulations, we identify a steep energy gradient in the equation of state at subsaturation densities, which acts as a source of heating with subsequent stiffening produced by thermal support. Accounting for related structural modifications in the tidal deformability reconciles our results with theoretical expectations. On the other hand, we show that gravitational waves are unlikely to distinguish the role of Lsym. In contrast to this, we find that the ejecta composition is significantly altered in our simulations, which employ an M1 moment scheme, when Lsym is varied. Based on our extracted dynamical ejecta properties, we compute r-process yields and find that they are distinct for the different Lsym, especially at lower mass numbers A 120. This suggests that electromagnetic counterparts are more likely to exhibit signatures; however, a direct connection to Lsym remains a challenge, given the complex interplay between details of the ejecta properties and the kilonova signal.