No Arabic abstract
In August 2017, the first detection of a binary neutron star merger, GW170817, made it possible to study neutron stars in compact binary systems using gravitational waves. Despite being the loudest (in terms of signal-to-noise ratio) gravitational wave detected to date, it was not possible to unequivocally determine that GW170817 was caused by the merger of two neutron stars instead of two black holes from the gravitational-wave data alone. That distinction was largely due to the accompanying electromagnetic counterpart. This raises the question: under what circumstances can gravitational-wave data alone, in the absence of an electromagnetic signal, be used to distinguish between different types of mergers? Here, we study whether a neutron-star--black-hole binary merger can be distinguished from a binary black hole merger using gravitational-wave data alone. We build on earlier results using chiral effective field theory to explore whether the data from LIGO and Virgo, LIGO A+, LIGO Voyager, or Cosmic Explorer could lead to such a distinction. The results suggest that the present LIGO-Virgo detector network will most likely be unable to distinguish between these systems even with the planned near-term upgrades. However, given an event with favorable parameters, third-generation instruments such as Cosmic Explorer will be capable of making this distinction. This result further strengthens the science case for third-generation detectors.
In models of minicharged dark matter associated with a hidden $U(1)$ symmetry, astrophysical black holes may acquire a dark charge, in such a way that the inspiral dynamics of binary black holes can be formally described by an Einstein-Maxwell theory. Charges enter the gravitational wave signal predominantly through a dipole term, but their effect is known to effectively first post-Newtonian order in the phase, which enables measuring the size of the charge-to-mass ratios, $|q_i/m_i|$, $i = 1,2$, of the individual black holes in a binary. We set up a Bayesian analysis to discover, or constrain, dark charges on binary black holes. After testing our framework in simulations, we apply it to selected binary black hole signals from the second Gravitational Wave Transient Catalog (GWTC-2), namely those with low masses so that most of the signal-to-noise ratio is in the inspiral regime. We find no evidence for charges on the black holes, and place typical 1-$sigma$ bounds on the charge-to-mass ratios of $|q_i/m_i| lesssim 0.2 - 0.3$.
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.
As current gravitational wave (GW) detectors increase in sensitivity, and particularly as new instruments are being planned, there is the possibility that ground-based GW detectors will observe GWs from highly eccentric neutron star binaries. We present the first detailed study of highly eccentric BNS systems with full (3+1)D numerical relativity simulations using consistent initial conditions, i.e., setups which are in agreement with the Einstein equations and with the equations of general relativistic hydrodynamics in equilibrium. Overall, our simulations cover two different equations of state (EOSs), two different spin configurations, and three to four different initial eccentricities for each pairing of EOS and spin. We extract from the simulated waveforms the frequency of the f-mode oscillations induced during close encounters before the merger of the two stars. The extracted frequency is in good agreement with f-mode oscillations of individual stars for the irrotational cases, which allows an independent measure of the supranuclear equation of state not accessible for binaries on quasi-circular orbits. The energy stored in these f-mode oscillations can be as large as $10^{-3}M_odot sim 10^{51}$ erg, even with a soft EOS. In order to estimate the stored energy, we also examine the effects of mode mixing due to the stars offset from the origin on the f-mode contribution to the GW signal. While in general (eccentric) neutron star mergers produce bright electromagnetic counterparts, we find that the luminosity decreases when the eccentricity becomes too large, due to a decrease of the ejecta mass. Finally, the use of consistent initial configurations also allows us to produce high-quality waveforms for different eccentricities which can be used as a testbed for waveform model development of highly eccentric binary neutron star systems.
Third-generation (3G) gravitational-wave detectors will observe thousands of coalescing neutron star binaries with unprecedented fidelity. Extracting the highest precision science from these signals is expected to be challenging owing to both high signal-to-noise ratios and long-duration signals. We demonstrate that current Bayesian inference paradigms can be extended to the analysis of binary neutron star signals without breaking the computational bank. We construct reduced order models for $sim 90,mathrm{minute}$ long gravitational-wave signals, covering the observing band ($5-2048,mathrm{Hz}$), speeding up inference by a factor of $sim 1.3times 10^4$ compared to the calculation times without reduced order models. The reduced order models incorporate key physics including the effects of tidal deformability, amplitude modulation due to the Earths rotation, and spin-induced orbital precession. We show how reduced order modeling can accelerate inference on data containing multiple, overlapping gravitational-wave signals, and determine the speedup as a function of the number of overlapping signals. Thus, we conclude that Bayesian inference is computationally tractable for the long-lived, overlapping, high signal-to-noise-ratio events present in 3G observatories.
We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession (wobbling) of the orbital plane, while three show a bobbing motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of $sim 40 rm km/s$ for the studied systems. Furthermore, for the chosen configurations, antialigned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $sim$ 16 gravitational wave cycles.