Draft version July 21, 2021 Typeset using LATEX twocolumn style in AASTeX63 Survival times of supramassive neutron stars resulting from binary neutron star mergers Paz Beniamini1, 2 and Wenbin Lu1, 3 1Theoretical Astrophysics, Walter Burke Institute for Theoretical Physics, Mail Code 350-17, Caltech, Pasadena, CA 91125, USA 2Astrophysics Research Center of the Open University (ARCO), The Open University of Israel, P.O Box 808, Ra'anana 43537, Israel 3Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA ABSTRACT A binary neutron star (BNS) merger can lead to various outcomes, from indefinitely stable neutron stars, through supramassive (SMNS) or hypermassive (HMNS) neutron stars supported only temporar- ily against gravity, to black holes formed promptly after the merger. Up-to-date constraints on the BNS total mass and the neutron star equation of state suggest that a long-lived SMNS may form in ∼ 0:45 − 0:9 of BNS mergers. A maximally rotating SMNS needs to lose ∼ 3 − 6 × 1052 erg of it's rotational energy before it collapses, on a fraction of the spin-down timescale. A SMNS formation imprints on the electromagnetic counterparts to the BNS merger. However, a comparison with ob- servations reveals tensions. First, the distribution of collapse times is too wide and that of released energies too narrow (and the energy itself too large) to explain the observed distributions of internal X-ray plateaus, invoked as evidence for SMNS-powered energy injection. Secondly, the immense energy injection into the blastwave should lead to extremely bright radio transients which previous studies found to be inconsistent with deep radio observations of short gamma-ray bursts. Furthermore, we show that upcoming all-sky radio surveys will constrain the extracted energy distribution, indepen- dently of a GRB jet formation. Our results can be self-consistently understood, provided that most BNS merger remnants collapse shortly after formation (even if their masses are low enough to allow for SMNS formation). This naturally occurs if the remnant retains half or less of its initial energy by the time it enters solid body rotation. Keywords: gamma-ray burst: general { gravitational waves { stars: jets { stars: neutron 1. INTRODUCTION longer than the rotational period of ∼ 1 ms, where the The electromagnetic (EM) appearance of binary neu- Shakura & Sunyaev(1973) dimensionless viscosity pa- −2 tron star (BNS) mergers depend strongly on the nature rameter is found to be α & 10 in the outer envelope and evolution of the post-merger remnant, which in turn of the remnant NS as given by magneto-hydrodynamic depends on the component masses as well as the equa- (MHD) turbulence (Kiuchi et al. 2018). However, sim- tion of state (EoS) of matter in the deep interior of neu- ulations have not been able to capture the full range tron stars (NSs) (Shibata & Hotokezaka 2019). Numer- of lengthscales and physical processes needed to un- ical simulations of the gravitational wave (GW) driven derstand the transport of angular momentum by e.g., dynamical merger process show that the initial remnant Kelvin-Helmholtz and magneto-rotational instabilities arXiv:2104.01181v2 [astro-ph.HE] 20 Jul 2021 is supported against gravitational collapse by strong dif- (e.g., Kiuchi et al. 2014; Ciolfi et al. 2019). A sim- ferential rotation and partially by thermal pressure (Shi- pler approach is to carry out (two-dimensional) viscous bata & Taniguchi 2006; Baiotti et al. 2008; Sekiguchi hydrodynamic simulations, under the assumption that et al. 2011; Hotokezaka et al. 2013; Kiuchi et al. 2014; sub-grid MHD processes operate efficiently to generate a Kaplan et al. 2014). Subsequently, the system evolves macroscopic viscosity (e.g., Shibata et al. 2017b; Radice −1 2017; Fujibayashi et al. 2018). It is found that the outer on the viscous timescale that is roughly α . 100 times envelope viscously spreads into a torus, which contains ∼ 0:1M and a large fraction of the angular momen- Corresponding author: Paz Beniamini tum. However, the evolution of the rotational profile of [email protected] the neutron star (e.g., how it approaches uniform rota- tion) and the share of angular momentum between the 2 Beniamini & Lu NS and the torus are highly uncertain as they depend et al. 2010; Metzger et al. 2011; Dall'Osso et al. 2011). sensitively on the choice of viscosity prescriptions. On On the other hand, the ejecta acquires a large kinetic en- longer timescales ∼ 1 s, neutrino cooling removes ther- ergy comparable to the rotational energy of SMNS, and 2 mal energy (∼ 0:05M c ), increases the central density when decelerated by the surrounding medium, it pro- of the NS, and may cause marginally stable systems to duces bright multi-band afterglow emission (Gao et al. collapse. If at the end of this phase the NS has a suf- 2013; Metzger & Bower 2014). ficiently low mass and sufficiently high angular momen- In this paper, we study the distributions of the sur- tum, it will likely settle into uniform rotation, and be- vival time and the emitted energy from SMNSs prior to come a supramassive neutron star (SMNS), supported the collapse, crucial ingredients for determining the EM against gravity by its fast rotation. counterparts of such systems (see e.g. Ravi & Lasky Observations of NSs near two solar masses in a number 2014). The goal is to make predictions based on known of sources (Demorest et al. 2010; Antoniadis et al. 2013; information (Galactic BNS statistics, LIGO observa- Cromartie et al. 2020) suggest that the pressure above tions, current EoS constraints; §2,3) and then compare nuclear saturation density of 2:8 × 1014 g cm−3 must be with observations (§4). We will show that these dif- sufficiently high such that the maximum mass of non- ferent components, brought together, reveal a tension. 1 rotating NSs, Mmax, is greater than 2M . Centrifugal The implications which could provide insights onto the support due to uniform rotation allows the maximum early stages of the merger remnant's evolution and its mass to be up to about 20% higher than that of the non- stability are discussed in §5. rotating spherical configuration (e.g., Cook et al. 1994; Breu & Rezzolla 2016). This means that a fraction of 2. METHOD 2 BNS merger remnants, at least the ones . 2:4M (and possibly even higher mass ones), could be spinning suf- As we are interested in systems surviving for time- ficiently rapidly as SMNSs, which undergo secular spin- scales longer than the GRB prompt duration we focus on down on longer timescales 1 s. cold and uniformly rotating NSs, assuming that the dif- Magnetic spin-down from a long-lived ( 1 s) SMNS ferential rotation has subsided (on timescale of ∼ 0:1 s) or stable NS provides a powerful source of energy injec- and neutrino cooling has ended (on timescale of ∼ 1 s). tion, which has an important impact on the EM coun- We apply realistic equations of state (EOS) using the 3 terparts of BNS mergers. For instance, the baryonic rns code (Stergioulas & Friedman 1995) to simulate 4 ejecta may be strongly heated by the non-thermal ra- asymmetric models of uniformly rotating cold NSs and diation from the pulsar wind nebulae and produce a consider different energy loss mechanisms. The latter is UV-optical transient that is much brighter than the tra- typically dominated by dipole spin-down which in turn ditional radioactive-decay powered kilonova/macronova depends on the magnetic field strength on the merger (Kasen & Bildsten 2010; Yu et al. 2013; Metzger & Piro remnant's surface. 2014). The non-thermal radiation in the nebula may Our calculation proceeds as follows. We consider first escape the ejecta when the bound-free optical depth is the observed sample of 11 Galactic binary systems with less than unity, generating X-ray emission at the level of well determined individual gravitational masses (listed the spin-down luminosity lasting for a spin-down time or in table1) to simulate (gravitational) chirp masses and until the SMNS collapses into a black hole (BH) (Zhang secondary to primary mass ratios: Mch; q ≡ M2=M1 ≤ 1 5 2013; Metzger & Piro 2014; Murase et al. 2018). It has according to observations . The observed sample can be also been proposed that the long-lasting plateau seen in the X-ray lightcurve of some short gamma-ray bursts 3 http://www.gravity.phys.uwm.edu/rns/ (GRBs) are produced by the nebula emission (Rowlinson 4 To avoid models that are unstable against radial perturbations we only use models that have a central energy density below a critical value defined by the requirement that it leads to the 1 It has been argued that electromagnetic observations of maximum mass of a non-rotating star (Stergioulas & Friedman GW170817 place an upper limit of Mmax . 2:3M (Margalit 1995). & Metzger 2017; Granot et al. 2017; Shibata et al. 2017a; Rez- 5 The cosmological distribution of BNS masses may be different zolla et al. 2018). between the Galactic population and the cosmological one. For 2 For instance, the PSR J1946+2052 BNS system (total mass 2:5± example, it is reasonable to expect that the NS masses could be 0:04M , Stovall et al. 2018), are expected to produce a SMNS, sensitive to the metallicity, which for cosmological events would after accounting for mass loss of a few percent M or more due span a wider range than for the relatively young Galactic BNS to baryonic kilonova ejecta (e.g., Radice et al. 2018b) and then ∼ population. However, this will tend to widen the distribution 0:1M from GW and neutrinos (Hotokezaka et al. 2013; Bernuzzi of survival times and therefore strengthen our results discussed et al. 2016). below.