No Arabic abstract
Extending previous work by a number of authors, we have recently presented a new approach in which the detection of gravitational waves from merging neutron star binaries can be used to determine the equation of state of matter at nuclear density and hence the structure of neutron stars. In particular, after performing a large number of numerical-relativity simulations of binaries with nuclear equations of state, we have found that the post-merger emission is characterized by two distinct and robust spectral features. While the high-frequency peak was already shown to be associated with the oscillations of the hypermassive neutron star produced by the merger and to depend on the equation of state, we have highlighted that the low-frequency peak is related to the merger process and to the total compactness of the stars in the binary. This relation is essentially universal and provides a powerful tool to set tight constraints on the equation of state. We here provide additional information on the extensive analysis performed, illustrating the methods used, the tests considered, as well as the robustness of the results. We also discuss additional relations that can be deduced when exploring the data and how these correlate with various properties of the binary. Finally, we present a simple mechanical toy model that explains the main spectral features of the post-merger signal and can even reproduce analytically the complex waveforms emitted right after the merger.
A number of works have shown that important information on the equation of state of matter at nuclear density can be extracted from the gravitational waves emitted by merging neutron-star binaries. We present a comprehensive analysis of the gravitational-wave signal emitted during the inspiral, merger and post-merger of 56 neutron-star binaries. This sample of binaries, arguably the largest studied to date with realistic equations of state, spans across six different nuclear-physics equations of state and ten masses, allowing us to sharpen a number of results recently obtained on the spectral properties of the gravitational-wave signal. Overall we find that: (i) for binaries with masses differing no more than $20%$, the frequency at gravitational-wave amplitudes maximum is related quasi-universally with the tidal deformability of the two stars; (ii) the spectral properties vary during the post-merger phase, with a transient phase lasting a few millisecond after the merger and followed by a quasi-stationary phase; (iii) when distinguishing the spectral peaks between these two phases, a number of ambiguities in the identification of the peaks disappear, leaving a simple and robust picture; (iv) using properly identified frequencies, quasi-universal relations are found between the spectral features and the properties of the neutron stars; (v) for the most salient peaks analytic fitting functions can be obtained in terms of the stellar tidal deformability or compactness. Altogether, these results support the idea that the equation of state of nuclear matter can be constrained tightly when a signal in gravitational waves from binary neutron stars is detected.
We present a robust method to characterize the gravitational wave emission from the remnant of a neutron star coalescence. Our approach makes only minimal assumptions about the morphology of the signal and provides a full posterior probability distribution of the underlying waveform. We apply our method on simulated data from a network of advanced ground-based detectors and demonstrate the gravitational wave signal reconstruction. We study the reconstruction quality for different binary configurations and equations of state for the colliding neutron stars. We show how our method can be used to constrain the yet-uncertain equation of state of neutron star matter. The constraints on the equation of state we derive are complimentary to measurements of the tidal deformation of the colliding neutron stars during the late inspiral phase. In the case of a non-detection of a post-merger signal following a binary neutron star inspiral we show that we can place upper limits on the energy emitted.
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 Virgos 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.
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 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.
Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.