No Arabic abstract
The global gravitational-wave detector network achieves higher detection rates, better parameter estimates, and more accurate sky localisation, as the number of detectors, $mathcal{I}$ increases. This paper quantifies network performance as a function of $mathcal{I}$ for BayesWave, a source-agnostic, wavelet-based, Bayesian algorithm which distinguishes between true astrophysical signals and instrumental glitches. Detection confidence is quantified using the signal-to-glitch Bayes factor, $mathcal{B}_{mathcal{S},mathcal{G}}$. An analytic scaling is derived for $mathcal{B}_{mathcal{S},mathcal{G}}$ versus $mathcal{I}$, the number of wavelets, and the network signal-to-noise ratio, SNR$_text{net}$, which is confirmed empirically via injections into detector noise of the Hanford-Livingston (HL), Hanford-Livingston-Virgo (HLV), and Hanford-Livingston-KAGRA-Virgo (HLKV) networks at projected sensitivities for the fourth observing run (O4). The empirical and analytic scalings are consistent; $mathcal{B}_{mathcal{S},mathcal{G}}$ increases with $mathcal{I}$. The accuracy of waveform reconstruction is quantified using the overlap between injected and recovered waveform, $mathcal{O}_text{net}$. The HLV and HLKV network recovers $87%$ and $86%$ of the injected waveforms with $mathcal{O}_text{net}>0.8$ respectively, compared to $81%$ with the HL network. The accuracy of BayesWave sky localisation is $approx 10$ times better for the HLV network than the HL network, as measured by the search area, $mathcal{A}$, and the sky areas contained within $50%$ and $90%$ confidence intervals. Marginal improvement in sky localisation is also observed with the addition of KAGRA.
We describe updates and improvements to the BayesWave gravitational wave transient analysis pipeline, and provide examples of how the algorithm is used to analyze data from ground-based gravitational wave detectors. BayesWave models gravitational wave signals in a morphology-independent manner through a sum of frame functions, such as Morlet-Gabor wavelets or chirplets. BayesWave models the instrument noise using a combination of a parametrized Gaussian noise component and non-stationary and non-Gaussian noise transients. Both the signal model and noise model employ trans-dimensional sampling, with the complexity of the model adapting to the requirements of the data. The flexibility of the algorithm makes it suitable for a variety of analyses, including reconstructing generic unmodeled signals; cross checks against modeled analyses for compact binaries; as well as separating coherent signals from incoherent instrumental noise transients (glitches). The BayesWave model has been extended to account for gravitational wave signals with generic polarization content and the simultaneous presence of signals and glitches in the data. We describe updates in the BayesWave prior distributions, sampling proposals, and burn-in stage that provide significantly improved sampling efficiency. We present standard review checks indicating the robustness and convergence of the BayesWave trans-dimensional sampler.
We provide a comprehensive multi-aspect study on the performance of a pipeline used by the LIGO-Virgo Collaboration for estimating parameters of gravitational-wave bursts. We add simulated signals with four different morphologies (sine-Gaussians, Gaussians, white-noise bursts, and binary black hole signals) to simulated noise samples representing noise of the two Advanced LIGO detectors during their first observing run. We recover them with the BayesWave (BW) pipeline to study its accuracy in sky localization, waveform reconstruction, and estimation of model-independent waveform parameters. BW localizes sources with a level of accuracy comparable for all four morphologies, with the median separation of actual and estimated sky locations ranging from 25.1$^{circ}$ to 30.3$^{circ}$. This is a reasonable accuracy in the two-detector case, and is comparable to accuracies of other localization methods studied previously. As BW reconstructs generic transient signals with sine-Gaussian wavelets, it is unsurprising that BW performs the best in reconstructing sine-Gaussian and Gaussian waveforms. BWs accuracy in waveform reconstruction increases steeply with network signal-to-noise ratio (SNR$_{rm net}$), reaching a $85%$ and $95%$ match between the reconstructed and actual waveform below SNR$_{rm net} approx 20$ and SNR$_{rm net} approx 50$, respectively, for all morphologies. BWs accuracy in estimating central moments of waveforms is only limited by statistical errors in the frequency domain, and is affected by systematic errors too in the time domain as BW cannot reconstruct low-amplitude parts of signals overwhelmed by noise. The figures of merit we introduce can be used in future characterizations of parameter estimation pipelines.
Autonomous gravitational-wave searches -- fully automated analyses of data that run without human intervention or assistance -- are desirable for a number of reasons. They are necessary for the rapid identification of gravitational-wave burst candidates, which in turn will allow for follow-up observations by other observatories and the maximum exploitation of their scientific potential. A fully automated analysis would also circumvent the traditional by hand setup and tuning of burst searches that is both labourious and time consuming. We demonstrate a fully automated search with X-Pipeline, a software package for the coherent analysis of data from networks of interferometers for detecting bursts associated with GRBs and other astrophysical triggers. We discuss the methods X-Pipeline uses for automated running, including background estimation, efficiency studies, unbiased optimal tuning of search thresholds, and prediction of upper limits. These are all done automatically via Monte Carlo with multiple independent data samples, and without requiring human intervention. As a demonstration of the power of this approach, we apply X-Pipeline to LIGO data to search for gravitational-wave emission associated with GRB 031108. We find that X-Pipeline is sensitive to signals approximately a factor of 2 weaker in amplitude than those detectable by the cross-correlation technique used in LIGO searches to date. We conclude with the prospects for running X-Pipeline as a fully autonomous, near real-time triggered burst search in the next LSC-Virgo Science Run.
A central challenge in Gravitational Wave Astronomy is identifying weak signals in the presence of non-stationary and non-Gaussian noise. The separation of gravitational wave signals from noise requires good models for both. When accurate signal models are available, such as for binary Neutron star systems, it is possible to make robust detection statements even when the noise is poorly understood. In contrast, searches for un-modeled transient signals are strongly impacted by the methods used to characterize the noise. Here we take a Bayesian approach and introduce a multi-component, variable dimension, parameterized noise model that explicitly accounts for non-stationarity and non-Gaussianity in data from interferometric gravitational wave detectors. Instrumental transients (glitches) and burst sources of gravitational waves are modeled using a Morlet-Gabor continuous wavelet frame. The number and placement of the wavelets is determined by a trans-dimensional Reversible Jump Markov Chain Monte Carlo algorithm. The Gaussian component of the noise and sharp line features in the noise spectrum are modeled using the BayesLine algorithm, which operates in concert with the wavelet model.
The multi-band template analysis (MBTA) pipeline is a low-latency coincident analysis pipeline for the detection of gravitational waves (GWs) from compact binary coalescences. MBTA runs with a low computational cost, and can identify candidate GW events online with a sub-minute latency. The low computational running cost of MBTA also makes it useful for data quality studies. Events detected by MBTA online can be used to alert astronomical partners for electromagnetic follow-up. We outline the current status of MBTA and give details of recent pipeline upgrades and validation tests that were performed in preparation for the first advanced detector observing period. The MBTA pipeline is ready for the outset of the advanced detector era and the exciting prospects it will bring.