Search for gravitational waves from a long-lived remnant of the binary neutron star merger GW170817

One unanswered question about the binary neutron star coalescence GW170817 is the nature of its post-merger remnant. A previous search for post-merger gravitational waves targeted high-frequency signals from a possible neutron star remnant with a maximum signal duration of 500 s. Here we revisit the neutron star remnant scenario with a focus on longer signal durations up until the end of the Second Advanced LIGO-Virgo Observing run, 8.5 days after the coalescence of GW170817. The main physical scenario for such emission is the power-law spindown of a massive magnetar-like remnant. We use four independent search algorithms with varying degrees of restrictiveness on the signal waveformand different ways of dealing with noise artefacts. In agreement with theoretical estimates, we find no significant signal candidates. Through simulated signals, we quantify that with the current detector sensitivity, nowhere in the studied parameter space are we sensitive to a signal from more than 1 Mpc away, compared to the actual distance of 40 Mpc. This study however serves as a prototype for post-merger analyses in future observing runs with expected higher sensitivity.

Spectral properties of the post-merger gravitational-wave signal from binary neutron stars

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.

Inferring the post-merger gravitational wave emission from binary neutron star coalescences

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.

Computing Fast and Reliable Gravitational Waveforms of Binary Neutron Star Merger Remnants

Gravitational waves have been detected from the inspiral of a binary neutron-star, GW170817, which allowed constraints to be placed on the neutron star equation of state. The equation of state can be further constrained if gravitational waves from a post-merger remnant are detected. Post-merger waveforms are currently generated by numerical-relativity simulations, which are computationally expensive. Here we introduce a hierarchical model trained on numerical-relativity simulations, which can generate reliable post-merger spectra in a fraction of a second. Our spectra have mean fitting factors of 0.95, which compares to fitting factors of 0.76 and 0.85 between different numerical-relativity codes that simulate the same physical system. This method is the first step towards generating large template banks of spectra for use in post-merger detection and parameter estimation.

Reanalysis of the binary neutron star mergers GW170817 and GW190425 using numerical-relativity calibrated waveform models

We reanalyze gravitational waves from binary-neutron-star mergers GW170817 and GW190425 using a numerical-relativity (NR) calibrated waveform model, the TF2+_Kyoto model, which includes nonlinear tidal terms. For GW170817, by imposing a uniform prior on the binary tidal deformability $tilde{Lambda}$, the symmetric $90%$ credible interval of $tilde{Lambda}$ is estimated to be $481^{+436}_{-359}$ and $402^{+465}_{-279}$ for the case of $f_mathrm{max}=1000$ and $2048~mathrm{Hz}$, respectively, where $f_mathrm{max}$ is the maximum frequency in the analysis. We also reanalyze the event with other waveform models: two post-Newtonian waveform models (TF2_PNTidal and TF2+_PNTidal), the TF2+_NRTidal model that is another NR calibrated waveform model, and its upgrade, the TF2+_NRTidalv2 model. While estimates of parameters other than $tilde{Lambda}$ are broadly consistent among various waveform models, our results indicate that estimates of $tilde{Lambda}$ depend on waveform models. However, the difference is smaller than the statistical error. For GW190425, we can only obtain little information on the binary tidal deformability. The systematic difference among the NR calibrated waveform models will become significant to measure $tilde{Lambda}$ as the number of detectors and events increase and sensitivities of detectors are improved.