The anticipated enhancements in detector sensitivity and the corresponding increase in the number of gravitational wave detections will make it possible to estimate parameters of compact binaries with greater accuracy assuming general relativity(GR),
and also to carry out sharper tests of GR itself. Crucial to these procedures are accurate gravitational waveform models. The systematic errors of the models must stay below statistical errors to prevent biases in parameter estimation and to carry out meaningful tests of GR. Comparisons of the models against numerical relativity (NR) waveforms provide an excellent measure of systematic errors. A complementary approach is to use balance laws provided by Einsteins equations to measure faithfulness of a candidate waveform against exact GR. Each balance law focuses on a physical observable and measures the accuracy of the candidate waveform vis a vis that observable. Therefore, this analysis can provide new physical insights into sources of errors. In this paper we focus on the angular momentum balance law, using post-Newtonian theory to calculate the initial angular momentum, surrogate fits to obtain the remnant spin and waveforms from models to calculate the flux. The consistency check provided by the angular momentum balance law brings out the marked improvement in the passage from texttt{IMRPhenomPv2} to texttt{IMRPhenomXPHM} and from texttt{SEOBNRv3} to texttt{SEOBNRv4PHM} and shows that the most rece
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 wa
ve 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.
Observations of gravitational waves from compact binary mergers have enabled unique tests of general relativity in the dynamical and non-linear regimes. One of the most important such tests are constraints on the post-Newtonian (PN) corrections to th
e phase of the gravitational wave signal. The values of these PN coefficients can be calculated within standard general relativity, and these values are different in many alternate theories of gravity. It is clearly of great interest to constrain these deviations based on gravitational wave observations. In the majority of such tests which have been carried out, and which yield by far the most stringent constraints, it is common to vary these PN coefficients individually. While this might in principle be useful for detecting certain deviations from standard general relativity, it is a serious limitation. For example, we would expect alternate theories of gravity to generically have additional parameters. The corrections to the PN coefficients would be expected to depend on these additional non-GR parameters whence, we expect that the various PN coefficients to be highly correlated. We present an alternate analysis here using data from the binary neutron star coalescence GW170817. Our analysis uses an appropriate linear combination of non-GR parameters that represent absolute deviations from the corresponding post-Newtonian inspiral coefficients in the TaylorF2 approximant phase. These combinations represent uncorrelated non-GR parameters which correspond to principal directions of their covariance matrix in the parameter subspace. Our results illustrate good agreement with GR. In particular, the integral non-GR phase is $Psi_{mbox{non-GR}} = (0.447pm253)times10^{-1}$ and the deviation from GR percentile is $p^{mbox{Dev-GR}}_{n}=25.85%$.
We present the results of a search in LIGO O2 public data for continuous gravitational waves from the neutron star in the low-mass X-ray binary Scorpius X-1. We search for signals with $approx$ constant frequency in the range 40-180 Hz. Thanks to the
efficiency of our search pipeline we can use a long coherence time and achieve unprecedented sensitivity, significantly improving on existing results. This is the first search that has been able to probe gravitational wave amplitudes that could balance the accretion torque at the neutron star radius. Our search excludes emission at this level between 67.5 Hz and 131.5 Hz, for an inclination angle $44^circ pm 6^circ$ derived from radio observations (Fomalont et al. 2001), and assuming that the spin axis is perpendicular to the orbital plane. If the torque arm is $approx $ 26 km -- a conservative estimate of the alfven radius -- our results are more constraining than the indirect limit across the band. This allows us to exclude certain mass-radius combinations and to place upper limits on the strength of the stars magnetic field. We also correct a mistake that appears in the literature in the equation that gives the gravitational wave amplitude at the torque balance (Abbott et al. 2017b, 2019a) and we re-interpret the associated latest LIGO/Virgo results in light of this.
Full, non-linear general relativity predicts a memory effect for gravitational waves. For compact binary coalescence, the total gravitational memory serves as an inferred observable, conceptually on the same footing as the mass and the spin of the fi
nal black hole. Given candidate waveforms for any LIGO event, then, one can calculate the posterior probability distribution functions for the total gravitational memory, and use them to compare and contrast the waveforms. In this paper we present these posterior distributions for the binary black hole merger events reported in the first Gravitational Wave Transient Catalog (GWTC-1), using the Phenomenological and Effective-One-Body waveforms. On the whole, the two sets of posterior distributions agree with each other quite well though we find larger discrepancies for the $ell=2, m=1$ mode of the memory. This signals a possible source of systematic errors that was not captured by the posterior distributions of other inferred observables. Thus, the posterior distributions of various angular modes of total memory can serve as diagnostic tools to further improve the waveforms. Analyses such as this would be valuable especially for future events as the sensitivity of ground based detectors improves, and for LISA which could measure the total gravitational memory directly.
We have shown previously that a merger of marginally outer trapped surfaces (MOTSs) occurs in a binary black hole merger and that there is a continuous sequence of MOTSs which connects the initial two black holes to the final one. In this paper, we c
onfirm this scenario numerically and we detail further improvements in the numerical methods for locating MOTSs. With these improvements, we confirm the merger scenario and demonstrate the existence of self-intersecting MOTSs formed in the immediate aftermath of the merger. These results will allow us to track physical quantities across the non-linear merger process and to potentially infer properties of the merger from gravitational wave observations.
We find strong numerical evidence for a new phenomenon in a binary black hole spacetime, namely the merger of marginally outer trapped surfaces (MOTSs). By simulating the head-on collision of two non-spinning unequal mass black holes, we observe that
the MOTS associated with the final black hole merges with the two initially disjoint surfaces corresponding to the two initial black holes. This yields a connected sequence of MOTSs interpolating between the initial and final state all the way through the non-linear binary black hole merger process. In addition, we show the existence of a MOTS with self-intersections formed immediately after the merger. This scenario now allows us to track physical quantities (such as mass, angular momentum, higher multipoles, and fluxes) across the merger, which can be potentially compared with the gravitational wave signal in the wave-zone, and with observations by gravitational wave detectors. This also suggests a possibility of proving the Penrose inequality mathematically for generic astrophysical binary back hole configurations.
Marginally outer trapped surfaces (MOTSs, or marginal surfaces in short) are routinely used in numerical simulations of black hole spacetimes. They are an invaluable tool for locating and characterizing black holes quasi-locally in real time while th
e simulation is ongoing. It is often believed that a MOTS can behave unpredictably under time evolution; an existing MOTS can disappear, and a new one can appear without any apparent reason. In this paper we show that in fact the behavior of a MOTS is perfectly predictable and its behavior is dictated by a single real parameter, the emph{stability parameter}, which can be monitored during the course of a numerical simulation. We demonstrate the utility of the stability parameter to fully understand the variety of marginal surfaces that can be present in binary black hole initial data. We also develop a new horizon finder capable of locating very highly distorted marginal surfaces and we show that even in these cases, the stability parameter perfectly predicts the existence and stability of marginal surfaces.
In this paper we design a search for continuous gravitational waves from three supernova remnants: Vela Jr., Cassiopeia A (Cas A) and G347.3. These systems might harbor rapidly rotating neutron stars emitting quasi-periodic gravitational radiation de
tectable by the advanced LIGO detectors. Our search is designed to use the volunteer computing project Einstein@Home for a few months and assumes the sensitivity and duty cycles of the advanced LIGO detectors during their first science run. For all three supernova remnants, the sky-positions of their central compact objects are well known but the frequency and spin-down rates of the neutron stars are unknown which makes the searches computationally limited. In a previous paper we have proposed a general framework for deciding on what target we should spend computational resources and in what proportion, what frequency and spin-down ranges we should search for every target, and with what search set-up. Here we further expand this framework and apply it to design a search directed at detecting continuous gravitational wave signals from the most promising three supernova remnants identified as such in the previous work. Our optimization procedure yields broad frequency and spin-down searches for all three objects, at an unprecedented level of sensitivity: The smallest detectable gravitational wave strain $h_0$ for Cas A is expected to be 2 times smaller than the most sensitive upper-limits published to date, and our proposed search, which was set-up and ran on the volunteer computing project Einstein@Home, covers a much larger frequency range.
Matched filtering is a commonly used technique in gravitational wave searches for signals from compact binary systems and from rapidly rotating neutron stars. A common issue in these searches is dealing with four extrinsic parameters which do not aff
ect the phase evolution of the system: the overall amplitude, initial phase, and two angles determining the overall orientation of the system. The F-statistic maximizes the likelihood function analytically over these parameters, while the B-statistic marginalizes over them. The B-statistic, while potentially more powerful and capable of incorporating astrophysical priors, is not as widely used because of the computational difficulty of performing the marginalization. In this paper we address this difficulty and show how the marginalization can be done analytically by combining the four parameters into a set of complex amplitudes. The results of this paper are applicable to both transient non-precessing binary coalescence events, and to long lived signals from rapidly rotating neutron stars.
