Problematic systematics in neutron-star merger simulations

Next-generation gravitational-wave instruments are expected to constrain the equation of state of dense nuclear matter by observing binaries involving neutron stars. We highlight a problematic systematic error in finite-temperature merger simulations, where shock heating associated with the neutron-star surface gives rise to elevated temperatures. We demonstrate the severe implications of this artificial heating by computing static and dynamical tidal parameters for neutron stars immersed in simulation temperature profiles. The induced systematic errors must be addressed if we want to build robust gravitational-wave signal models for neutron-star, or indeed neutron star-black hole, binaries.