Identifying thermal effects in neutron star merger remnants with model-agnostic waveform reconstructions and third-generation detectors

We explore the prospects for identifying differences in simulated gravitational-wave signals of binary neutron star (BNS) mergers associated with the way thermal effects are incorporated in the numerical-relativity modelling. We consider a hybrid approach in which the equation of state (EoS) comprises a cold, zero temperature, piecewise-polytropic part and a thermal part described by an ideal gas, and a tabulated approach based on self-consistent, microphysical, finite-temperature EoS. We use time-domain waveforms corresponding to BNS merger simulations with four different EoS. Those are injected into Gaussian noise given by the sensitivity of the third-generation detector Einstein Telescope and reconstructed using BayesWave, a Bayesian data-analysis algorithm that recovers the signals through a model-agnostic approach. The two representations of thermal effects result in frequency shifts of the dominant peaks in the spectra of the post-merger signals, for both the quadrupole fundamental mode and the late-time inertial modes. For some of the EoS investigated those differences are large enough to be told apart, especially in the early post-merger phase when the signal amplitude is the loudest.