Source author record

Reinhard Prix

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

21 published item(s)

preprint2021arXiv

Results from an Einstein@Home search for continuous gravitational waves from G347.3 at low frequencies in LIGO O2 data

We present results of a search for periodic gravitational wave signals with frequency between 20 and 400 Hz, from the neutron star in the supernova remnant G347.3-0.5, using LIGO O2 public data. The search is deployed on the volunteer computing project Einstein@Home, with thousands of participants donating compute cycles to make this endevour possible. We find no significant signal candidate and set the most constraining upper limits to date on the amplitude of gravitational wave signals from the target, corresponding to deformations below $10^{-6}$ in a large part of the band. At the frequency of best strain sensitivity, near $166$ Hz, we set 90\%\ confidence upper limits on the gravitational wave intrinsic amplitude of $h_0^{90\%}\approx 7.0\times10^{-26}$. Over most of the frequency range our upper limits are a factor of 20 smaller than the indirect age-based upper limit.

preprint2020arXiv

Deep-Learning Continuous Gravitational Waves: Multiple detectors and realistic noise

The sensitivity of wide-parameter-space searches for continuous gravitational waves is limited by computational cost. Recently it was shown that Deep Neural Networks (DNNs) can perform all-sky searches directly on (single-detector) strain data, potentially providing a low-computing-cost search method that could lead to a better overall sensitivity. Here we expand on this study in two respects: (i) using (simulated) strain data from two detectors simultaneously, and (ii) training for directed (i.e.\ single sky-position) searches in addition to all-sky searches. For a data timespan of $T = 10^5\, s$, the all-sky two-detector DNN is about $7\%$ less sensitive (in amplitude $h_0$) at low frequency ($f=20\,Hz$), and about $51\,\%$ less sensitive at high frequency ($f=1000\,Hz$) compared to fully-coherent matched-filtering (using WEAVE). In the directed case the sensitivity gap compared to matched-filtering ranges from about $7-14\%$ at $f=20\,Hz$ to about $37-49\%$ at $f=1500\,Hz$. Furthermore we assess the DNN's ability to generalize in signal frequency, spindown and sky-position, and we test its robustness to realistic data conditions, namely gaps in the data and using real LIGO detector noise. We find that the DNN performance is not adversely affected by gaps in the test data or by using a relatively undisturbed band of LIGO detector data instead of Gaussian noise. However, when using a more disturbed LIGO band for the tests, the DNN's detection performance is substantially degraded due to the increase in false alarms, as expected.

preprint2020arXiv

Search for Continuous Gravitational Waves from the Central Compact Objects in Supernova Remnants Cassiopeia A, Vela Jr. and G347.3-0.5

We perform a sub-threshold follow-up search for continuous nearly-monochromatic gravitational waves from the central compact objects associated with the supernova remnants Vela Jr., Cassiopeia A, and SNR G347.3$-$0.5. Across the three targets, we investigate the most promising ~ 10,000 combinations of gravitational wave frequency and frequency derivative values, based on the results from an Einstein@Home search of the LIGO O1 observing run data, dedicated to these objects. The selection threshold is set so that a signal could be confirmed using the newly released O2 run LIGO data. In order to achieve best sensitivity we perform two separate follow-up searches, on two distinct stretches of the O2 data. Only one candidate survives the first O2 follow-up investigation, associated with the central compact object in SNR G347.3-0.5, but it is not conclusively confirmed. In order to assess a possible astrophysical origin we use archival X-ray observations and search for amplitude modulations of a pulsed signal at the putative rotation frequency of the neutron star and its harmonics. This is the first extensive electromagnetic follow-up of a continuous gravitational wave candidate performed to date. No significant associated signal is identified. New X-ray observations contemporaneous with the LIGO O3 run will enable a more sensitive search for an electromagnetic counterpart. A focused gravitational wave search in O3 data based on the parameters provided here should be easily able to shed light on the nature of this outlier. Noise investigations on the LIGO instruments could also reveal the presence of a coherent contamination.

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}$.

preprint2015arXiv

