Refine
Has Fulltext
- no (4)
Document Type
- Article (4)
Language
- English (4)
Is part of the Bibliography
- yes (4)
Institute
- Institut für Physik und Astronomie (4) (remove)
The third observing run of the LIGO-Virgo Collaboration has resulted in many gravitational wave detections, including the binary neutron star merger GW190425. However, none of these events have been accompanied with an electromagnetic transient found during extensive follow-up searches. In this article, we perform new numerical-relativity simulations of binary neutron star and black hole-neutron star systems that have a chirp mass consistent with GW190425. Assuming that the GW190425's sky location was covered with sufficient accuracy during the electromagnetic follow-up searches, we investigate whether the nondetection of the kilonova is compatible with the source parameters estimated through the gravitational -wave analysis and how one can use this information to place constraints on the properties of the system. Our simulations suggest that GW190425 is incompatible with an unequal mass binary neutron star merger with a mass ratio q < 0.8 when considering stiff or moderately stiff equations of state if the binary was face on and covered by the observation. Our analysis shows that a detailed observational result for kilonovae will be useful to constrain the mass ratio of binary neutron stars in future events.
With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and antialigned spins within a range of dimensionless spins of chi similar to [-0.28, 0.58]. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 arc not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.
Numerical-relativity simulations are essential for studying the last stages of the binary neutron star coalescence. Unfortunately, for stable simulations there is the need to add an artificial low-density atmosphere. Here we discuss a new framework in which we can effectively set the density surrounding the neutron stars to zero to ensure a more accurate simulation. We test our method with a number of single star test cases and for an equal-mass binary neutron star simulation. While the bulk motion of the system is not influenced, and hence there is no improvement with respect to the emitted gravitational-wave signal, we find that the new approach is superior with respect to mass conservation and it allows a much better tracking of outward moving material. This will allow a more accurate simulation of the ejected material and supports the interpretation of present and future multimessenger observations with more accurate numerical-relativity simulations.
The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical relativity simulations for four different physical setups with mass ratios q = 1.25, 1.50, 1.75, 2.00, and total gravitational mass M = 2.7 M???. Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately second-order converging results for the dominant (2,2) mode, but also the subdominant (2,1), (3,3), and (4,4) modes, while generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model used to describe tidal effects is a valid description for high-mass-ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.