Gravitational wave observations of compact binary coalescences provide precision probes of strong-field gravity. There is thus now a standard set of null tests of general relativity (GR) applied to LIGO-Virgo detections and many more such tests propo
sed. However, the relation between all these tests is not yet well understood. We start to investigate this by applying a set of standard tests to simulated observations of binary black holes in GR and with phenomenological deviations from GR. The phenomenological deviations include self-consistent modifications to the energy flux in an effective-one-body (EOB) model, the deviations used in the second post-Newtonian (2PN) TIGER and FTA parameterized tests, and the dispersive propagation due to a massive graviton. We consider four types of tests: residuals, inspiral-merger-ringdown consistency, parameterized (TIGER and FTA), and modified dispersion relation. We also check the consistency of the unmodeled reconstruction of the waveforms with the waveform recovered using GR templates. These tests are applied to simulated observations similar to GW150914 with both large and small deviations from GR and similar to GW170608 just with small deviations from GR. We find that while very large deviations from GR are picked up with high significance by almost all tests, more moderate deviations are picked up by only a few tests, and some deviations are not recognized as GR violations by any test at the moderate signal-to-noise ratios we consider. Moreover, the tests that identify various deviations with high significance are not necessarily the expected ones. We also find that the 2PN (1PN) TIGER and FTA tests recover much smaller deviations than the true values in the modified EOB (massive graviton) case. Additionally, we find that of the GR deviations we consider, the residuals test is only able to detect extreme deviations from GR. (Abridged)
Two important parameters inferred from the gravitational wave signals of binaries of precessing black holes are the spin tilt angles, i.e., the angles at which the black holes spin axes are inclined with respect to the binarys orbital angular momentu
m. The LIGO-Virgo parameter estimation analyses currently provide spin tilts at a fiducial reference frequency, often the lowest frequency used in the data analysis. However, the most astrophysically interesting quantities are the spin tilts when the binary was formed, which can be significantly different from those at the reference frequency for strongly precessing binaries. The spin tilts at formally infinite separation are a good approximation to the tilts at formation in many formation channels and can be computed efficiently for binary black holes using precession-averaged evolution. Here, we present a new code for computing the tilts at infinity that combines the precession-averaged evolution with orbit-averaged evolution at high frequencies and illustrate its application to GW190521 and other binary black hole detections from O3a. We have empirically determined the transition frequency between the orbit-averaged and precession-averaged evolution to produce tilts at infinity with a given accuracy. We also have regularized the precession-averaged equations in order to obtain good accuracy for the very close-to-equal-mass binary parameters encountered in practice. This additionally allows us to investigate the singular equal-mass limit of the precession-averaged expressions, where we find an approximate scaling of $1/(1 - q)$ with the mass ratio $q$.
LIGO and Virgo have recently observed a number of gravitational wave (GW) signals that are fully consistent with being emitted by binary black holes described by general relativity. However, there are theoretical proposals of exotic objects that can
be massive and compact enough to be easily confused with black holes. Nevertheless, these objects differ from black holes in having nonzero tidal deformabilities, which can allow one to distinguish binaries containing such objects from binary black holes using GW observations. Using full Bayesian parameter estimation, we investigate the possibility of constraining the parameter space of such black hole mimickers with upcoming GW observations. Employing perfect fluid stars with a polytropic equation of state as a simple model that can encompass a variety of possible black hole mimickers, we show how the observed masses and tidal deformabilities of a binary constrain the equation of state. We also show how such constraints can be used to rule out some simple models of boson stars.
In the adiabatic post-Newtonian (PN) approximation, the phase evolution of gravitational waves (GWs) from inspiralling compact binaries in quasicircular orbits is computed by equating the change in binding energy with the GW flux. This energy balance
equation can be solved in different ways, which result in multiple approximants of the PN waveforms. Due to the poor convergence of the PN expansion, these approximants tend to differ from each other during the late inspiral. Which of these approximants should be chosen as templates for detection and parameter estimation of GWs from inspiraling compact binaries is not obvious. In this paper, we present estimates of the effective higher order (beyond the currently available 4PN and 3.5PN) non-spinning terms in the PN expansion of the binding energy and the GW flux that minimize the difference of multiple PN approximants (TaylorT1, TaylorT2, TaylorT4, TaylorF2) with effective one body waveforms calibrated to numerical relativity (EOBNR). We show that PN approximants constructed using the effective higher order terms show significantly better agreement (as compared to 3.5PN) with the inspiral part of the EOBNR. For non-spinning binaries with component masses $m_{1,2} in [1.4 M_odot, 15 M_odot]$, most of the approximants have a match (faithfulness) of better than 99% with both EOBNR and each other.
Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spi
ns, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $sim10^{-2}M_odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]
