Title: Distinguishing binary black hole precessional morphologies with gravitational wave observations
The precessional motion of binary black holes can be classified into one of three morphologies, based on the evolution of the angle between the components of the spins in the orbital plane: Circulating, librating around 0, and librating around π. These different morphologies can be related to the binary’s formation channel and are imprinted in the binary’s gravitational wave signal. In this paper, we develop a Bayesian model selection method to determine the preferred spin morphology of a detected binary black hole. The method involves a fast calculation of the morphology which allows us to restrict to a specific morphology in the Bayesian stochastic sampling. We investigate the prospects for distinguishing between the different morphologies using gravitational waves in the Advanced LIGO/Advanced Virgo network with their plus-era sensitivities. For this, we consider fiducial high- and low-mass binaries having different spin magnitudes and signal-to-noise ratios (SNRs). We find that in the cases with high spin and high SNR, the true morphology is strongly favored with log10 Bayes factors ≳4 compared to both alternative morphologies when the binary’s parameters are not close to the boundary between morphologies. However, when the binary parameters are close to the boundary between morphologies, only one alternative morphology is strongly disfavored. In the low-spin, high-SNR cases, the true morphology is still favored with a log10 Bayes factor ∼2 compared to one alternative morphology, while in the low-SNR cases the log10 Bayes factors are at most ∼1 for many binaries. We also consider the gravitational wave signal from GW200129_065458 that has some evidence for precession (modulo data quality issues) and find that there is no preference for a specific morphology. Our method for restricting the prior to a given morphology is publicly available through an easy-to-use python package called bbh_spin_morphology_prior.  more » « less
Award ID(s):
2205920 2308887
PAR ID:
10496766
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review D
Volume:
108
Issue:
10
ISSN:
2470-0010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The spin orientations of spinning binary black hole (BBH) mergers detected by ground-based gravitational wave detectors such as LIGO and Virgo can provide important clues about the formation of such binaries. However, these spin tilts, i.e., the angles between the spin vector of each black hole and the binary’s orbital angular momentum vector, can change due to precessional effects as the black holes evolve from a large separation to their merger. The tilts inferred at a frequency in the sensitive band of the detectors by comparing the signal with theoretical waveforms can thus be significantly different from the tilts when the binary originally formed. These tilts at the binary’s formation are well approximated in many scenarios by evolving the BBH backward in time to a formally infinite separation. Using the tilts at infinite separation also places all binaries on an equal footing in analyzing their population properties. In this paper, we perform parameter estimation for simulated BBHs and investigate the differences between the tilts one infers directly close to merger and those obtained by evolving back to infinite separation. We select simulated observations such that their configurations show particularly large differences in their orientations close to merger and at infinity. While these differences may be buried in the statistical noise for current detections, we show that in future plus-era (A+ and Virgo+) detectors, they can be easily distinguished in some cases. We also consider the tilts at infinity for BBHs in various spin morphologies and at the endpoint of the up-down instability. In particular, we find that we are able to easily identify the up-down instability cases as such from the tilts at infinity. 
    more » « less
  2. Abstract We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 10 14 , and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding p = 10 −3 (≈3 σ ) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields p = 5 × 10 −5 to 1.9 × 10 −4 (≈3.5 σ –4 σ ). Assuming a fiducial f −2/3 characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is 2.4 − 0.6 + 0.7 × 10 − 15 (median + 90% credible interval) at a reference frequency of 1 yr −1 . The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal. 
    more » « less
  3. The most general bound binary black hole system has an eccentric orbit and precessing spins. The detection of such a system with significant eccentricity close to the merger would be a clear signature of dynamical formation. In order to study such systems, it is important to be able to evolve their spins and eccentricity from the larger separations at which the binary formed to the smaller separations at which it is detected, or vice versa. Knowledge of the precessional evolution of the binary’s orbital angular momentum can also be used to twist up aligned-spin eccentric waveform models to create a spin-precessing eccentric waveform model. In this paper, we present a new publicly available code to evolve eccentric, precessing binary black holes using orbit-averaged post-Newtonian (PN) equations from the literature. The spin-precession dynamics is 2PN accurate, i.e., with the leading spin-orbit and spin-spin corrections. The evolution of orbital parameters (orbital frequency, eccentricity, and periastron precession), which follow the quasi-Keplerian parametrization, is 3PN accurate in the point particle terms and includes the leading order spin-orbit and spin-spin effects. All the spin-spin terms include the quadrupole-monopole interaction. The eccentricity enhancement functions in the fluxes use the high-accuracy hyperasymptotic expansions from Loutrel and Yunes []. We discuss various features of the code and study the evolution of the orbital and spin-precession parameters of eccentric, precessing binary black holes. In particular, we study the dependence of the spin morphologies on eccentricity, where we find that the transition point from one spin morphology to another can depend nonmonotonically on eccentricity, and the fraction of binaries in a given morphology at a given point in the evolution of a population depends on the instantaneous eccentricity. 
    more » « less
  4. Gravitational waves from black hole binary mergers carry information about the component spins, but inference is sensitive to analysis assumptions, which may be broken by terrestrial noise transients known as glitches. Using a variety of simulated glitches and gravitational wave signals, we study the conditions under which glitches can bias spin measurements. We confirm the theoretical expectation that inference and subtraction of glitches invariably leaves behind residual power due to statistical uncertainty, no matter the strength (signal-to-noise ratio, SNR) of the original glitch. Next we show that low-SNR glitches — including those below the threshold for flagging data-quality issues — can still significantly bias spin inference. Such biases occur for a range of glitch morphologies, even in cases where glitches and signals are not precisely aligned in phase. Furthermore, we find that residuals of glitch subtraction can result in biases as well. Our results suggest that joint inference of the glitch and gravitational wave parameters, with appropriate models and priors, is required to address these uncertainties inherent in glitch mitigation via subtraction. 
    more » « less
  5. We determine the ability of Cosmic Explorer, a proposed third-generation gravitational-wave observatory, to detect eccentric binary neutron stars and to measure their eccentricity. We find that for a matched-filter search, template banks constructed using binaries in quasicircular orbits are effectual for eccentric neutron star binaries with e<0.004 (e<0.003)is the binary’s eccentricity at a gravitational-wave frequency of 7 Hz. We show that stochastic template placement can be used to construct a matched-filter search for binaries with larger eccentricities and construct an effectual template bank for binaries with e<0.05. We show that the computational cost of both the search for binaries in quasicircular orbits and eccentric orbits is not significantly larger for Cosmic Explorer than for Advanced LIGO and is accessible with present-day computational resources. We investigate Cosmic Explorer’s ability to distinguish between circular and eccentric binaries. We estimate that for a binary with a signal-to-noise ratio of 20 (800), Cosmic Explorer can distinguish between a circular binary and a binary with eccentricity e>~1e-2 (1e-3) at 90% confidence. 
    more » « less