Revisiting the Constraint on the Equation of State of Neutron Stars Based on Binary Neutron Star Mergers

The merger of a neutron star (NS)–NS binary can form different productions of compact remnants, among which a supramassive NS (SMNS) can create an internal plateau, and the following steep decay marks the collapse of the SMNS. The proportion of the SMNS and the corresponding collapse time are often used to constrain the NS equation of state (EOS). This paper revisits this topic by considering the effect of an accretion disk on a compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks ∼100 s) of SMNSs formed from an NS–NS merger is almost unaffected by the initial surface magnetic field (Bs,i) of the NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on Bs,i significantly. Coupling the constraints from the SMNS fraction, we exclude some EOSs and obtain three candidate EOSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse time and the luminosity of the internal plateau (in the short gamma-ray bursts) for observations obtained based on the three candidate EOSs, it is shown that only the EOS of ENG is favored. Our sample, based on the ENG EOS and a mass distribution motivated by Galactic systems, suggests that approximately 99% of NS–NS mergers collapse to form a black hole (BH) within 107s. This includes scenarios promptly forming a BH (36.5%), an SMNS (60.7%), or a stable NS that transitions into a BH or an SMNS following accretion (2.1%). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814, are more likely to be BHs.