Source author record

David Keitel

David Keitel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

gr-qc astro-ph.HE astro-ph.IM astro-ph.CO math.ST physics.data-an Statistics Theory

Catalog footprint

What is connected

13works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Empirically estimating the distribution of the loudest candidate from a gravitational-wave search

Searches for gravitational-wave signals are often based on maximizing a detection statistic over a bank of waveform templates, covering a given parameter space with a variable level of correlation. Results are often evaluated using a noise-hypothesis test, where the background is characterized by the sampling distribution of the loudest template. In the context of continuous gravitational-wave searches, properly describing said distribution is an open problem: current approaches focus on a particular detection statistic and neglect template-bank correlations. We introduce a new approach using extreme value theory to describe the distribution of the loudest template's detection statistic in an arbitrary template bank. Our new proposal automatically generalizes to a wider class of detection statistics, including (but not limited to) line-robust statistics and transient continuous-wave signal hypotheses, and improves the estimation of the expected maximum detection statistic at a negligible computing cost. The performance of our proposal is demonstrated on simulated data as well as by applying it to different kinds of (transient) continuous-wave searches using O2 Advanced LIGO data. We release an accompanying Python software package, distromax, implementing our new developments.

preprint2022arXiv

Search setup for long-duration transient gravitational waves from glitching pulsars during LIGO-Virgo third observing run

Pulsars are spinning neutron stars which emit an electromagnetic beam. We expect pulsars to slowly decrease their rotational frequency. However, sudden increases of the rotational frequency have been observed from different pulsars. These events are called "glitches" and they are followed by a relaxation phase with timescales from days to months. Gravitational wave (GW) emission may follow these peculiar events. We give an overview of the setup for an analysis of GW data from the Advanced LIGO and Virgo third observing run (O3) arXiv:2112.10990 searching for transient GW signals lasting hours to months after glitches in known pulsars. The search method consists of placing a template grid in frequency-spindown space with fixed grid spacings. Then, for each point we compute the transient F-statistic which is maximized over a set of transient parameters like the duration and start time of the potential signals. A threshold on the detection statistic is then set, and we search for peaks over the parameter space for each candidate.

preprint2021arXiv

A detailed analysis of GW190521 with phenomenological waveform models

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF.

preprint2021arXiv

New twists in compact binary waveform modelling: a fast time domain model for precession

We present IMRPhenomTPHM, a phenomenological model for the gravitational wave signals emitted by the coalescence of quasi-circular precessing binary black holes systems. The model is based on the "twisting up" approximation, which maps non-precessing signals to precessing ones in terms of a time dependent rotation described by three Euler angles, and which has been utilized in several frequency domain waveform models that have become standard tools in gravitational wave data analysis. Our model is however constructed in the time domain, which allows several improvements over the frequency domain models: we do not use the stationary phase approximation, we employ a simple approximation for the precessing Euler angles for the ringdown signal, and we implement a new method for computing the Euler angles through the evolution of the spin dynamics of the system, which is more accurate and also computationally efficient.

preprint2020arXiv

Enhancing Gravitational-Wave Science with Machine Learning

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

preprint2020arXiv

Towards the routine use of subdominant harmonics in gravitational-wave inference: re-analysis of GW190412 with generation X waveform models

We re-analyse the gravitational-wave event GW190412 with state-of-the-art phenomenological waveform models. This event, which has been associated with a black hole merger, is interesting due to the significant contribution from subdominant harmonics. We use both frequency-domain and time-domain waveform models. The PhenomX waveform models constitute the fourth generation of frequency-domain phenomenological waveforms for black hole binary coalescence; they have more recently been complemented by the time-domain PhenomT models, which open up new strategies to model precession and eccentricity, and to perform tests of general relativity with the phenomenological waveforms approach. Both PhenomX and PhenomT have been constructed with similar techniques and accuracy goals, and due to their computational efficiency this "generation X" model family allows the routine use of subdominant spherical harmonics in Bayesian inference. We show the good agreement between these and other state-of-the-art waveform models for GW190412, and discuss the improvements over the previous generation of phenomenological waveform models. We also discuss practical aspects of Bayesian inference such as run convergence, variations of sampling parameters, and computational cost.

