We analyze the properties of the gravitational wave signal emitted after the merger of a binary neutron star system when the remnant survives for more than a 80 ms (and up to 140ms). We employ four different piecewise polytropic equations of state supplemented by an ideal fluid thermal component. We find that the post-merger phase can be subdivided into three phases: an early post-merger phase (where the quadrupole mode and a few subdominant features are active), the intermediate post-merger phase (where only the quadrupole mode is active) and the late post-merger phase (where convective instabilities trigger inertial modes). The inertial modes have frequencies somewhat smaller than the quadrupole modes. In one model, we find an interesting association of a corotation of the quadrupole mode in parts of the star with a revival of its amplitude. The gravitational wave emission of inertial modes in the late post-merger phase is concentrated in a narrow frequency region and is potentially detectable by the planned third-generation detectors. This allows for the possibility of probing not only the cold part of the equation of state, but also its dependence on finite temperature. In view of these results, it will be important to investigate the impact of various type of viscosities on the potential excitation of inertial modes in binary neutron star merger remnants.

We analyse the phenomenological implications of the two-families scenario on the merger of compact stars. That scenario is based on the coexistence of both hadronic stars and strange quark stars. After discussing the classification of the possible mergers, we turn to detailed numerical simulations of the merger of two hadronic stars, i.e., "first family" stars in which delta resonances and hyperons are present, and we show results for the threshold mass of such binaries, for the mass dynamically ejected and the mass of the disk surrounding the post-merger object. We compare these results with those obtained within the one-family scenario and we conclude that relevant signatures of the two-families scenario can be suggested, in particular: the possibility of a rapid collapse to a black hole for masses even smaller than the ones associated to GW170817; during the first milliseconds, oscillations of the postmerger remnant at frequencies higher than the ones obtained in the one-family scenario; a large value of the mass dynamically ejected and a small mass of the disk, for binaries of low total mass. Finally, based on a population synthesis analysis, we present estimates of the number of mergers for: two hadronic stars; hadronic star - strange quark star; two strange quark stars. We show that for unequal mass systems and intermediate values of the total mass, the merger of a hadronic star and a strange quark star is very likely (GW170817 has a possible interpretation into this category of mergers). On the other hand, mergers of two strange quark stars are strongly suppressed.

We present the first very long-term simulations (extending up to ~140 ms after merger) of binary neutron star mergers with piecewise polytropic equations of state and in full general relativity. Our simulations reveal that at a time of 30-50 ms after merger, parts of the star become convectively unstable, which triggers the excitation of inertial modes. The excited inertial modes are sustained up to several tens of milliseconds and are potentially observable by the planned third-generation gravitational-wave detectors at frequencies of a few kilohertz. Since inertial modes depend on the rotation rate of the star and they are triggered by a convective instability in the postmerger remnant, their detection in gravitational waves will provide a unique opportunity to probe the rotational and thermal state of the merger remnant. In addition, our findings have implications for the long-term evolution and stability of binary neutron star remnants

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between $0.75$ and $0.99$. Specifically, these systems are J1756-2251, J0737-3039A, J1906+0746, B1534+12, J0453+1559 and B1913+16. We follow the dynamics of the merger from the late stage of the inspiral process up to $\sim$ 20 ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems ($q=0.75$, J0453+1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapse within the simulated time.

We present results from three-dimensional general relativistic simulations of binary neutron star coalescences and mergers using public codes. We considered equal mass models where the baryon mass of the two Neutron Stars (NS) is $1.4M_{\odot}$, described by four different equations of state (EOS) for the cold nuclear matter (APR4, SLy, H4, and MS1; all parametrized as piecewise polytropes). We started the simulations from four different initial interbinary distances ($40, 44.3, 50$, and $60$ km), including up to the last 16 orbits before merger. That allows to show the effects on the gravitational wave phase evolution, radiated energy and angular momentum due to: the use of different EOSs, the orbital eccentricity present in the initial data and the initial separation (in the simulation) between the two stars. Our results show that eccentricity has a major role in the discrepancy between numerical and analytical waveforms until the very last few orbits, where "tidal" effects and missing high-order post-Newtonian coefficients also play a significant role. We test different methods for extrapolating the gravitational wave signal extracted at finite radii to null infinity. We show that an effective procedure for integrating the Newman-Penrose $\psi_4$ signal to obtain the gravitational wave strain $h$ is to apply a simple high-pass digital filter to $h$ after a time domain integration, where only the two physical motivated integration constants are introduced. That should be preferred to the more common procedures of introducing additional integration constants, integrating in the frequency domain or filtering $\psi_4$ before integration.

We present three-dimensional simulations of the dynamics of binary neutron star (BNS) mergers from the late stage of the inspiral process up to $\sim 20$ ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH). We investigate five equal-mass models of total gravitational mass $2.207$, $2.373$, $2.537$, $2.697$ and $2.854 M_\odot$, respectively, and four unequal mass models with $M_{\mathrm{ADM}}\simeq 2.53\ M_\odot$ and $q\simeq 0.94$, $0.88$, $0.82$, and $0.77$ (where $q = M^{(1)}/M^{(2)}$ is the mass ratio). We use a semi-realistic equation of state (EOS) namely, the seven-segment piece-wise polytropic SLyPP with a thermal component given by $\Gamma_{th} = 1.8$. We have also compared the resulting dynamics (for one model) using both, the BSSN-NOK and CCZ4 methods for the evolution of the gravitational sector, and also different reconstruction methods for the matter sector, namely PPM, WENO and MP5. Our results show agreement and high resolution, but superiority of BSSN-NOK supplemented by WENO reconstruction at lower resolutions. One of the important characteristics of the present investigation is that, for the first time, this has been done using only publicly available open source software, in particular, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. All of the source code and parameters used for the simulations have been made publicly available. This not only makes it possible to re-run and re-analyze our data; it also enables others to directly build upon this work for future research.

We present results about the effect of the use of a stiffer equation of state, namely the ideal-fluid $\Gamma=2.75$ ones, on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change on the critical value of the instability parameter $\beta$ for the emergence of the instability when the adiabatic index $\Gamma$ is changed from 2 to 2.75 in order to mimic the behavior of a realistic equation of state. In particular, we show that the threshold for the onset of the bar-mode instability is reduced by this change in the stiffness and give a precise quantification of the change in value of the critical parameter $\beta_c$. We also extend the analysis to lower values of $\beta$ and show that low-beta shear instabilities are present also in the case of matter described by a simple polytropic equation of state.