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

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 binary neutron star simulations resulting from such initial data and show for the first time up to optimal fifth-order convergence in the gravitational waveform phase.