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 su
pplemented 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 me
rgers, 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
In this work we analyze the gravitational wave signal from hypermassive neutron stars formed after the merger of binary neutron star systems, focusing on its spectral features. The gravitational wave signals are extracted from numerical relativity si
mulations of models already considered by De Pietri et al. [Phys. Rev. D 93, 064047 (2016)], Maione et al. [Classical Quantum Gravity 33, 175009 (2016)], and Feo et al. [Classical Quantum Gravity 34, 034001 (2017)], and allow us to study the effect of the total baryonic mass of such systems (from $2.4 M_{odot}$ to $3 M_{odot}$), the mass ratio (up to $q = 0.77$), and the neutron star equation of state, both in equal and highly unequal mass binaries. We use the peaks we find in the gravitational spectrum as an independent test of already published hypotheses of their physical origin and empirical relations linking them with the characteristics of the merging neutron stars. In particular, we highlight the effects of the mass ratio, which in the past was often neglected. We also analyze the temporal evolution of the emission frequencies. Finally, we introduce a modern variant of Pronys method to analyze the gravitational wave postmerger emission as a sum of complex exponentials, trying to overcome some drawbacks of both Fourier spectra and least-squares fitting. Overall, the spectral properties of the postmerger signal observed in our simulation are in agreement with those proposed by other groups. More specifically, we find that the analysis of Bauswein and Stergioulas [Phys. Rev. D 91, 124056 (2015)] is particularly effective for binaries with very low masses or with a small mass ratio and that the mechanical toy model of Takami et al. [Phys. Rev. D 91, 064001 (2015)] provides a comprehensive and accurate description of the early stages of the postmerger.
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+1
559 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}$, descr
ibed 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 bl
ack 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 $qsimeq 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 on the effect of the stiffness of the equation of state on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change in the threshold for the emerg
ence of the instability for a range of the adiabatic $Gamma$ index from 2.0 to 3.0, including two values chosen to mimic more realistic equations of state at high densities.
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 deter
mine 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.
We present Monte Carlo results for the thermodynamics of pure SU(N) gauge theories with $N=2,...,6$ in 2+1 dimensions. We focus on the confined phase region $T<T_c$ and study thermodynamics variables such as the trace of the energy-momentum tensor, p
ressure, energy and entropy density using the integral method. We also investigate scaling properties with $N$ of the different observables. We compare our results with a gas of free glueballs and the bosonic string predictions for the Hagedorn spectrum.