preprint2016arXiv

An Einstein@home search for continuous gravitational waves from Cassiopeia A

We report the results of a directed search for continuous gravitational-wave emission in a broad frequency range (between 50 and 1000 Hz) from the central compact object of the supernova remnant Cassiopeia A (Cas A). The data comes from the sixth science run of LIGO and the search is performed on the volunteer distributed computing network Einstein@Home. We find no significant signal candidate, and set the most constraining upper limits to date on the gravitational-wave emission from Cas A, which beat the indirect age-based upper limit across the entire search range. At around 170 Hz (the most sensitive frequency range), we set 90% confidence upper limits on the gravitational wave amplitude $h_0$ of $\sim\!\!~2.9\times 10^{-25}$, roughly twice as constraining as the upper limits from previous searches on Cas A. The upper limits can also be expressed as constraints on the ellipticity of Cas A; with a few reasonable assumptions, we show that at gravitational-wave frequencies greater than 300~Hz, we can exclude an ellipticity of $\gtrsim\!\!~10^{-5}$.

preprint2016arXiv

Robust semicoherent searches for continuous gravitational waves with noise and signal models including hours to days long transients

The vulnerability to single-detector instrumental artifacts in standard detection methods for long-duration quasimonochromatic gravitational waves from nonaxisymmetric rotating neutron stars [continuous waves (CWs)] was addressed in past work [D. Keitel et al., Phys. Rev. D 89, 064023 (2014).] by a Bayesian approach. An explicit model of persistent single-detector disturbances led to a generalized detection statistic with improved robustness against such artifacts. Since many strong outliers in semicoherent searches of LIGO data are caused by transient disturbances that last only a few hours, we extend the noise model to cover such limited-duration disturbances, and demonstrate increased robustness in realistic simulated data. Besides long-duration CWs, neutron stars could also emit transient signals which, for a limited time, also follow the CW signal model (tCWs). As a pragmatic alternative to specialized transient searches, we demonstrate how to make standard semicoherent CW searches more sensitive to transient signals. Considering tCWs in a single segment of a semicoherent search, Bayesian model selection yields a new detection statistic that does not add significant computational cost. On simulated data, we find that it increases sensitivity towards tCWs, even of varying durations, while not sacrificing sensitivity to classical CW signals, and still being robust to transient or persistent single-detector instrumental artifacts.

preprint2015arXiv

Line-robust statistics for continuous gravitational waves: safety in the case of unequal detector sensitivities

The multi-detector F-statistic is close to optimal for detecting continuous gravitational waves (CWs) in Gaussian noise. However, it is susceptible to false alarms from instrumental artefacts, for example quasi-monochromatic disturbances ('lines'), which resemble a CW signal more than Gaussian noise. In a recent paper [Keitel et al 2014, PRD 89 064023], a Bayesian model selection approach was used to derive line-robust detection statistics for CW signals, generalising both the F-statistic and the F-statistic consistency veto technique and yielding improved performance in line-affected data. Here we investigate a generalisation of the assumptions made in that paper: if a CW analysis uses data from two or more detectors with very different sensitivities, the line-robust statistics could be less effective. We investigate the boundaries within which they are still safe to use, in comparison with the F-statistic. Tests using synthetic draws show that the optimally-tuned version of the original line-robust statistic remains safe in most cases of practical interest. We also explore a simple idea on further improving the detection power and safety of these statistics, which we however find to be of limited practical use.

preprint2015arXiv

Search for continuous gravitational waves: Improving robustness versus instrumental artifacts

The standard multidetector F-statistic for continuous gravitational waves is susceptible to false alarms from instrumental artifacts, for example monochromatic sinusoidal disturbances ('lines'). This vulnerability to line artifacts arises because the F-statistic compares the signal hypothesis to a Gaussian-noise hypothesis, and hence is triggered by anything that resembles the signal hypothesis more than Gaussian noise. Various ad-hoc veto methods to deal with such line artifacts have been proposed and used in the past. Here we develop a Bayesian framework that includes an explicit alternative hypothesis to model disturbed data. We introduce a simple line model that defines lines as signal candidates appearing only in one detector. This allows us to explicitly compute the odds between the signal hypothesis and an extended noise hypothesis, resulting in a new detection statistic that is more robust to instrumental artifacts. We present and discuss results from Monte-Carlo tests on both simulated data and on detector data from the fifth LIGO science run. We find that the line-robust statistic retains the detection power of the standard F-statistic in Gaussian noise. In the presence of line artifacts it is more sensitive, even compared to the popular F-statistic consistency veto, over which it improves by as much as a factor of two in detectable signal strength.