Directed searches for continuous gravitational waves from binary systems: parameter-space metrics and optimal Scorpius X-1 sensitivity

We derive simple analytic expressions for the (coherent and semi-coherent) phase metrics of continuous-wave sources in low-eccentricity binary systems, both for the long-segment and short- segment regimes (compared to the orbital period). The resulting expressions correct and extend previous results found in the literature. We present results of extensive Monte-Carlo studies comparing metric mismatch predictions against the measured loss of detection statistic for binary parameter offsets. The agreement is generally found to be within ~ 10%-30%. As an application of the metric template expressions, we estimate the optimal achievable sensitivity of an Einstein@Home directed search for Scorpius X-1, under the assumption of sufficiently small spin wandering. We find that such a search, using data from the upcoming advanced detectors, would be able to beat the torque- balance level [1,2] up to a frequency of ~ 500 - 600 Hz, if orbital eccentricity is well-constrained, and up to a frequency of ~ 160 - 200 Hz for more conservative assumptions about the uncertainty on orbital eccentricity.

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

Postprocessing methods used in the search for continuous gravitational-wave signals from the Galactic Center

The search for continuous gravitational-wave signals requires the development of techniques that can effectively explore the low-significance regions of the candidate set. In this paper we present the methods that were developed for a search for continuous gravitational-wave signals from the Galactic Center [1]. First, we present a data-selection method that increases the sensitivity of the chosen data set by 20%-30% compared to the selection methods used in previous directed searches. Second, we introduce postprocessing methods that reliably rule out candidates that stem from random fluctuations or disturbances in the data. In the context of [J. Aasi et al., Phys. Rev. D 88, 102002 (2013)] their use enabled the investigation of candidates three standard deviations below the loudest expected candidate in Gaussian noise from the entire search. Such low-significance regions had not been explored in continuous gravitational-wave searches before. We finally present a new procedure for deriving upper limits on the gravitational-wave amplitude, which is several times faster with respect to the standard injection-and-search approach commonly used.

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.

preprint2014arXiv

Fully coherent follow-up of continuous gravitational-wave candidates: an application to Einstein@Home results

We characterize and present the details of the follow-up method used on the most significant outliers of the Hough Einstein@Home all-sky search for continuous gravitational waves arXiv:1207.7176. This follow-up method is based on the two-stage approach introduced in arXiv:1303.2471, consisting of a semicoherent refinement followed by a fully coherent zoom. We quantify the efficiency of the follow-up pipeline using simulated signals in Gaussian noise. This pipeline does not search beyond first-order frequency spindown, and therefore we also evaluate its robustness against second-order spindown. We present the details of the Hough Einstein@Home follow-up (arXiv:1207.7176) on three hardware-injected signals and on the 8 most significant outliers of unknown origin.

preprint2014arXiv

New Coordinates for the Amplitude Parameter Space of Continuous Gravitational Waves

The parameter space for continuous gravitational waves (GWs) can be divided into amplitude parameters (signal amplitude, inclination and polarization angles describing the orientation of the source, and an initial phase) and phase-evolution parameters. The division is useful in part because the Jaranowski-Krolak-Schutz (JKS) coordinates on the four-dimensional amplitude parameter space allow the GW signal to be written as a linear combination of four template waveforms with the JKS coordinates as coefficients. We define a new set of coordinates on the amplitude parameter space, with the same properties, which is more closely connected to the physical amplitude parameters. These naturally divide into two pairs of Cartesian-like coordinates on two-dimensional subspaces, one corresponding to left- and the other to right-circular polarization. We thus refer to these as CPF (circular polarization factored) coordinates. The corresponding two sets of polar coordinates (known as CPF-polar) can be related in a simple way to the physical parameters. We illustrate some simplifying applications for these various coordinate systems, such as: a calculation of Jacobians between various coordinate systems; an illustration of the signal coordinate singularities associated with left- and right-circular polarization, which correspond to the origins of the two two-dimensional subspaces; and an elucidation of the form of the log-likelihood ratio between hypotheses of Gaussian noise with and without a continuous GW signal. These are used to illustrate some of the prospects for approximate evaluation of a Bayesian detection statistic defined by marginalization over the physical parameter space. Additionally, in the presence of simplifying assumptions about the observing geometry, we are able to explicitly evaluate the integral for the Bayesian detection statistic, and compare it to the approximate results.

