In addition to the emission of gravitational waves (GWs) the coalescence and merger of two neutron stars will produce a variety of electromagnetic (EM) signals. In this work we combine a set of 172 ...numerical relativity simulations performed by different groups employing different numerical codes and we present fits for the mass, kinetic energy, and the velocities of the dynamical ejected material. The obtained fits have residual errors of the same order as the uncertainties of the numerical relativity simulations. Additionally, we comment on the geometry and composition of the ejecta and discuss the influence of the stars' individual rotation. The derived fits can be used to approximate the main properties of kilonovae (macronovae) and radio flares. Furthermore, in cases for which the ejecta mass is known with high precision, we present a way to determine the luminosity and lightcurve of the kilonovae. Overall, the correlation between the binary parameters and the EM signals allows one in the case of a GW detection to approximate possible EM counterparts when first estimates of the masses are available. After a possible kilonovae observation our results could also be used to restrict the region of the parameter space which has to be covered by numerical relativity simulations.
We construct closed-form gravitational waveforms (GWs) with tidal effects for the coalescence of binary neutron stars. The method relies on a new set of eccentricity-reduced and high-resolution ...numerical relativity (NR) simulations and is composed of three steps. First, tidal contributions to the GW phase are extracted from the time-domain NR data. Second, those contributions are employed to fix high-order coefficients in an effective and resummed post-Newtonian expression. Third, frequency-domain tidal approximants are built using the stationary phase approximation. Our tidal approximants are valid from the low frequencies to the strong-field regime. They can be analytically added to any binary black hole GW model to obtain a binary neutron star waveform, either in the time or in the frequency domain. This work provides simple, flexible, and accurate models ready to be used in both searches and parameter estimation of binary neutron star events.
Gravitational wave (GW) astronomy has consolidated its role as a new observational window to reveal the properties of compact binaries in the Universe. In particular, the discovery of the first ...binary neutron star coalescence, GW170817, led to a number of scientific breakthroughs as the possibility to place constraints on the equation of state of cold matter at supranuclear densities. These constraints and all scientific results based on them require accurate models describing the GW signal to extract the source properties from the measured signal. In this article, we study potential systematic biases during the extraction of parameters from non-spinning sources using different descriptions for both, the point-particle dynamics and tidal effects. We find that for the considered cases the mass and spin recovery show almost no systematic bias with respect to the chosen waveform model. However, the extracted tidal effects can be strongly biased, where we find generally that Post-Newtonian approximants predict neutron stars with larger deformability and radii than numerical relativity tuned models. Noteworthy, an increase in the post-Newtonian order in the tidal phasing does not lead to a monotonic change in the estimated properties. We find that for a signal with strength similar to GW170817, but observed with design sensitivity with the network of the two LIGO detectors and the Virgo detector, the estimated tidal parameters can differ by more than a factor of two depending on the employed tidal description of the waveform approximant. This shows the current need for the development of better waveform models to extract reliably the source properties from upcoming GW detections.
In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak ...frequency f2 related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that f(2) is correlated with the tidal coupling constant κ(2)^T that characterizes the binary tidal interactions during the late-inspiral merger. The relation f(2)(κ(2)^T) depends very weakly on the binary total mass, mass ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.
The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of ...these systems. We propose an effective-one-body model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to the merger. Our effective-one-body model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multiorbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semianalytical model that captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.
ABSTRACT The combined detection of a binary neutron star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum – GW170817/AT2017gfo/GRB170817A – ...marks a breakthrough in the field of multimessenger astronomy. Between the plethora of modelling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, and GRB170817 to obtain truly multimessenger constraints on the tidal deformability $\tilde{\Lambda } \in 302,860$, total binary mass M ∈ 2.722, 2.751 M⊙, the radius of a 1.4 solar mass neutron star $R \in 11.3,13.5 \,\,\rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of q ≤ 1.27, all at 90 per cent confidence. Our joint novel analysis uses new phenomenological descriptions of the dynamical ejecta, debris disc mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.
We present new (3+1)-dimensional numerical relativity simulations of the binary neutron star (BNS) mergers that take into account the NS spins. We consider different spin configurations, aligned or ...antialigned to the orbital angular momentum, for equal- and unequal-mass BNSs and for two equations of state. All the simulations employ quasiequilibrium circular initial data in the constant rotational velocity approach, i.e. they are consistent with the Einstein equations and in hydrodynamical equilibrium. We study the NS rotation effect on the energetics, the gravitational waves (GWs) and on the possible electromagnetic (EM) emission associated to dynamical mass ejecta. For dimensionless spin magnitudes of χ∼0.1 we find that both spin-orbit interactions and spin-induced quadrupole deformations affect the late-inspiral merger dynamics. The latter is, however, dominated by finite-size effects. Spin (tidal) effects contribute to GW phase differences up to ∼5 (20) radians accumulated during the last eight orbits to merger. Similarly, after merger the collapse time of the remnant and the GW spectrogram are affected by the NSs rotation. Spin effects in dynamical ejecta are clearly observed in unequal-mass systems in which mass ejection originates from the tidal tail of the companion. Consequently kilonovae and other EM counterparts are affected by spins. We find that spin aligned to the orbital angular momentum leads to brighter EM counterparts than antialigned spin with luminosities up to a factor of 2 higher.
The discovery of the compact binary coalescence in both gravitational waves and electromagnetic radiation marks a breakthrough in the field of multimessenger astronomy and has improved our knowledge ...in a number of research areas. However, an open question is the exact origin of the observables and if one can confirm reliably that GW170817 and its electromagnetic counterparts resulted from a binary neutron star merger. To answer the question if the observation of GW170817, GRB170817A, and AT2017gfo could be explained by the merger of a neutron star with a black hole, we perform a joint multimessenger analysis of the gravitational waves, the short gamma-ray burst, and the kilonova. Assuming a black hole–neutron star system, we derive multimessenger constraints for the tidal deformability of the neutron star of Λ>425 and for the mass ratio of q<2.03 at 90% confidence, with peaks in the likelihood near Λ=830 and q=1.0. Overall, we find that a black hole–neutron star merger could explain the observed signatures; however, our analysis shows that a binary neutron star origin of GW170817 seems more plausible.
We study the gravitational-wave peak luminosity and radiated energy of quasicircular neutron star mergers using a large sample of numerical relativity simulations with different binary parameters and ...input physics. The peak luminosity for all the binaries can be described in terms of the mass ratio and of the leading-order post-Newtonian tidal parameter solely. The mergers resulting in a prompt collapse to black hole have the largest peak luminosities. However, the largest amount of energy per unit mass is radiated by mergers that produce a hypermassive neutron star or a massive neutron star remnant. We quantify the gravitational-wave luminosity of binary neutron star merger events, and set upper limits on the radiated energy and the remnant angular momentum from these events. We find that there is an empirical universal relation connecting the total gravitational radiation and the angular momentum of the remnant. Our results constrain the final spin of the remnant black hole and also indicate that stable neutron star remnant forms with super-Keplerian angular momentum.
Gravitational wave astronomy has established its role in measuring the equation of state governing cold supranuclear matter. To date and in the near future, gravitational wave measurements from ...neutron star binaries are likely to be restricted to the inspiral. However, future upgrades and the next generation of gravitational wave detectors will enable us to detect the gravitational wave signatures emitted after the merger of two stars, at times when densities beyond those in single neutron stars are reached. Therefore, the postmerger gravitational wave signal enables studies of supranuclear matter at its extreme limit. To support this line of research, we present new and updated phenomenological relations between the binary properties and characteristic features of the postmerger evolution. Most notably, we derive an updated relation connecting the mass-weighted tidal deformability and the maximum neutron star mass to the dominant emission frequency of the postmerger spectrum. With the help of a configuration-independent Bayesian analysis using simplified Lorentzian model functions, we find that the main emission frequency of the postmerger remnant, for signal-to-noise ratios of 8 and above, can be extracted within a 1-sigma uncertainty of about 100 Hz for Advanced LIGO and Advanced Virgo's design sensitivities. In some cases, even a postmerger signal-to-noise ratio of 4 can be sufficient to determine the main emission frequency. This will enable us to measure binary and equation-of-state properties from the postmerger, to perform a consistency check between different parts of the binary neutron star coalescence, and to put our physical interpretation of neutron star mergers to the test.