Hierarchical analysis of gravitational-wave measurements of binary black hole spin-orbit misalignments

Binary black holes may form both through isolated binary evolution and through dynamical interactions in dense stellar environments. The formation channel leaves an imprint on the alignment between the black hole spins and the orbital angular momentum. Gravitational waves from these systems directly encode information about the spin--orbit misalignment angles, allowing them to be (weakly) constrained. Identifying sub-populations of spinning binary black holes will inform us about compact binary formation and evolution. We simulate a mixed population of binary black holes with spin--orbit misalignments modelled under a range of assumptions. We then develop a hierarchical analysis and apply it to mock gravitational-wave observations of these populations. Assuming a population with dimensionless spin magnitudes of $\chi = 0.7$, we show that tens of observations will make it possible to distinguish the presence of subpopulations of coalescing binary black holes based on their spin orientations. With $100$ observations it will be possible to infer the relative fraction of coalescing binary black holes with isotropic spin directions (corresponding to dynamical formation in our models) with a fractional uncertainty of $\sim 40\%$. Meanwhile, only $\sim 5$ observations are sufficient to distinguish between extreme models---all binary black holes either having exactly aligned spins or isotropic spin directions.Comment: 12 pages, 9 figures. Updated to match version published in MNRAS as 10.1093/mnras/stx176

Gravitational wave energy spectrum of a parabolic encounter

We derive an analytic expression for the energy spectrum of gravitational waves from a parabolic Keplerian binary by taking the limit of the Peters and Matthews spectrum for eccentric orbits. This demonstrates that the location of the peak of the energy spectrum depends primarily on the orbital periapse rather than the eccentricity. We compare this weak-field result to strong-field calculations and find it is reasonably accurate (~10%) provided that the azimuthal and radial orbital frequencies do not differ by more than ~10%. For equatorial orbits in the Kerr spacetime, this corresponds to periapse radii of rp > 20M. These results can be used to model radiation bursts from compact objects on highly eccentric orbits about massive black holes in the local Universe, which could be detected by LISA.Comment: 5 pages, 3 figures. Minor changes to match published version; figure 1 corrected; references adde

Parameter estimation on compact binary coalescences with abruptly terminating gravitational waveforms

Gravitational-wave astronomy seeks to extract information about astrophysical systems from the gravitational-wave signals they emit. For coalescing compact-binary sources this requires accurate model templates for the inspiral and, potentially, the subsequent merger and ringdown. Models with frequency-domain waveforms that terminate abruptly in the sensitive band of the detector are often used for parameter-estimation studies. We show that the abrupt waveform termination contains significant information that affects parameter-estimation accuracy. If the sharp cutoff is not physically motivated, this extra information can lead to misleadingly good accuracy claims. We also show that using waveforms with a cutoff as templates to recover complete signals can lead to biases in parameter estimates. We evaluate when the information content in the cutoff is likely to be important in both cases. We also point out that the standard Fisher matrix formalism, frequently employed for approximately predicting parameter-estimation accuracy, cannot properly incorporate an abrupt cutoff that is present in both signals and templates; this observation explains some previously unexpected results found in the literature. These effects emphasize the importance of using complete waveforms with accurate merger and ringdown phases for parameter estimation.Comment: Very minor changes to match published versio

Forward Modeling of Double Neutron Stars: Insights from Highly-Offset Short Gamma-Ray Bursts

We present a detailed analysis of two well-localized, highly offset short gamma-ray bursts---GRB~070809 and GRB~090515---investigating the kinematic evolution of their progenitors from compact object formation until merger. Calibrating to observations of their most probable host galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of hundreds of kilometers per second. Marginalizing over uncertainties in the stellar mass--halo mass relation, we find that the second-born neutron star in the GRB~070809 and GRB~090515 progenitor systems received a natal kick of $\gtrsim 200~\mathrm{km\,s}^{-1}$ at the 78\% and 91\% credible levels, respectively. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.Comment: 18 pages, 7 figures, 1 table. ApJ, in pres

Inference on gravitational waves from coalescences of stellar-mass compact objects and intermediate-mass black holes