preprint2013arXiv

Octahedron configuration for a displacement noise-cancelling gravitational wave detector in space

We study for the first time a three-dimensional octahedron constellation for a space-based gravitational wave detector, which we call the Octahedral Gravitational Observatory (OGO). With six spacecraft the constellation is able to remove laser frequency noise and acceleration disturbances from the gravitational wave signal without needing LISA-like drag-free control, thereby simplifying the payloads and placing less stringent demands on the thrusters. We generalize LISA's time-delay interferometry to displacement-noise free interferometry (DFI) by deriving a set of generators for those combinations of the data streams that cancel laser and acceleration noise. However, the three-dimensional configuration makes orbit selection complicated. So far, only a halo orbit near the Lagrangian point L1 has been found to be stable enough, and this allows only short arms up to 1400 km. We derive the sensitivity curve of OGO with this arm length, resulting in a peak sensitivity of about $2\times10^{-23}\,\mathrm{Hz}^{-1/2}$ near 100 Hz. We compare this version of OGO to the present generation of ground-based detectors and to some future detectors. We also investigate the scientific potentials of such a detector, which include observing gravitational waves from compact binary coalescences, the stochastic background and pulsars as well as the possibility to test alternative theories of gravity. We find a mediocre performance level for this short-arm-length detector, between those of initial and advanced ground-based detectors. Thus, actually building a space-based detector of this specific configuration does not seem very efficient. However, when alternative orbits that allow for longer detector arms can be found, a detector with much improved science output could be constructed using the octahedron configuration and DFI solutions demonstrated in this paper. (abridged)

preprint2012arXiv

An F-statistic based multi-detector veto for detector artifacts in continuous-wave gravitational wave data

Continuous gravitational waves (CW) are expected from spinning neutron stars with non-axisymmetric deformations. A network of interferometric detectors (LIGO, Virgo and GEO600) is looking for these signals. They are predicted to be very weak and retrievable only by integration over long observation times. One of the standard methods of CW data analysis is the multi-detector F-statistic. In a typical search, the F-statistic is computed over a range in frequency, spin-down and sky position, and the candidates with highest F values are kept for further analysis. However, this detection statistic is susceptible to a class of noise artifacts, strong monochromatic lines in a single detector. By assuming an extended noise model - standard Gaussian noise plus single-detector lines - we can use a Bayesian odds ratio to derive a generalized detection statistic, the line veto (LV-) statistic. In the absence of lines, it behaves similarly to the F-statistic, but it is more robust against line artifacts. In the past, ad-hoc post-processing vetoes have been implemented in searches to remove these artifacts. Here we provide a systematic framework to develop and benchmark this class of vetoes. We present our results from testing this LV-statistic on simulated data.

preprint2011arXiv

Constrained probability distributions of correlation functions

Context: Two-point correlation functions are used throughout cosmology as a measure for the statistics of random fields. When used in Bayesian parameter estimation, their likelihood function is usually replaced by a Gaussian approximation. However, this has been shown to be insufficient. Aims: For the case of Gaussian random fields, we search for an exact probability distribution of correlation functions, which could improve the accuracy of future data analyses. Methods: We use a fully analytic approach, first expanding the random field in its Fourier modes, and then calculating the characteristic function. Finally, we derive the probability distribution function using integration by residues. We use a numerical implementation of the full analytic formula to discuss the behaviour of this function. Results: We derive the univariate and bivariate probability distribution function of the correlation functions of a Gaussian random field, and outline how higher joint distributions could be calculated. We give the results in the form of mode expansions, but in one special case we also find a closed-form expression. We calculate the moments of the distribution and, in the univariate case, we discuss the Edgeworth expansion approximation. We also comment on the difficulties in a fast and exact numerical implementation of our results, and on possible future applications.