preprint2013arXiv

Flat parameter-space metric for all-sky searches for gravitational-wave pulsars

All-sky, broadband, coherent searches for gravitational-wave pulsars are computationally limited. It is therefore important to make efficient use of available computational resources, notably by minimizing the number of templates used to cover the signal parameter space of sky position and frequency evolution. For searches over the sky, however, the required template density (determined by the parameter-space metric) is different at each sky position, which makes it difficult in practice to achieve an efficient covering. Previous work on this problem has found various choices of sky and frequency coordinates that render the parameter-space metric approximately constant, but which are limited to coherent integration times of either less than a few days, or greater than several months. These limitations restrict the sensitivity achievable by hierarchical all-sky searches, and hinder the development of follow-up pipelines for interesting gravitational-wave pulsar candidates. We present a new flat parameter-space metric approximation, and associated sky and frequency coordinates, that do not suffer from these limitations. Furthermore, the new metric is numerically well-conditioned, which facilitates its practical use.

preprint2013arXiv

Fully coherent follow-up of continuous gravitational-wave candidates

The search for continuous gravitational waves from unknown isolated sources is computationally limited due to the enormous parameter space that needs to be covered and the weakness of the expected signals. Therefore semi-coherent search strategies have been developed and applied in distributed computing environments such as Einstein@Home, in order to narrow down the parameter space and identify interesting candidates. However, in order to optimally confirm or dismiss a candidate as a possible gravitational-wave signal, a fully-coherent follow-up using all the available data is required. We present a general method and implementation of a direct (2-stage) transition to a fully-coherent follow-up on semi-coherent candidates. This method is based on a grid-less Mesh Adaptive Direct Search (MADS) algorithm using the F-statistic. We demonstrate the detection power and computing cost of this follow-up procedure using extensive Monte-Carlo simulations on (simulated) semi-coherent candidates from a directed as well as from an all-sky search setup.

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.

preprint2012arXiv

Search for Continuous Gravitational Waves: Optimal StackSlide method at fixed computing cost

Coherent wide parameter-space searches for continuous gravitational waves are typically limited in sensitivity by their prohibitive computing cost. Therefore semi-coherent methods (such as StackSlide) can often achieve a better sensitivity. We develop an analytical method for finding optimal StackSlide parameters at fixed computing cost under ideal conditions of gapless data with Gaussian stationary noise. This solution separates two regimes: an unbounded regime, where it is always optimal to use all the data, and a bounded regime with a finite optimal observation time. Our analysis of the sensitivity scaling reveals that both the fine- and coarse-grid mismatches contribute equally to the average StackSlide mismatch, an effect that had been overlooked in previous studies. We discuss various practical examples for the application of this optimization framework, illustrating the potential gains in sensitivity compared to previous searches.

preprint2011arXiv

Bayesian versus frequentist upper limits

While gravitational waves have not yet been measured directly, data analysis from detection experiments commonly includes an upper limit statement. Such upper limits may be derived via a frequentist or Bayesian approach; the theoretical implications are very different, and on the technical side, one notable difference is that one case requires maximization of the likelihood function over parameter space, while the other requires integration. Using a simple example (detection of a sinusoidal signal in white Gaussian noise), we investigate the differences in performance and interpretation, and the effect of the "trials factor", or "look-elsewhere effect".

preprint2011arXiv

Search method for long-duration gravitational-wave transients from neutron stars