Gravitational waves from coalescences of neutron stars or stellar-mass black holes into intermediate-mass black holes (IMBHs) of $\gtrsim 100$ solar masses represent one of the exciting possible sources for advanced gravitational-wave detectors. These sources can provide definitive evidence for the existence of IMBHs, probe globular-cluster dynamics, and potentially serve as tests of general relativity. We analyse the accuracy with which we can measure the masses and spins of the IMBH and its companion in intermediate-mass ratio coalescences. We find that we can identify an IMBH with a mass above $100 ~ M_\odot$ with $95\%$ confidence provided the massive body exceeds $130 ~ M_\odot$. For source masses above $\sim200 ~ M_\odot$, the best measured parameter is the frequency of the quasi-normal ringdown. Consequently, the total mass is measured better than the chirp mass for massive binaries, but the total mass is still partly degenerate with spin, which cannot be accurately measured. Low-frequency detector sensitivity is particularly important for massive sources, since sensitivity to the inspiral phase is critical for measuring the mass of the stellar-mass companion. We show that we can accurately infer source parameters for cosmologically redshifted signals by applying appropriate corrections. We investigate the impact of uncertainty in the model gravitational waveforms and conclude that our main results are likely robust to systematics.Comment: 9 pages, 11 figure

Importance of transient resonances in extreme-mass-ratio inspirals

The inspiral of stellar-mass compact objects, like neutron stars or stellar-mass black holes, into supermassive black holes provides a wealth of information about the strong gravitational-field regime via the emission of gravitational waves. In order to detect and analyse these signals, accurate waveform templates which include the effects of the compact object's gravitational self-force are required. For computational efficiency, adiabatic templates are often used. These accurately reproduce orbit-averaged trajectories arising from the first-order self-force, but neglect other effects, such as transient resonances, where the radial and poloidal fundamental frequencies become commensurate. During such resonances the flux of gravitational waves can be diminished or enhanced, leading to a shift in the compact object's trajectory and the phase of the waveform. We present an evolution scheme for studying the effects of transient resonances and apply this to an astrophysically motivated population. We find that a large proportion of systems encounter a low-order resonance in the later stages of inspiral; however, the resulting effect on signal-to-noise recovery is small as a consequence of the low eccentricity of the inspirals. Neglecting the effects of transient resonances leads to a loss of 4% of detectable signals

Accuracy of inference on the physics of binary evolution from gravitational-wave observations

The properties of the population of merging binary black holes encode some of the uncertain physics of the evolution of massive stars in binaries. The binary black hole merger rate and chirp mass distribution are being measured by ground-based gravitational-wave detectors. We consider isolated binary evolution and explore how accurately the physical model can be constrained with such observations by applying the Fisher information matrix to the merging black hole population simulated with the rapid binary population synthesis code COMPAS. We investigate variations in four COMPAS parameters: common envelope efficiency, kick velocity dispersion, and mass loss rates during the luminous blue variable and Wolf--Rayet stellar evolutionary phases. We find that 1000 observations would constrain these model parameters to a fractional accuracy of a few percent. Given the empirically determined binary black hole merger rate, we can expect gravitational-wave observations alone to place strong constraints on the physics of stellar and binary evolution within a few years.Comment: 12 pages, 9 figures; version accepted by Monthly Notices of the Royal Astronomical Societ

Testing general relativity using golden black-hole binaries

The coalescences of stellar-mass black-hole binaries through their inspiral, merger, and ringdown are among the most promising sources for ground-based gravitational-wave (GW) detectors. If a GW signal is observed with sufficient signal-to-noise ratio, the masses and spins of the black holes can be estimated from just the inspiral part of the signal. Using these estimates of the initial parameters of the binary, the mass and spin of the final black hole can be uniquely predicted making use of general-relativistic numerical simulations. In addition, the mass and spin of the final black hole can be independently estimated from the merger--ringdown part of the signal. If the binary black hole dynamics is correctly described by general relativity (GR), these independent estimates have to be consistent with each other. We present a Bayesian implementation of such a test of general relativity, which allows us to combine the constraints from multiple observations. Using kludge modified GR waveforms, we demonstrate that this test can detect sufficiently large deviations from GR, and outline the expected constraints from upcoming GW observations using the second-generation of ground-based GW detectors.Comment: 5 pages, 2 fig