Inferring spin tilts at formation from gravitational wave observations of binary black holes: Interfacing precession-averaged and orbit-averaged spin evolution

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 momentum. 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$.

Testing the multipole structure and conservative dynamics of compact binaries using gravitational wave observations: The spinning case

In an earlier work [S. Kastha et al., PRD {bf 98}, 124033 (2018)], we developed the {it parametrized multipolar gravitational wave phasing formula} to test general relativity, for the non-spinning compact binaries in quasi-circular orbit. In this paper, we extend the method and include the important effect of spins in the inspiral dynamics. Furthermore, we consider parametric scaling of PN coefficients of the conserved energy for the compact binary, resulting in the parametrized phasing formula for non-precessing spinning compact binaries in quasi-circular orbit. We also compute the projected accuracies with which the second and third generation ground-based gravitational wave detector networks as well as the planned space-based detector LISA will be able to measure the multipole deformation parameters and the binding energy parameters. Based on different source configurations, we find that a network of third-generation detectors would have comparable ability to that of LISA in constraining the conservative and dissipative dynamics of the compact binary systems. This parametrized multipolar waveform would be extremely useful not only in deriving the first upper limits on any deviations of the multipole and the binding energy coefficients from general relativity using the gravitational wave detections, but also for science case studies of next generation gravitational wave detectors.

Testing the multipole structure of compact binaries using gravitational wave observations

We propose a novel method to test the consistency of the multipole moments of compact binary systems with the predictions of General Relativity (GR). The multipole moments of a compact binary system, known in terms of symmetric and trace-free tensors, are used to calculate the gravitational waveforms from compact binaries within the post-Newtonian (PN) formalism. For nonspinning compact binaries, we derive the gravitational wave phasing formula, in the frequency domain, parametrizing each PN order term in terms of the multipole moments which contribute to that order. Using GW observations, this {it{parametrized multipolar phasing}} would allow us to derive the bounds on possible departures from the multipole structure of GR and hence constrain the parameter space of alternative theories of gravity. We compute the projected accuracies with which the second generation ground-based detectors, such as Advanced Laser Interferometer Gravitational-wave Observatory (LIGO), the third generation detectors such as Einstein Telescope and Cosmic Explorer, as well as space-based detector Laser Interferometer Space Antenna (LISA) will be able to measure these multipole parameters. We find that while Advanced LIGO can measure the first two or three multipole coefficients with good accuracy, Cosmic Explorer and Einstein Telescope may be able to measure the first four multipole coefficients which enter the phasing formula. Intermediate mass ratio inspirals, with mass ratio of several tens, in the frequency band of planned space-based LISA mission should be able to measure all the seven multipole coefficients which appear in the 3.5PN phasing formula. Our finding highlights the importance of this class of sources for probing the strong-field gravity regime. The proposed test will facilitate the first probe of the multipolar structure of Einsteins general relativity.