Gravitational-wave parameter estimation for compact binary signals typically relies on sequential estimation of the properties of the detector Gaussian noise and of the binary parameters. This ...procedure assumes that the noise variance, expressed through its power spectral density, is perfectly known in advance. We assess the impact of this approximation on the estimated parameters by means of an analysis that simultaneously estimates the noise and compact binary parameters, thus allowing us to marginalize over uncertainty in the noise properties. We compare the traditional sequential estimation method and the new full marginalization method using events from the GWTC-3 catalog. We find that the recovered signals and inferred parameters agree to within their statistical measurement uncertainty. At current detector sensitivities, uncertainty about the noise power spectral density is a subdominant effect compared to other sources of uncertainty.
An isotropic stochastic background of nanohertz gravitational waves creates excess residual power in pulsar-timing-array datasets, with characteristic inter-pulsar correlations described by the ...Hellings-Downs function. These correlations appear as nondiagonal terms in the noise covariance matrix, which must be inverted to obtain the pulsar-timing-array likelihood. Searches for the stochastic background, which require many likelihood evaluations, are therefore quite computationally expensive. We propose a more efficient method: we first compute approximate posteriors by ignoring cross correlations, and then reweight them to exact posteriors via importance sampling. We show that this technique results in accurate posteriors and marginal likelihood ratios, because the approximate and exact posteriors are similar, which makes reweighting especially accurate. The Bayes ratio between the marginal likelihoods of the exact and approximate models, commonly used as a detection statistic, is also estimated reliably by our method, up to ratios of at least \(10^6\).
Measurement of spin-precession in black hole binary mergers observed with gravitational waves is an exciting milestone as it relates to both general relativistic dynamics and astrophysical binary ...formation scenarios. In this study, we revisit the evidence for spin-precession in GW200129 and localize its origin to data in LIGO Livingston in the 20--50\,Hz frequency range where the signal amplitude is lower than expected from a non-precessing binary given all the other data. These data are subject to known data quality issues as a glitch was subtracted from the detector's strain data. The lack of evidence for spin-precession in LIGO Hanford leads to a noticeable inconsistency between the inferred binary mass ratio and precessing spin in the two LIGO detectors, something not expected from solely different Gaussian noise realizations. We revisit the LIGO Livingston glitch mitigation and show that the difference between a spin-precessing and a non-precessing interpretation for GW200129 is smaller than the statistical and systematic uncertainty of the glitch subtraction, finding that the support for spin-precession depends sensitively on the glitch modeling. We also investigate the signal-to-noise ratio \(\sim7\) trigger in the less sensitive Virgo detector. Though not influencing the spin-precession studies, the Virgo trigger is grossly inconsistent with the ones in LIGO Hanford and LIGO Livingston as it points to a much heavier system. We interpret the Virgo data in the context of further data quality issues. While our results do not disprove the presence of spin-precession in GW200129, we argue that any such inference is contingent upon the statistical and systematic uncertainty of the glitch mitigation. Our study highlights the role of data quality investigations when inferring subtle effects such as spin-precession for short signals such as the ones produced by high-mass systems.
The increasing sensitivity of gravitational-wave detectors has brought about an increase in the rate of astrophysical signal detections as well as the rate of "glitches"; transient and non-Gaussian ...detector noise. Temporal overlap of signals and glitches in the detector presents a challenge for inference analyses that typically assume the presence of only Gaussian detector noise. In this study we perform an extensive exploration of the efficacy of a recently proposed method that models the glitch with sine-Gaussian wavelets while simultaneously modeling the signal with compact-binary waveform templates. We explore a wide range of glitch families and signal morphologies and demonstrate that the joint modeling of glitches and signals (with wavelets and templates respectively) can reliably separate the two. We find that the glitches that most affect parameter estimation are also the glitches that are well modeled by such wavelets due to their compact time-frequency signature. As a further test, we investigate the robustness of this analysis against waveform systematics like those arising from the exclusion of higher-order modes and spin-precession effects. Our analysis provides an estimate of the signal parameters; the glitch waveform to be subtracted from the data; and an assessment of whether some detected excess power consists of a glitch, signal, or both. We analyze the low-significance triggers (191225_215715 and 200114_020818) and find that they are both consistent with glitches overlapping high-mass signals.