No Arabic abstract
We investigate the observability of higher harmonics in gravitational wave signals emitted during the coalescence of binary black holes. We decompose each mode into an overall amplitude, dependent upon the masses and spins of the system, and an orientation-dependent term, dependent upon the inclination and polarization of the source. Using this decomposition, we investigate the significance of higher modes over the parameter space and show that the $ell = 3$, $m = 3$ mode is most significant across much of the sensitive band of ground-based interferometric detectors, with the $ell = 4$, $m = 4$ having a significant contribution at high masses. We introduce the higher mode signal-to-noise ratio (SNR), and show that a simple threshold on this SNR can be used as a criterion for observation of higher harmonics. Finally, we investigate observability in a population of binaries and observe that higher harmonics will only be observable in a few percent of binaries, typically those with unequal masses and viewed close to edge-on.
Current template-based gravitational wave searches for compact binary coalescences (CBC) use waveform models that neglect the higher order modes content of the gravitational radiation emitted, considering only the quadrupolar $(ell,|m|)=(2,2)$ modes. We study the effect of such a neglection for the case of aligned-spin CBC searches for equal-spin (and non-spinning) binary black holes in the context of t
The gravitational waveform of a merging stellar-mass binary is described at leading order by a quadrupolar mode. However, the complete waveform includes higher-order modes, which encode valuable information not accessible from the leading-order mode alone. Despite this, the majority of astrophysical inferences so far obtained with observations of gravitational waves employ only the leading order mode because calculations with higher-order modes are often computationally challenging. We show how to efficiently incorporate higher-order modes into astrophysical inference calculations with a two step procedure. First, we carry out Bayesian parameter estimation using a computationally cheap leading-order-mode waveform, which provides an initial estimate of binary parameters. Second, we weight the initial estimate using higher-order mode waveforms in order to fold in the extra information from the full waveform. We use mock data to demonstrate the effectiveness of this method. We apply the method to each binary black hole event in the first gravitational-wave transient catalog GWTC-1 to obtain posterior distributions and Bayesian evidence with higher-order modes. Performing Bayesian model selection on the events in GWTC-1, we find only a weak preference for waveforms with higher order modes. We discuss how this method can be generalized to a variety of other applications.
A frequentist asymptotic expansion method for error estimation is employed for a network of gravitational wave detectors to assess the amount of information that can be extracted from gravitational wave observations. Mathematically we derive lower bounds in the errors that any parameter estimator will have in the absence of prior knowledge to distinguish between the post-Einsteinian (ppE) description of coalescing binary systems and that of general relativity. When such errors are smaller than the parameter value, there is possibility to detect these violations from GR. A parameter space with inclusion of dominant dephasing ppE parameters $(beta, b)$ is used for a study of first- and second-order (co)variance expansions, focusing on the inspiral stage of a nonspinning binary system of zero eccentricity detectible through Adv. LIGO and Adv. Virgo. Our procedure is an improvement of the Cram{e}r-Rao Lower Bound. When Bayesian errors are lower than our bound it means that they depend critically on the priors. The analysis indicates the possibility of constraining deviations from GR in inspiral SNR ($rho sim 15-17$) regimes that are achievable in upcoming scientific runs (GW150914 had an inspiral SNR $sim 12$). The errors on $beta$ also increase errors of other parameters such as the chirp mass $mathcal{M}$ and symmetric mass ratio $eta$. Application is done to existing alternative theories of gravity, which include modified dispersion relation of the waveform, non-spinning models of quadratic modified gravity, and dipole gravitational radiation (i.e., Brans-Dicke type) modifications.
Gravitational wave observations of compact binaries allow us to test general relativity (and modifications thereof) in the strong and highly-dynamical field regime of gravity. Here we confront two extensions to general relativity, dynamical Chern-Simons and Einstein-dilaton-Gauss-Bonnet theories, against the gravitational wave sources from the GWTC-1 and GWTC-2 catalogs by the LIGO-Virgo Collaboration. By stacking the posterior of individual events, we strengthen the constraint on the square root of the coupling parameter in Einstein-dilaton-Gauss-Bonnet gravity to $sqrt{alpha_{rm tiny EdGB}} < 1.7$ km, but we are unable to place meaningful constraints on dynamical Chern-Simons gravity. Importantly, we also show that our bounds are robust to (i) the choice of general-relativity base waveform model, upon which we add modifications, (ii) unknown higher post-Newtonian order terms in the modifications to general relativity, (iii) the small-coupling approximation, and (iv) uncertainties on the nature of the constituent compact objects.
Estimates of the source parameters of gravitational-wave (GW) events produced by compact binary mergers rely on theoretical models for the GW signal. We present the first frequency-domain model for inspiral, merger and ringdown of the GW signal from precessing binary-black-hole systems that also includes multipoles beyond the leading-order quadrupole. Our model, {tt PhenomPv3HM}, is a combination of the higher-multipole non-precessing model {tt PhenomHM} and the spin-precessing model {tt PhenomPv3} that includes two-spin precession via a dynamical rotation of the GW multipoles. We validate the new model by comparing to a large set of precessing numerical-relativity simulations and find excellent agreement across the majority of the parameter space they cover. For mass ratios $<5$ the mismatch improves, on average, from $sim6%$ to $sim 2%$ compared to {tt PhenomPv3} when we include higher multipoles in the model. However, we find mismatches $sim8%$ for the mass-ratio $6$ and highly spinning simulation. As a first application of the new model we have analysed the binary black hole event GW170729. We find larger values for the primary black hole mass of $58.25^{+11.73}_{-12.53} , M_odot$ (90% credible interval). The lower limit ($sim 46 , M_odot$) is comparable to the proposed maximum black hole mass predicted by different stellar evolution models due to the pulsation pair-instability supernova (PPISN) mechanism. If we assume that the primary ac{BH} in GW170729 formed through a PPISN then out of the four PPISN models we considered only the model of Woosley (2017) is consistent with our mass measurements at the 90% level.