We introduce a search method for a new class of gravitational-wave signals, namely long-duration O(hours - weeks) transients from spinning neutron stars. We discuss the astrophysical motivation from glitch relaxation models and we derive a rough estimate for the maximal expected signal strength based on the superfluid excess rotational energy. The transient signal model considered here extends the traditional class of infinite-duration continuous-wave signals by a finite start-time and duration. We derive a multi-detector Bayes factor for these signals in Gaussian noise using $\F$-statistic amplitude priors, which simplifies the detection statistic and allows for an efficient implementation. We consider both a fully coherent statistic, which is computationally limited to directed searches for known pulsars, and a cheaper semi-coherent variant, suitable for wide parameter-space searches for transients from unknown neutron stars. We have tested our method by Monte-Carlo simulation, and we find that it outperforms orthodox maximum-likelihood approaches both in sensitivity and in parameter-estimation quality.

preprint2010arXiv

Searching for Galactic White Dwarf Binaries in Mock LISA Data using an F-Statistic Template Bank

We describe an F-statistic search for continuous gravitational waves from galactic white-dwarf binaries in simulated LISA Data. Our search method employs a hierarchical template-grid based exploration of the parameter space. In the first stage, candidate sources are identified in searches using different simulated laser signal combinations (known as TDI variables). Since each source generates a primary maximum near its true "Doppler parameters" (intrinsic frequency and sky position) as well as numerous secondary maxima of the F-statistic in Doppler parameter space, a search for multiple sources needs to distinguish between true signals and secondary maxima associated with other, "louder" signals. Our method does this by applying a coincidence test to reject candidates which are not found at nearby parameter space positions in searches using each of the three TDI variables. For signals surviving the coincidence test, we perform a fully coherent search over a refined parameter grid to provide an accurate parameter estimation for the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters with errors comparable to statistical expectations, although there is still some systematic excess with respect to statistical errors expected from Gaussian noise. An experimental iterative pipeline with seven rounds of signal subtraction and re-analysis of the residuals allows us to increase the number of signals recovered to a total of 3419 with 29 false alarms.

preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

preprint2009arXiv

Targeted search for continuous gravitational waves: Bayesian versus maximum-likelihood statistics

We investigate the Bayesian framework for detection of continuous gravitational waves (GWs) in the context of targeted searches, where the phase evolution of the GW signal is assumed to be known, while the four amplitude parameters are unknown. We show that the orthodox maximum-likelihood statistic (known as F-statistic) can be rediscovered as a Bayes factor with an unphysical prior in amplitude parameter space. We introduce an alternative detection statistic ("B-statistic") using the Bayes factor with a more natural amplitude prior, namely an isotropic probability distribution for the orientation of GW sources. Monte-Carlo simulations of targeted searches show that the resulting Bayesian B-statistic is more powerful in the Neyman-Pearson sense (i.e. has a higher expected detection probability at equal false-alarm probability) than the frequentist F-statistic.

preprint2007arXiv

Search for Continuous Gravitational Waves: simple criterion for optimal detector networks

We derive a simple algebraic criterion to select the optimal detector network for a coherent wide parameter-space (all-sky) search for continuous gravitational waves. Optimality in this context is defined as providing the highest (average) sensitivity per computing cost. This criterion is a direct consequence of the properties of the multi-detector F-statistic metric, which has been derived recently. Interestingly, the choice of the optimal network only depends on the noise-levels and duty-cycles of the respective detectors, and not on the available computing power.

preprint2006arXiv

The search for continuous gravitational waves: metric of the multi-detector F-statistic

We develop a general formalism for the parameter-space metric of the multi-detector F-statistic, which is a matched-filtering detection statistic for continuous gravitational waves. We find that there exists a whole family of F-statistic metrics, parametrized by the (unknown) amplitude parameters of the gravitational wave. The multi-detector metric is shown to be expressible in terms of noise-weighted averages of single-detector contributions, which implies that the number of templates required to cover the parameter space does not scale with the number of detectors. Contrary to using a longer observation time, combining detectors of similar sensitivity is therefore the computationally cheapest way to improve the sensitivity of coherent wide-parameter searches for continuous gravitational waves. We explicitly compute the F-statistic metric family for signals from isolated spinning neutron stars, and we numerically evaluate the quality of different metric approximations in a Monte-Carlo study. The metric predictions are tested against the measured mismatches and we identify regimes in which the local metric is no longer a good description of the parameter-space structure.