Entropy-limited higher-order central scheme for neutron star merger simulations

Numerical relativity simulations are the only way to calculate exact gravitational waveforms from binary neutron star mergers and design templates for gravitational-wave astronomy. The accuracy of these numerical calculations is critical in quantifying tidal effects near merger that currently one main sources uncertainty waveforms. In this work, we explore use an entropy-based flux-limiting scheme high-order, convergent spacetimes. effectively tracks stellar surface physical shocks using residual entropy equation thus allowing unlimited central flux schemes regions smooth flow. We perform first with such a method demonstrate up fourth-order convergence waveform phase. reduces phase error factor five when compared state-of-the-art high-order characteristic can be employed producing faithful modelling.