Eccentricity reduction of binary neutron star initial data with the entropy based flux limiting scheme

TL;DR

This paper extends an entropy-based flux limiting scheme to construct low-eccentricity initial data for binary neutron star simulations, achieving high-order convergence in gravitational waveform phase.

Contribution

It introduces the use of the EFL method for initial data construction, enabling more accurate and efficient binary neutron star merger simulations.

Findings

Achieved up to fifth-order convergence in gravitational waveform phase.

Produced several new binary neutron star simulations with eccentricity reduction.

Demonstrated the effectiveness of EFL in initial data setup.

Abstract

The construction of high-resolution shock-capturing schemes is vital in producing highly accurate gravitational waveforms from neutron star binaries. The entropy based flux limiting (EFL) scheme is able to perform fast converging binary neutron star merger simulations reaching up to fourth-order convergence in the gravitational waveform phase. In these results the EFL method was used only in the dynamical evolution of initial data constructed with the Lorene library. Here, we extend the use of the EFL method to the construction of eccentricity reduced initial data for neutron star binaries and present several new BNS simulations resulting from such initial data and show for the first time up to optimal fifth-order convergence in the gravitational waveform phase.